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ABSTRACT 


This  report  presents  time -domain  maximum  likelihood  adaptive  beam- 
forming results  for  data  from  the  Alaska  Long-Period  Array  and  the  Korean 
Short-Period  Array.  Adaptive-beamforming  detection  gain  relative  to  beam- 
steering is  investigated  for  unmixed  long-period  seismic  events  in  the  pre- 
sence of  background  noise  and  for  simulated  mixed  events  where  two  data  sam- 
ples, each  containing  a signal,  are  summed  to  create  a composite  sample  con- 
taining an  interfering  event. 

In  the  adaptive  beamforming  studies  for  background  noise,  the  perfor- 
mance of  two  different  adaptive  algorithms,  the  effect  of  using  closely-spaced 
partial  arrays  instead  of  the  full  ALP  A array,  and  the  effect  of  different  fre- 
quency filters  applied  to  the  input  channels  before  array  beamforming  are  ex- 
amined. 

In  the  interfering-event  analysis,  the  relative  on-azimuth  to  off- 
azimuth  event  strength  at  which  an  on-azimuth  event  detection  is  possible  de- 
termines the  detection  gain  of  adaptive  multichannel  filtering  over  time-shift- 
and-sum  beamforming. 
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A.  PURPOSE  OF  THIS  STUDY 
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The  adaptive  processing  task  of  the  Vela  Network  Evaluation  and  Auto- 
matic Processing  Research  Program  has  as  its  objectives 

• Continued  improvement  of  adaptive-processing  gains  relative 
to  beamsteering  for  unmixed  long -period  seismic  events  in  the 
presence  of  background  noise. 

• Evaluation  of  potential  adaptive-beamforming  detection  improve- 
ment over  beamsteering  for  both  long-period  and  short-period 
signals  buried  in  off-azimuth  interfering  events. 

This  report  deals  with  results  obtained  from  operating  a maximum-likelihood 
adaptive  beamforming  system  on  ALPA  long-period  data  and  Korean  short- 
period  data.  To  synthesize  interfering  events  by  adding  and  scaling  two  dif- 
ferent recorded  data  samples,  a new  computer  program  employing  floating- 
point arithmetic  was  developed  to  achieve  the  objectives  of  this  study.  This 
program  has  the  capability  of  processing  long-period  data  from  NORSAR  and 
JLASA  as  well  as  ALPA  data  and  Korean  short-period  data. 

B.  DESCRIPTION  OF  ARRAYS  USED 

Figure  1-1  diagrams  the  geometrical  layout  of  the  Alaska  Long-Period 
Array  (ALPA),  which  is  a 19-element  hexagonal  long -period  array  with  20- 
kilometer  spacing  between  sites.  Table  1-1  gives  the  ALPA  site  locations.  In 
this  report,  the  sites  at  ALPA  are  referenced  in  terms  of  their  transmission 
order.  Thus,  site  10  refers  to  site  3-45  in  the  official  nomenclature. 
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ADAPTIVE  MULTICHANNEL.  FILTERING 

Multichannel  filtering  is  a form  of  array  processing  in  which  multiple 
channel  inputs  undergo  individual  frequency- shaping  and  phase-shift  filtering 
prior  to  the  channel-summation  operation  which  produces  the  beamformer 
output.  Figure  1-2  is  a schematic  diagram  of  multichannel  filtering.  This 
illustration  incorporates  the  option  to  preprocess  the  transducer  outputs  be- 
fore they  are  input  to  the  multichannel  beamformer.  Examples  of  preproces- 
sing are  frequency  filtering  (most  commonly  with  identical  frequency  respon- 
ses on  all  channels)  and  time  shifting  to  align  waves  emanating  from  a parti- 
cular direction.  The  preprocessed  transducer  outputs  become  the  input  chan- 
nels to  a multichannel  filter  set,  where  individual  filters  (generally  different 
from  channel  to  channel)  are  applied  to  the  input  channels.  These  filters  are 
implemented  as  convolution  filters  in  time-domain  processing  or  as  complex- 
valued multiplicative  filters  in  frequency-domain  processing.  The  multichan- 
nel filter  output  is  created  by  summing  the  individual  filtered  channel  outputs. 

In  systems  where  second-order  statistics  (crosscorrelation  functions 
and  crosspower  spectra)  are  used  to  describe  interrelationships  among  the  in- 
put channels,  there  are  two  basic  forms  of  multichannel  filtering.  In  Wiener- 
Kolmogorov  multichannel  filtering,  the  average  squared  error  between  the  de- 
sired signal  and  the  multichannel  filter  output  is  minimized.  To  minimize  the 
mean  square  error,  the  crosscorrelation  functions  or  crosspower  spectra  be- 
tween the  input  channels  and  the  desired  signal  are  required.  In  maximum 
likelihood  multichannel  filtering,  the  average  squared  output  from  the  multi- 
channel filter  set  is  minimized  subject  to  signal-preservation  constraints  which 
place  some  suitably-chosen  frequency  response  on  the  signal.  For  maximum 
likelihood  multichannel  filtering,  unlike  Wiener-Kolmogorov  filtering,  only  the 
direction  of  the  signal  needs  to  be  specified,  but  not  the  signal-to-noise  ratio. 

Multichannel  filtering  can  be  employed  with  fixed  or  time-varying  filter 
sets.  When  the  filters  are  updated  as  new  data  inputs  enter  the  multichannel 


1-4 


FIGURE  1-2 

SCHEMATIC  DIAGRAM  OF  MULTICHANNEL  FILTERING 


processor,  the  process  is  called  adaptive  filtering.  Adaptive-filtering  algo- 
rithms with  significant  computational  advantages  over  fixed  multichannel  fil- 
tering are  available.  When  the  inputs  to  the  multichannel  processor  are  time- 
stationary (in  the  wide  sense),  these  algorithms  yield  filter  sets  which  con- 
verge in  the  mean  to  the  corresponding  fixed  multichannel  filter  sets.  After 
adaptive  filter  sets  reach  the  vicinity  of  the  corresponding  fixed  filter  sets, 
they  fluctuate  about  the  fixed-filter  solution  in  the  presence  of  time-stationary 
data:  the  adaptive  filters  converge  in  the  mean  in  the  sense  that  the  average 
position  of  the  fluctuating  adaptive  filters  is  identical  to  the  fixed-filter  solu- 
tion. When  the  statistics  of  the  data  entering  the  multichannel  processor  slow- 
ly change  with  time,  adaptive  filtering  can  react  to  the  changes  in  a semi- 
continuous  manner.  If  fixed  filtering  is  used  in  this  situation,  newly-designed 
filters  change  in  a more  abrupt  fashion.  When,  as  in  this  case,  the  statistics 
of  the  data  shift  with  time,  the  adaptive-filter  solution  lags  behind  the  fixed- 
filter  solution  corresponding  to  the  instantaneous  statistics.  The  extent  of  the 
lag  can  be  controlled  by  changing  the  adaptation  rate.  The  choice  of  an  adap- 
tation rate  involves  a tradeoff  between  misadjustment  (higher-than-optimum 
error  or  power  due  to  the  adaptive-filter  fluctuations)  and  the  lag  behind  the 
optimum  instantaneous  fixed-filter  solution.  A different  kind  of  lag  occurs 
when  fixed  filter  sets  are  periodically  redesigned:  statistics  must  be  accumu- 
lated over  a design  interval  so  that,  as  a result,  the  fixed-filter  solution  can- 
not be  implemented  until  the  next  design  interval. 

In  the  conventional  technique  of  array  processing,  simple  time  delays 
or  phase  shifts  are  applied  to  the  input  channels  before  summing  to  generate 
the  beam  output.  Optimum  multichannel  filtering  introduces  considerable  new 
flexibility  into  the  beamforming  process.  Since  it  is  possible  to  weight  the  in- 
put channels  differently,  channels  with  higher  signal-to-noise  ratios  can  be 
emphasized  at  the  expense  of  noisier  channels.  When  well  instrumented  arrays 
are  utilized,  this  capability  is  generally  of  minor  importance.  A far  more 
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consequential  feature  of  adaptive  filtering  is  the  ability  to  form  array  antenna 
patterns  which  optimally  pass  a signal  while  simultaneously  rejecting  propa- 
gating noise.  Deep  nulls  can  be  aimed  toward  off-beam  noise  sources.  When 
strong  off-azimuth  noise  sources  are  present,  the  creation  of  such  nulls  is  an 
automatic  result  of  the  optimality  of  the  multichannel  processor.  The  conven- 
tional time-shift-and-sum  or  phase-shift  processor,  in  contrast,  has  a beam 
pattern  determined  solely  by  the  steer  direction  and  the  array  geometry. 

In  most  cases,  the  potential  improvement  of  optimum  multichannel  fil- 
tering relative  to  beamsteering  is  determined  by  the  coherence  of  the  noise 
field  across  an  array:  the  greater  the  similarity  of  the  noise  field  from  chan- 
nel to  channel,  the  greater  is  the  optimum-multichannel-filter  improvement 
over  beamsteering.  When,  on  the  other  hand,  noise  is  completely  uncorrelated 
between  sensors  and  identical  signal  and  noise  power  levels  are  encountered  at 
all  array  sites,  there  is  no  potential  for  improvement:  in  this  case,  the  opti- 
mum filter  set  is  a beamsteer  processor.  The  decision  to  employ  or  not  to 
employ  an  optimum-filter  technique  of  processing  depends  critically  on  mea- 
surements of  the  noise  field  at  any  given  array.  Once  these  measurements  are 
available,  the  additional  cost  of  implementing  an  optimum -filter  system  can  be 
quantitatively  weighed  against  the  advantages  of  greater  noise  suppression  re- 
lative to  the  conventional  beamsteer  processing  technique. 

D.  DESCRIPTION  OF  THE  MAXIMUM  LIKELIHOOD  ADAPTIVE  FILTER 
ALGORITHM 

The  adaptive -filter  output  y(t)  at  time  t is  formed  by  applying  a con- 
volution filter  to  each  channel  and  summing  the  outputs  of  all  channels: 

M N 

y(t)  = ai(j)xi(t-j) 

i=l  j = -N 
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where  a.(t)  is  the  filter  weight  for  the  i-th  channel  at  a lag  of  j sample 
points,  x.(t-j)  is  the  value  of  the  channel  i at  time  t-j,  M is  the  number 
of  channels,  and  2N+1  is  the  total  length  of  the  filter  in  points.  Prior  to 
forming  the  filter  output,  each  channel  is  time-shifted  to  time-align  energy 
arriving  from  the  desired  steer  direction. 

The  adaptive  filter  weights  are  updated  by  the  following  algorithm: 
new  old 

aj0)  = a.(j)  + X(t)  y(t)  [x(t-j)  - x.(t-j)] 

where 

M 

= if  Y 

i=l 

and  X(t)  controls  the  adaptation  rate  at  time  t . This  update  algorithm  in- 
corporates the  maximum  likelihood  constraints. 

E.  GEOMETRICAL  INTERPRETATION  OF  THE  ABF  ALGORITHM 

In  vector  form,  the  adaptive -beamforming  filter  update  equation  may 
be  written 

A(t  + \t)  - A(t)  = X(t)  XTA(X-X)  , 

where  the  superscript  T denotes  transposition,  and  where  the  filter  weight 
vector  A , the  data  vector  X , and  the  beamsteer  output  vector  X are,  res- 
pectively, 
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The  objective  of  maximum  likelihood  adaptive  beamforming  is  to  re- 
duce the  average  squared  filter  output 

2 T T T T 

y (t)  = (A  X)(X  A)  = A XX  A 

subject  to  a set  of  signal-preservation  constraints  on  the  filter  vector  A . 

After  preshifting  the  input  channels  to  time-align  energy  from  the  look  direc- 
tion, these  constraints  can  be  written 

M 

^ a.(j)  = d(j)  (j  = -N, . . . , -1,  0,  1, ....  N), 

i=l 

where  the  constants  d(j)  specify  a convolution  filter  having  the  desired  fre- 
quency response  on  a signal  from  the  steer  direction.  When  sach  a signal  s(t) 
appears  in  identical  form  on  all  channels,  the  signal  output  from  the  beam- 


d(j)  s(t-j)  . 


N 

3<t.j)  = 

j = -N 


For  the  adaptive  beamforming  employed  in  this  study,  a white  frequency  res 
ponse  is  specified  by  setting 


d(j)  = 8. 

jo 


(j  = -N, . . . , -1,  0,  1 


N)  , 


where 


5. 

jo 


is  the  Kronecker  delta  operator 


a.  .(*  “ j = ° 

JO  ( 0 if  j 4 o . 

2 

To  reduce  the  average  squared  filter  output  y (t)  , the  method  of  steep- 
est descent  (with  two  modifications)  is  used.  In  the  unmodified  form  of  the 
method  of  steepest  descent,  the  filter  vector  A moves  in  the  direction  opposite 
to  the  gradient  of  the  average  squared  filter  output  y2(t)  : 

A(t  + At)  - A(t)  = -cV(ATXXTA)  = -2cXXTA. 


The  first  modification  is  to  replace  the  crosscorrelation  matrix  XX^ 
with  the  rank-one  matrix  XX  formed  from  the  instantaneous  vector  X at 
time  t : 

A(t  + At)  - A(t)  = -2cX(XTA)=  -2  c X y(t)  . 

This  modification  of  Widrow  (Widrow,  1966)  approximates,  in  effect,  the 
crosscorrelation-matrix  time  averaging  through  successive  applications  of  the 
filter  update  algorithm.  The  approximation  becomes  increasingly  accurate  as 
the  rate  of  change  of  the  filter  vector  A is  slowed  by  reducing  the  scalar  con- 
vergence parameter  c . The  filter  vector  A converges  in  the  mean  under 


suitably  prescribed  conditions  (Daniell,  1968)  to  the  vector  obtained  by  using 

7F 


the  crosscorrelation  matrix  XX  in  the  update  equation.  Ultimately  the  filter 
vector  A oscillates  about  its  mean.  The  size  of  the  oscillations  can  be  con- 
trolled by  varying  the  parameter  c . The  reason  for  the  Widrow  modification 
is  a reduction  in  the  computational  operations  required  for  the  filter  update 
from  a number  of  proportional  to  the  square  of  the  dimension  of  the  vectors 
X and  A to  a number  linearly  proportional  to  their  dimension. 


The  second  modification  of  the  steepest-descent  method  is  to  alter  the 
direction  of  the  vector  -cV  y^(t)  so  that  the  ensuing  update  vector 

L T 2 1 

A(t  + At)  - A(t)  is  the  vector  nearest  to  -cV  y (t)  which  satisfies  the  con- 
straints on  the  filter  update  vector.  Since  the  sum  across  channels 


M 


H ai(j) 


i=l 


of  the  filter  vector  A is  a fixed  value  d(j)  at  any  lag  value  j , the  filter  up- 
date vector  A(t  + At)  - A(t)  must  sum  to  zero  at  each  lag: 


M 


z[": 


new 


old 


1 


(j)  - a.(j)  J =0  (j  = -N, . . . , -1 , 0,  1,  „ . . , N) 


i=l 


The  filter  update  vector  must  be  perpendicular  to  each  of  the  (2N+1)  unit 


vectors  U specified  by  their  components 
K 


[H.(j)]k  = yvS  , 


where  8.,  is  the  Kronecker  delta  operator.  Each  of  the  vectors  U,  is  zero 
jk  k 


except  in  the  k-th  lag  position,  where  all  components  are  equal  to  1 / \/ivi  . 
In  vector  form,  the  filter  update  vector  must  satisfy  the  (2N+1)  constraints 


U l [A(t  + At)  - A(t)]  =0  (j  = 0,  1, 


. N) 


I 


1 


I I 


•I 


The  vector  X , since  it  has  identical  components  at  any  lag  value  j , is  a 


linear  combination  of  the  vectors  U and  is  perpendicular  to  all  possible  up- 

K 


date  vectors  satisfying  the  constraint  conditions.  The  vector  (X-X)  , on  the 
other  hand,  satisfies  the  constraint  conditions: 


M N 


Vm  l£<X-X)  = Yj  X ®jk  [xi(j)  - 5(i)] 


i=l  j = -N 


M 


= X [xi<k>  - x<k>] = X xi<k) 

- i=l 


i=l 


M 


Mx(k) 


Mx(k)  - Mx(k)  = 0 . 


Thus  the  vector  X can  be  resolved  into  two  mutually  orthogonal  components 
X (perpendicular  to  the  constraint  space  for  the  filter  update  vector)  and  X-X 
(lying  within  the  constraint  space).  The  negative  -cV^y^(t)J  of  the  scaled 
gradient  of  y (t)  is  a scalar  multiple  - [2cy(t)}x  of  the  vector  X.  The 
nearest  point  to  -cV^y^(t)J  on  the  constraint  space  is  the  vector 
[2cy(t)j  (X-X)  formed  by  subtracting  the  component  -[2cy(t)}x  perpendi- 
cular to  the  contraint  space  from  the  scaled  negative  gradient  vector 
-[2cy(t)]X  . The  final  form  of  the  filter  update  equation  is,  therefore, 


= ^2cy(t)] 


A(t  -I-  At)  - A(t)  = 2 c y (t)  (X-X) 


This  situation  is  illustrated  in  Figure  1-3.  The  resultant  filter  update  vector 
is  the  projection  of  -cV  £y^(t)J  onto  the  constraint  space. 


* 
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F.  A SIMPLE  ILLUSTRATIVE  EXAMPLE 

Suppose  that  three  channels  are  input  to  an  adaptive  beamformer  with 
one  lag  per  channel.  Signals  are  defined  to  be  simultaneous  spikes  on  all 
three  channels.  Noise,  on  the  other  hand,  appears  in  the  form  of  unit- 
amplitude  spikes  on  the  second  channel.  When  signal  and  noise  occur,  they 
are  characterized  by  the  respective  vectors 


x^t)  = s 

rr 

II 

s = 

x2(t)  = s 

and  N = 

x2<*>  = 1 

x3(t)  = s 

_*3<‘>  = °_ 

Figure  1-4  depicts  time  series  inputs  for  this  simple  example.  The  adaptive 
filter  set  is  initialized  with  beamsteer  weights: 


a^O)  = 1/3 


A = 


a2(0)  = 1/3 


a3(0)  = 1/3 


An  optimum  filter  set  is 


A = 


a^O)  = 1/2 
a2(0)  = 0 
a3(0)  = 1/2 


Whenever  a noise  spike  first  appears  on  channel  2,  the  filter  vector  A is  up 
dated  according  to  the  equation 


A(t  + At)  = A(t)  + 2cy{t)  (N-N) 


Channel  1 


Channel  3 


Channel  1 


Channel  3 — — — — 

(b)  Noise  Model 


FIGURE  1-4 

A SIMULATED  DISPLAY  OF  THE  TIME  SERIES  INPUT 
CORRESPONDING  TO  THE  ILLUSTRATIVE  EXAMPLE 


The  adaptive  update  equation  shifts  its  weighting  from  channel  2 to  channels  1 
and  3.  The  second  time  the  noise  spike  appears  on  channel  2,  the  filter  output 


is 


[0  1 o] 


1 2 ~ I 

T + 1TC 


The  amplitude  of  the  noise  output  is  diminished  whenever  0 < c < 3/2.  After 
k updates,  the  adaptive  filter  vector  is 


1 1 , . 4 k 

T ~ T (i  • Tc) 


4 .k 

T c 


T (i- 


4 k 

Tc) 


When  the  (k+l)-st  noise  spike  appears  on  channel  2,  the  output  of  the  adaptive 


beamformer  is 


y(t)  = 


T<> 


4 k 

Tc> 


Each  time  a new  noise  spike  appears  on  channel  2,  the  adaptive  filter  output  is 
( 1 - 4c/ 3)  times  its  previous  value.  After  each  new  noise  spike,  the  difference 
between  the  optimum  filter  vector  and  the  old  filter  vector  is  reduced  by  the 
same  factor  (1  - 4c/ 3).  In  this  simple  example  where  the  noise  crosscorrela- 
tion matrix  NN^  has  a single  non-zero  eigenvalue,  the  average  squared  filter- 
ed noise  output  y^(t)  is  minimized  when  c = 3/4.  In  this  case,  noise  is  com- 
pletely eliminated  starting  with  the  second  noise  spike  on  channel  2.  When 
0 < c < 3/4,  the  noise  output  always  has  the  same  sign  as  the  spike  on  channel  2. 
When  c > 3/4  , however,  the  noise  output  alternately  differs  and  agrees  in  sign 
with  the  input  spike  on  channel  2.  The  magnitude  of  the  output  spike  always  di- 
minishes as  long  as  0 < c < 3/2  . When  c > 3/2,  on  the  other  hand,  the  adap- 
tive filter  vector  diverges  and  the  ABF  output  increases  in  magnitude  with  the 
arrival  of  each  new  noise  spike. 

In  thj.s  idealized  example,  the  signal  does  not  affect  the  filter  update 

since 
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In  the  event  that  a signal  spike  appears  simultaneously  on  all  channels,  the  out- 
put from  the  beamformer  is  s since  the  filter  weights  rum  to  one. 


G.  ORGANIZATION  OF  REPORT 

Section  II  studies  the  effect  of  some  changes  made  since  the  previous 
adaptive  processing  study  in  the  procedures  used  to  implement  time-domain 
maximum  likelihood  adaptive  multichannel  filtering.  The  criterion  for  evaluat- 
ing these  changes  is  the  adaptive-beamforming  signal-to-noise  gain  relative  to 
beamsteering  for  unmixed  long-period  seismic  events  in  the  presence  of  back- 
ground noise.  The  first  topic  examined  is  the  effect  of  replacing  the  previously- 
employed  adaptive  update  algorithm,  which  speeds  up  in  the  presence  of  a sig- 
nal, with  a new  adaptive  update  algorithm  which  slows  down  in  the  same  situa- 
tion. Next,  using  the  new  adaptive  algorithm,  adaptive  beamforming  improve- 
ment over  beamsteering  for  the  full  ALPA  array  is  compared  with  that  obtained 
with  closely-spaced  six-channel  partial  ALPA  arrays.  Finally,  the  effectiveness 
of  a new  prefilter  for  the  individual  vertical-component  channels  input  to  the 
beamformers  is  assessed  in  terms  of  its  detection  capability.  This  assessment 
of  the  new  prefilter  utilizes  adaptive  filter  lengths  of  both  31  and  15  points  per 
channel  in  order  to  determine  whether  any  meaningful  processing  loss  occurs 
with  the  shorter  filter  length,  which  halves  the  computational  load.  This  sec- 
tion is  intended  to  aid  in  evaluating  the  potential  usefulness  of  adaptive  filtering 
systems  for  processing  long-period  arrays  similar  to  ALPA. 

Section  111  examines  the  possibility  that  adaptive  beamforming  can  de- 
tect signals  buried  in  oft- azimuth  interfering  events  when  the  conventional 


A 
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time-shift-and- sum  processor  is  unable  to  do  so.  To  simulate  mixed  events, 
two  separately-recorded  data  samples  are  scaled  and  summed  to  create  a com- 
posite sample  used  in  the  adaptive-filter  update  procedure.  Beam  outputs  are 
computed  for  both  data  samples  individually  to  ascertain  the  effect  of  the  adap- 
tive beamforming  process  on  each  sample  in  the  composite  data  sample.  Then 
the  beam  outputs  for  both  samples  in  the  composite  sample  are  summed  in  or- 
der to  observe  the  result  of  combining  the  two  samples.  The  detection  capa- 
bility of  adaptive  processing  is  determined  from  a visual  examination  of  the 
composite  beam-output  trace.  This  investigation  of  the  ability  of  adaptive 
multichannel  filtering  to  detect  signals  masked  by  off-azimuth  interfering  events 
employs  both  long-period  ALP  A data  and  Korean  short-period  data.  The  pre- 
liminary results  presented  in  this  section  should  assist  those  people  interested 
in  improved  separation  of  mixed  events. 

Section  IV  presents  the  conclusions  of  this  study.  For  those  with  in- 
sufficient time  to  analyze  in  detail  the  results  in  Sections  II  and  III  of  this  re- 
port, this  section  gives  the  highlights  of  this  investigation. 
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SECTION  II 

COMPARISON  OF  PROCESSING  PROCEDURES 


A.  INTRODUCTION 


This  section  presents  a comparison  of  various  techniques  for  imple- 
menting adaptive  beamforming.  The  previous  adaptive  processing  study 
(Barnard,  1973)  employed  the  adaptive  filter  update  algorithm 


A(t  + At) 


2*iXTA  (X-X) 

— T — ' 

(X-X)  (X-X) 


which  speeds  up  in  the  presence  of  an  on-azimuth  signal.  The  first  topic  of 
investigation  is  the  effect  on  signal-to-noise  ratio  improvement  relative  to 
beamsteering  when  the  old  adaptive  algorithm  is  replaced  by  a new  adaptive 
filter  update  algorithm 


A(t  + At) 


A + 


2jiXTA  (X-X) 
T 

XX 


» 


which  slows  down  in  the  presence  of  an  on-azimuth  signal.  The  next  subject 
of  examination  is  the  performance  of  the  new  adaptive  algorithm  compared 
with  beamsteering  when  all  available  sensors  at  ALPA  are  input  to  the  time- 
shift-and-sum  and  adaptive  beamformers,  Of  particular  interest  are  the  dif- 
ferences between  the  full-array  results  and  the  results  obtained  with  closely- 
spaced  six-channel  arrays  such  as  the  ones  used  to  compare  the  new  adaptive 
algorithm  with  the  old  adaptive  algorithm.  Finally,  the  effect  of  a new  spe- 
cially designed  signal-enhancement  prefilter  on  detection  performance  and 
adaptive  processing  signal-to-noise  gain  over  beamsteering  is  compared  with 
that  of  the  old  prefilter  applied  during  the  previous  contract  to  the  site 


vertical-component  data  entering  the  beamformers.  The  new  prefilter  is 
evaluated  using  adaptive  filter  lengths  of  31  and  15  points  per  channel  in  or- 
der to  determine  whether  the  shorter  filter  length  produces  any  meaningful 
processing  loss. 

All  results  in  this  section  apply  to  individual  unmixed  seismic  events 
in  the  presence  of  background  noise.  To  maintain  compatibility  with  the  pro- 
cessing results  obtained  previously  for  the  old  adaptive  algorithm,  adaptive 
filter  sets  are  designed  individually  for  the  data  sample  containing  back- 
ground noise  and  for  the  data  sample  used  in  signal-degradation  measure- 
ments. The  criterion  of  performance  in  evaluating  the  various  techniques 
for  implementing  adaptive  beamforming  is  the  estimated  signal-to-noise  ratio 
improvement  relative  to  beamsteering.  The  method  of  approach  used  to  esti- 
mate signal-to-noise  ratio  improvement  is  to  measure  the  difference  between 
noise  reduction  and  signal  degradation.  Both  noise  reduction  and  signal  de- 
gradation are  computed  in  terms  of  the  drop  in  power  for  the  adaptive  filter 
output  relative  to  the  beamsteer  output.  Since  it  is  impossible  in  the  case  of 
real  data  to  distinguish  between  signal  power  and  noise  power  when  a signal 
is  present,  the  data  samples  for  signal  measurements  are  different  from 
those  used  for  noise  measurements.  Noise  reduction  is  computed  over  data 
samples  approximately  four  hours  long.  Signal  degradation,  on  the  other 
hand,  is  computed  over  gates  approximately  five  minutes  long.  In  the  case 
of  signal  degradation  measurements,  the  effect  of  adaptive  filtering  on  con- 
taminating noise  is  lumped  together  with  the  effect  on  the  signal  proper.  As 
long  as  the  adaptive  filter  set  reduces  the  contaminating  noise  power  at  least 
ar  much  as  the  signal  power,  the  signal  degradation  is  not  underestimated 
aud  the  signal-to-noise  ratio  improvement  is  not  overestimated.  Noise  re- 
duction and  signal  degradation  are  calculated  for  the  same  pair  of  data  sam- 
ples over  a range  of  convergence  rates.  From  these  values,  the  signal-to- 
noise  gain  is  estimated  as  a function  of  convergence  rate. 
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Subsection  B is  a mathematical  discussion  of  the  reasons  why  the  new 
adaptive  algorithm  should  outperform  the  old  adaptive  algorithm.  Subsection 
C is  an  empirical  comparison  of  the  two  adaptive  algorithms.  In  Subsection 
D,  full-array  processing  results  and  six-channel  processing  results  are  com- 
pared for  the  new  adaptive  algorithm  when  the  site  vertical  components,  be- 
fore input  to  the  beamformers,  are  passed  through  the  old  prefilter  used 
during  the  previous  study,  A heuristic  development  of  the  techniques  used 
in  designing  a new  signal-enhancement  prefilter  for  detection  purposes  ap- 
pears in  Subsection  E.  The  effect  of  applying  the  new  prefilter  to  the  vertical- 
component  data  is  examined  in  Subsection  F for  adaptive  filter  lengths  of  31 
and  15  points  per  channel.  Finally,  Subsection  G summarizes  the  results 
obtained  with  the  adaptive  beamforming  techniques  studied  in  this  section. 

B.  CHARACTERISTICS  OF  SELECTED  MAXIMUM  LjKELIHCOD 

ADAPTIVE  ALGORITHMS 

The  basic  maximum  likelihood  adaptive  algorithm 
A(t  + At)  = A+  2^XTA(X-X) 

T f T1 

attempts  to  minimize  the  mean  square  adaptive  filter  output  A ELXX  JA 

. . ,,  . M * 

subject  to  the  signal-preservation  constraint  conditions  ^ a-(j)  = ° • 

i=  1 1 °J 

(j  = -N,  . . . , -1,  0,  1,...,  N).  When  a constant,  data-independent  conver- 
gence parameter  jU  is  used  to  process  seismic  data,  this  basic  algorithm 
has  two  noteworthy  disadvantages. 

First,  because  of  seismic -event  amplitudes  often  spanning  several 
orders  of  magnitude,  the  convergence  parameter  must  be  selected  on  the 
basis  of  the  specific  situations  encountered  in  a particular  data  sample.  The 
convergence-parameter  range  within  which  no  divergence  occurs  is  a func- 
tion of  those  sections  of  data  where  the  power  is  greatest.  Likewise,  the 
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convergence-parameter  range  where  signal-to-noise  gain  is  maximized  or- 
dinarily differs  from  data  sample  to  data  sample.  Even  after  the  develop- 
ment of  suitable  convergence-rate  normalization  and  power-measurement 
procedures,  the  adaptive  filter  vector  A in  many  cases  moves  exceedingly 
slowly  in  the  presence  of  the  only  stable  feature  of  the  noise  field-background 
noise. 


The  second  important  drawback  of  the  basic  algorithm  is  the  undue 
emphasis  it  places  on  transient  bursts  of  energy  from  seismic  events.  Fre- 
quently, in  the  case  of  long-period  data,  events  are  the  dominant  power 
source.  After  an  event  dies  out,  a new  event  from  the  same  general  direc- 
tion cannot,  in  general,  be  expected  for  some  time.  Immediately  after  an 
event,  therefore,  the  data  is  no  longer  quasi-stationary:  the  exponentially- 
smoothed  time  average  for  the  matrix  XX  does  not  approximate  the  prob- 
abilistic matrix  expectation  eCxX  3 at  the  instant  of  filter  application.  As 
a result,  adaptive  beamforming  cannot  achieve  near-optimum  noise  suppres- 
sion until  exponential  damping  reduces  the  discrepancy  between  the  time  aver- 
age  and  the  probabilistic  expectation  to  an  acceptable  level.  After  large  events, 
this  exponential-decay  process  may  be  lengthy.  In  seismic  processing,  the 
basic  algorithm  can  be  improved  by  emphasizing  the  stable  background  noise 
field  at  the  expense  of  seismic  events. 

In  the  previous  adaptive  processing  study  (Barnard,  1973),  the  modi- 
fied adaptive  algorithm 


A(t  + At)  = A+  2^X  A (X-X) 

(X-X)T(X-X) 

was  employed.  This  algorithm  attempts  to  minimize  the  quantity 
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subject  to  the  maximum  likelihood  constraints  2 a.(j)  = 8 ..  A disadvan- 

J 1=1  i J oj 

tage  of  this  algorithm  is  that  each  update  reduces  the  power  output  for  the  old 

2 

data  vector  by  a factor  (1  - 2 y)  regardless  of  the  fact  that  the  data  vector 

X may  be  consistent  with  a signal: 


XTA(t  + At) 


2/iXT(X-X)  1 
(X-X)T(X-X)  J 

2M(X-X)T(X-X)  1 
(X-X)T(X-X)  J 


= XTA  (1  - 2n), 


so  that 

AT(t +At)  XXTA(t  +At)  = ATXXTA(1  - 2/i)2. 


Furthermore,  to  prevent  the  squared  magnitude 

i A(t  + At)  - ai2  = J?-  y-ffl- 

(X-X)  (X-X) 

of  the  filter  update  vector  for  this  algorithm  from  becoming  enormous  when 
a strong  signal  traverses  the  array,  the  adaptive  filter  set  must  be  frozen 
whenever  a signal  is  detected.  A filter-freeze  procedure  using  the  Fisher 
detector  algorithm 

- _ (M-1)XTX 

F - T — 

(X-X)  (X-X) 

was  used  in  the  previous  study  to  curb  excessive  movement  of  the  filter  vec- 
tor A. 


It  0 
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In  the  current  9tudy,  the  adaptive  algorithm 

A(t  + At)  = A + 2MXTA(X-X) 

T 

X X 

is  utilized  in  an  attempt  to  minimize  the  quantity 


subject  to  the  same  maximum  likelihood  constraints.  With  this  algorithm, 
the  squared  magnitude 


|A(t+At)-A|2  = 44Zy2(t)  (X-X)T(X-X) 

T 2 

(xxxr 

of  the  filter  update  vector  decreases  when  a strong  signal  propagates  across 
the  array.  If  the  updated  filter  vector  were  applied  to  the  old  data  vector, 
moreover,  the  adaptive  filter  output  would  be 


XTA(t+At)  = XTA  |l+  2M  X (X-X) 

' xTx 

2m(X-X)T(X-X)1 

XTX  J , 

so  that  the  squared  filter  output  would  be  reduced  by  a factor  of 
[l  - 2M{X-X)T(X-X)/XTX]2  instead  of  (1  - . In  conformity  with  the 

design  goals  of  the  maximum  likelihood  multichannel  filter,  there  would  be 
no  attenuation  at  all  in  the  case  of  an  ideal  signal.  With  weak  signals  having 
signal-to-noise  ratios  too  low  to  trigger  a filter  freeze,  this  algorithm  ap- 
pears to  have  a definite  advantage  over  the  algorithm  previously  employed: 
the  convergence  parameter  fi  for  the  new  algorithm  would  probably  be  sig- 
nificantly greater  than  for  the  old  algorithm  if  signal  degradation  were  the 
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same  for  both  algorithms,  so  that  greater  noise  reduction  and  hence  greater 
signal-to-noise  gain  might  be  achieved  by  the  new  algorithm.  The  new  algo- 
rithm, like  the  old  one,  would  also  be  less  sensitive  to  power  bursts  from 
seismic  events,  so  that  non-repeating  seismic  events  would  be  forgotten 
more  rapidly  than  with  the  basic  maximum  likelihood  algorithm 


A(t  + At)  = A+  2/iXTA(X-X). 


Normalization  of  the  data  vector  X by  its  absolute  value  | X | in  the  new 
algorithm  instead  of  the  absolute  value  | X-X  I , furthermore,  is  sensibly 
motivated:  the  adaptive  filter  responds  to  the  phase  characteristics  of  the 
data  vector  X and  not  to  its  magnitude.  It  is  difficult  to  describe  precisely 
the  effect  of  the  normalizing  factor  | X-X  I for  the  old  algorithm. 

C.  COMPARISON  OF  THE  NEW  ADAPTIVE  ALGORITHM  WITH  THE 
OLD 

This  subsection  contains  an  empirical  comparison  of  the  adaptive 


algoJt  ithm 


A (t  + At)  = A + 


2jiXTA(X-X) 

T 

X X 


used  in  this  study  with  the  algorithm 

2jiXTA(X-X) 

A (t  + At)  = A+  — YZ 

(X-X)  (X-X) 

employed  during  the  previous  contract.  In  order  to  make  a meaningful  com- 
parison of  the  two  adaptive  algorithms,  the  same  two  data  samples  which  pro- 
vided the  critical  results  of  the  previous  investigation  serve  as  input  to  the 
new  adaptive  algorithm.  The  principal  criterion  for  evaluating  algorithm 
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performance  is  the  signal-to -noise  ratio  improvement  of  the  adaptive  beam- 
former  relative  to  beamsteering.  The  approach  used  to  estimate  signal-to- 
noise  ratio  improvement  is  to  measure  the  difference  between  noise  reduction 
and  signal  degradation.  Both  noise  reduction  and  signal  degradation  are  com- 
puted in  terms  of  the  drop  in  power  for  the  adaptive  filter  output  relative  to 
the  beamsteer  output.  For  consistency,  adaptive  filter  sets  are  designed  in- 
dividually (as  before)  for  the  data  sample  containing  background  noise  and  for 
the  data  sample  used  in  signal-degradation  measurements.  The  results  pre- 
sented give  signal-to- noise  ratio  improvement  as  a function  of  convergence 
rate.  Results  for  the  old  algorithm  come  from  the  previous  study,  while  re- 
sults for  the  new  algorithm  originate  from  work  performed  during  the  current 
contract. 

To  obtain  noise  reduction  measurements  for  the  old  algorithm,  a 
noise  sample  spanning  the  period  0800  to  1145  on  day  238  of  1970  is  utilized. 
The  steer  direction  corresponds  to  an  azimuth  of  270°  and  an  apparent  velo- 
city of  3.5  km/sec.  Noise  reduction  is  computed  as  a function  of  the  conver- 
gence parameter  /!  over  the  entire  frequency  band  0 to  0.  5 Hz.  The  vertical 
components  of  sites  1,  2,  3,  6,  8,  and  9 from  AL.PA  are  passed  through  a 
prefilter  with  the  frequency  response  shown  in  Figure  II- 1 before  input  to  the 
adaptive  beamformer.  (All  beamforming  results  in  this  'iubsection  and  the 
next  employ  data  passing  through  the  prefilter.  ) The  PDE  bulletin  lists  no 
events  between  0639  and  1502  on  August  26  of  1970.  A small  signal  reached 
ALPA  at  1021,  however,  from  an  azimuth  just  under  60°.  The  energy  from 
this  signal  is  only  a small  fraction  of  the  total  energy  within  the  entire  sam- 
ple and  should  have  only  a minor  effect  on  the  reported  noise  reduction  values. 
Table  II- 1 and  Figure  II-2  present  the  noise  reduction  in  dB  for  the  old  algo- 
rithm as  a function  of  the  convergence  parameter  fj,  . 

To  obtain  signal  degradation  measurements  for  the  old  algorithm, 
three  signals  were  selected  out  of  a data  sample  from  2000  to  2345  on  day 
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TABLE  II- 1 

ADAPTIVE  FILTERING  NOISE  REDUCTION  VERSUS 
CONVERGENCE  RATE  FOR  OLD  ALGORITHM 
(ALPA,  6 SITES,  1970  DAY  238  0800-1  145) 


Convergence  Rate 

Noise  Reduction 

(M) 

(dB) 

0.  002 

1.099 

0.  003 

1.  359 

0.  004 

1.521 

0.  005 

1.639 

0.  007 

1. 813 

0.  010 

1.  983 

0.  015 

2.  186 

0.  020 

2.  346 

0.  030 

2.635 

0.  050 

3.  135 

0.  070 

3.  562 

0,  100 

4.  056 

0.  150 

4.547 

0.  200 

4.786 

0.  250 

4.  892 

0.  300 

4.  940 

0.  350 

4.  945 

0.  400 

4.  939 

0.  500 

4.  896 

MMMNMIIWWVB 
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276  of  1971  (October  3).  Several  events  from  the  PDE  bulletin  and  LASA 
bulletin  had  arrival  times  within  this  period.  Table  II-2  gives  the  PDE  events 
which  arrived  during  the  sample  interval.  In  addition  to  these  events,  two 
events  located  within  three  degrees  of  the  Panama  event  were  included  in  the 
LASA  bulletin  for  this  time  interval.  The  first  of  these  had  a measured 
LASA  bodywave  magnitude  of  3.6  and  arrived  at  LASA  47  minutes  and  34  sec- 
onds before  the  PDE  event  from  Panama.  The  second  had  a measured  LASA 
bodywave  magnitude  of  3.9  and  arrived  at  LASA  6 minutes  and  38  seconds  be- 
fore the  PDE  event  from  Panama.  In  processing  the  three  signals  selected, 
sites  8,  9,  12,  15,  16,  and  17  are  input  to  the  time-shift-and-sum  and  adap- 
tive filter  beams. 

A weak  signal  arrived  at  ALPA  at  approximately  2138.  The  signal  - 
to-noise  ratio  on  the  beamsteer  output  is  about  6 dB.  No  event  corresponding 
to  this  signal  can  be  found  on  either  the  PDE  bulletin  or  LASA  bulletin.  Fre- 
quency-wavenumber spectra.  Fisher-detector  scans,  and  beamsteer-output 
power  comparisons  indicate  an  arrival  azimuth  between  300°  and  305°.  The 
beamsteer  and  adaptive-filter  beams  are  aimed  toward  302.  5°.  Table  II-3 
and  Figure  II-3  depict  signal  degradation  for  the  old  algorithm  as  a function 
of  convergence  rate  over  the  frequency  band  0 to  0.  5 Hz.  The  negative  de- 
gradation at  (J=  0.  002  means  that  the  adaptive  filter  trace  contains  more 
energy  than  the  beamsteer  output  trace  over  the  512-point  computation  gate. 
Combining  these  signal  degradation  measurements  with  the  noise  reduction 
measurements  from  day  2 38  of  1970  yields  an  estimate  of  the  signal-to-noise 
gain  which  would  have  occurred  on  day  238  of  1970  if  the  weak  signal  had  ar- 
rived at  ALPA  during  the  period  0800-1145.  Table  II-4  and  Figure  11-4  con- 
tain the  resulting  signal-to-noise  gain  estimate  as  a function  of  convergence 
rate.  Best  signal-to-noise  ratio  improvement  is  obtained  at  H = 0.004. 

The  remaining  two  signals  used  in  signal-degradation  measurements 
are  large  signals  from  near  the  east  coast  of  Kamchatka.  The  first  of  these 
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TABLE  II- 2 


I 

I 

I 


:i 


fi 


i. 


PDE  EVENTS  ARRIVING  AT  ALPA 
BETWEEN  2000  AND  2345  ON  OCTOBER  3,  1971 


Origin  Time 
(Hr  -Min-Sec) 

Lat 

(Deg) 

I'Ong 

(Deg) 

Region 

mb 

Azimuth 

(Deg) 

20:42:46.  7 

3.  5N 

82.  9 W 

South  of  Panama 

4.  7 

113 

20:  54:48.  5 

55.  7N 

162.  IE 

Near  East  Coast 
of  Kamchatka 

4.  5 

27  3 

21:  54: 12.  9 

55.  8N 

162.  2E 

Near  East  Coast 
of  Kamchatka 

4.  9 

27  3 

22:34:54.  8 

4.  2S 

152.  7E 

New  Britain 
Region 

4.  8 

240 

23:17:26.4 

51.  8N 

173.  3 W 

Andreanof  Islands, 
Aleutian  Is 

4.  3 

237 
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TABLE  II-3 


ADAPTIVE  FILTERING  SIGNAL  DEGRADATION  VERSUS 
CONVERGENCE  RATE  WITH  OLD  ALGORITHM  FOR 
A WEAK  SIGNAL  FROM  300°- 305° 

(ALPA,  6 SITES,  1971  DAY  276  21.  37.01-21.45.  32) 


C 


onvergence  Rate 


Signal  Degradation 


(M) 


(dB) 


0.  002 
0.  003 
0.  004 
0.  005 
0.  007 
0.  010 
0.  015 
0.  020 
0„  030 
0.  050 
0.  070 
0.  100 
0.  150 
0.  200 
0,  250 


-0. 024 
0.  104 
0.  255 
0o  408 
0.650 

0.  945 

1.  301 
1.  569 

1.  970 

2.  562 
30  049 

3.  713 

4.  505 
4.973 

5.  176 


<> 
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FIGURE  II- 3 


SIGNAL  DEGRADATION  AS  A FUNCTION  OF  CONVERGENCE 


JiVJi’il-*-'  ^ * — 

RATE  WITH  OLD  ALGORITHM  FOR  A WEAK  SIGNAL  FROM 
300°  - 305°  {ALPA,  6 SITES,  1971  DAY  276  21 . 37.  01 -21 . 45.  32) 
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TABLE  H-4 

ADAPTIVE  FILTERING  SIGNAL-TO-NOISE  GAIN  VERSUS 
CONVERGENCE  RATE  WITH  OLD  ALGORITHM  FOR 
A WEAK  SIGNAL  FROM  300°- 305° 

(ALPA,  6 SITES,  1971  DAY  276  21.  37.  01  - 21.45.  32, 
USING  NOISE  FROM  1970  DAY  238  0800-1145) 


Convergence  Rate 


Signal-to-Noise  Gain 


0.  002 
0.003 
0.  004 
0.005 
0.007 
0.010 
0.015 
0.020 
0.0  30 
0.  050 
0.  070 
0.  100 
0.  150 
0.200 
0.250 


1.  123 
1. 255 
1.  266 
1. 231 
1.  163 
1. 038 
0.  885 
0.  777 
0.  665 
0.  573 
0.  513 
0.  343 
0.  042 
-0. 187 
-0.  284 
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FIGURE  II -4 


SIGNAL -TO-NOISE  GAIN  AS  A FUNCTION  OF  CONVERGENCE 
RATE  WITH  OLD  ALGORITHM  FOR  A WEAK  SIGNAL  FROM 
300°-305°  (ALPA,  6 SITES,  1971  DAY  276  21.37.01-21.45.32, 
USING  NOISE  FROM  1970  DAY  238  0800-1145) 
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has  a bodywave  magnitude  of  4.  5,  an  origin  time  of  20:54:48.  5,  a latitude  of 
55.7  N,  and  a longitude  of  162.  1 E.  The  azimuth  with  respect  to  ALP  A is 
273°.  Accordingly,  the  beams  are  trained  in  this  direction.  Eighteen  dB 
is  the  estimated  signal-to-noise  ratio  on  the  beamsteer  output.  Table  11-5 


and  Figure  11-5  show  the  signal  degradation  for  this  signal  using  the  old  algo- 
rithm. The  computation  gate  for  this  signal  is  256  seconds  instead  of  512 


seconds  (as  in  the  case  of  the  weak  signal  from  300°-305°).  The  resulting 


estimated  signal-to-noise  ratio  improvement  using  noise  reduction  measure- 
ments from  day  238  of  1970  appears  in  Table  11-6  and  Figure  11-6.  Best 
ABF  performance  again  occurs  at  H = 0.  004. 


The  last  signal  used  to  measure  signal  degradation  is  a magnitude  4.  9 
event  from  almost  exactly  the  same  location.  The  PDE  bulletin  gives  the 
origin  time  as  21:54:12.9,  the  latitude  as  55.8  N,  and  the  longitude  as  162.2  E. 
Again  the  beams  are  directed  toward  a 273°  azimuth.  For  this  event,  the 
signal-to-noise  ratio  on  the  beamsteer  output  is  close  to  24  dB.  Signal  de- 
gradation over  a 256-point  interval  for  the  old  algorithm  is  portrayed  in  Table 
11-7  and  Figure  11-7.  It  is  considerably  greater  than  for  the  event  with  an  18- 
dB  signal-to-noise  ratio.  Using  the  day  238  noise  sample  for  noise  reduction 
measurements,  the  resulting  signal-to-noise  gain  for  the  old  algorithm  is  re- 
presented in  Table  11-8  and  Figure  11-8.  Max!mum  improvement  is  obtained 
at  the  slowest  convergence  rate. 


The  principal  result  in  processing  the  three  signals  with  the  old  algo- 
rithm is  that  the  signal  degradation  increases  as  the  signal-to-noise  ratio  in- 
creases. With  the  old  algorithm,  at  least,  the  maximum  likelihood  constraints 
do  not  guarantee  preservation  of  the  signal.  When  no  noise  is  present,  the 


constraint  conditions  and  differences  in  successive  realizations  of  the  data 


vector  X ultimately  limit  the  signal-degradation  levels  which  are  possible. 
As  the  influence  of  noise  becomes  significant,  the  resultant  retardation  of  the 
effective  adaptation  rate  reduces  the  measured  signal-degradation  values. 
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TABLE  U-5 

ADAPTIVE  FILTERING  SIGNAL  DEGRADATION  VERSUS 
CONVERGENCE  RATE  WITH  OLD  ALGORITHM  FOR 
A STRONG  SIGNAL  FROM  KAMCHATKA 
(ALPA,  6 SITES,  1971  DAY  276  21.07.  18-2i.l  1.33) 


Convergence  Rate 

m 

Signal  Degradation 
(dB) 

0.  002 

0.  269 

0C  003 

0.416 

0.  004 

0.556 

0.  005 

0.693 

0.  007 

0.  944 

0.  010 

1.279 

0.  015 

1.743 

0.  020 

2.  137 

0.  030 

2.  718 

0.  050 

3.  423 

0.  070 

3.775 

0.  100 

3.  989 

0.  150 

4.  119 

0.  200 

4.  203 

0.  250 

4.  243 

Convergence  Rate  /i 


FIGURE  II- 5 

SIGNAL  DEGRADATION  AS  A FUNCTION  OF  CONVERGENCE 
RATE  WITH  OLD  ALGORITHM  FOR  A STRONG  SIGNAL  FROM 
KAMCHATKA  (ALPA,  6 SITES,  1971  DAY  276  21. 07.  18-21.  11.  33) 


TABLE  11-6 


I 

I 
I 

j ADAPTIVE  FILTERING  SIGNAL-TO- NOISE  GAIN  VERSUS  CONVERGENCE 

| RATE  WITH  OLD  ALGORITHM  FOR  A STRONG  SIGNAL  FROM 

KAMCHATKA  (ALP A,  6 SITES,  197  1 DAY  276  2 1.  07.  18- 
- 21.11.33,  USING  NOISE  FROM  1970  DAY  238  0800-1  145) 


«* 

Convergence  Rate 
(M) 

Signal-Lo-Noise  Gain 
(dB) 

— 

0.  002 

0.  830 

» i 

0.  003 

0.943 

-- 

0.  004 

0.  965 

• • 

0.  005 

0.  946 

t * 

i i 

0.  007 

0.  869 

a 

0.  010 

0.  704 

0.  015 

0.443 

i. 

0.  020 

0.  209 

0.  030 

-0.  083 

m ■ 

0.  050 

-0.  288 

T" 

0.  070 

-0.  213 

• m 

0.  100 

0.  067 

* • 

0.  150 

0.428 

• « 

0.  200 

0.  589 

* *i 

0.  250 

0.  649 
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Gain  (dB) 


TABLE  II- 7 

ADAPTIVE  FILTERING  SIGNAL  DEGRADATION  VERSUS 
CONVERGENCE  RATE  WITH  OLD  ALGORITHM  ^OR 
A VERY  STRONG  SIGNAL  FROM  KAMCHAIKA 
(ALPA,  6 SITES,  1971  DAY  276  22.  07.02-22,  11.  17) 


Convergence  Rate 

Signal  Degradation 

(M) 

(dB) 

0.002 

0.  963 

0.  003 

1.  324 

0.  004 

1.655 

0.  005 

1.  977 

0.  007 

2.571 

0.  010 

3.  360 

0.  015 

4.450 

0.  020 

5.294 

0.  030 

6.  367 

0.050 

7.  179 

0.  070 

7.418 

0.  100 

7.461 

0.  150 

7.465 

0.  200 

7.495 

0.  250 

7.497 

De 


a 


o] 


Convergence  Rate  /< 


li  BilP 


FIGURE  II -7 

SIGNAL  DEGRADATION  AS  A FUNCTION  OF  CONVERGENCE 
RATE  WITH  OLD  ALGORITHM  FOR  A VERY  STRONG  SIGNAL 
FROM  KAMCHATKA  (ALPA,  6 SITES,  1971  DAY  276  22.  07.  02-22.  1 1.17) 


TABLE  H-8 


I 


ADAPTIVE  FILTERING  SIGNAL-TO-NOISE  GAIN  VERSUS 
CONVERGENCE  RATE  WITH  OLD  ALGORITHM  FOR 
A VERY  STRONG  SIGNAL  FROM  KAMCHATKA 
(ALPA,  6 SITES,  1971  DAY  276  22.07.  02-22,11.17, 
USING  NOISE  FROM  1970  DAY  238  0800-1  145) 


Convergence  Rate 

Signal-to-Noise  Gain 

00 

(dB) 

0.  002 

0.  136 

0.  003 

0.  035 

0.  004 

-0.134 

0.  005 

-0.  338 

0.007 

-0.758 

0.  010 

-1. 377 

0.  015 

-2.  264 

0.  020 

-2.  948 

0.  030 

-3.  732 

0.  050 

-4. 044 

0.  070 

-3.  856 

0.  100 

-3.405 

0.  150 

-2.  918 

0.  200 

-2.  709 

0.  250 

-2. 605 
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FIGURE  II -8 

SIGNAL -TO-NOISE  GAIN  AS  A FUNCTION  OF  CONVERGENCE 
RATE  WITH  OLD  ALGORITHM  FOR  A VERY  STRONG  SIGNAL  FROM 
KAMCHATKA  (ALPA,  6 SITES,  1971  DAY  276  22.  07.  02-22.  11.17, 
USING  NOISE  FROM  1970  DAY  238  0800-1145) 
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Sections  IV  and  V of  the  previous  report  (Barnard,  1973)  contain  a 
mors  comprehensive  discussion  of  the  results  for  the  old  algorithm  and  pre- 
sent several  Calcomp  plots  of  the  beamsteer  and  adaptive  filter  output  traces 
from  the  data  samples  processed  so  far  in  this  subsection.  To  conserve 
space,  they  do  not  appear  in  this  report. 

In  the  computer  program  implementing  the  old  algorithm,  the  com- 
ponents of  the  adaptive  filter  vector  A are  expressed  in  16-bit  integer  form. 

- T - 

In  addition,  the  denominator  term  (X-X)  (X-X)  is  approximated  by  the  quan- 
M 

tity  (2N+1)  2 P.(t)  , where  P.(t)  in>  an  exponentially-smoothed  power 

i=l  1 1 

average  for  the  i-th  channel.  P.(t)  is  computed  by  the  formula 


P.(t)  = (1  - € ) P.(  t - At)  + c[x(t)  - x.(t)]  , 

where  e is  an  input  parameter.  Initially,  P^(t)  is  zero,  and  several  values 

of  P.(t)  are  computed  before  the  filter  vector  A is  allowed  to  update.  By 

way  of  contrast,  in  the  computer  program  for  the  current  study,  computations 

utilize  standard  IBM  360  32-bit  floating-point  arithmetic,  and  the  denominator 
T 

term  X X is  the  true  dot  product  of  the  instantaneous  data  vector  X with 
itself.  The  inability  of  the  old  algorithm  at  high  convergence  rates  to  achieve 
noise  reduction  values  as  great  as  those  for  the  new  algorithm  may  be  due  to 
the  computational  differences  described  here. 

Noise  reduction  measurements  for  the  new  algorithm  come  from  the 
same  data  sample  used  in  noise  reduction  measurements  for  the  old  algorithm. 
The  sample  extends  from  0800  to  1145  on  day  238  of  1970.  The  only  difference 
in  processing  is  that  an  apparent  velocity  of  3.  6 km/ sec  is  specified  instead 
of  the  3.5  km/ sec  velocity  employed  with  the  old  algorithm.  The  3,6  km/ sec 
velocity  is  probably  a slightly  more  realistic  value  and  is  the  default  value  in 
the  current  computer  program.  The  3.5  km/ sec  velocity  was  used  in  the  old 
computer  program  to  obtain  the  results  for  the  old  adaptive  algorithm.  The 
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slight  difference  in  velocity  should  have  only  a minor  effect  on  the  processing 
results.  Table  II- 9 and  Figure  II- 9 give  the  noise  reduction  measurements 
for  the  new  algorithm  as  a function  of  the  convergence  parameter  /i  . The 
adaptation  rate  is  always  slower  with  the  new  algorithm  than  with  the  old 
algorithm  at  any  particular  value  of  the  convergence  parameter  H . The  rea- 
son is  that  the  denominator  term 

(X-X)T(X-X)  = xTx  - xTx 

f 

of  the  old  algorithm  is  always  smaller  than  the  corresponding  term  X X of 
the  new  algorithm.  For  spatially  uncorrelated  noise  of  equal  strength  on  all 
sensors,  the  beamsteer  output  power  is  1/M  times  the  average  power  across 
the  sensors  (where  M is  the  number  of  sensors).  Consequently,  the  effec- 
tive adaptation  rate  of  the  new  algorithm  is  (M-l)/M  times  the  rate  of  the 
old  algorithm  in  that  case.  For  propagating  plane  wave  energy  reasonably 
far  from  the  steer  direction,  the  same  is  true  in  an  average  sense.  In  any 
specific  case,  however,  the  result  depends  on  the  array  geometry,  the  fre- 
quency content  of  the  incoming  energy,  the  steer  direction,  and  the  direction 
of  the  incoming  energy.  For  an  ideal  signal  (where  X = X),  the  filter  vector 
A stops  completely  with  the  new  algorithm,  whereas  the  step  size 
lA(t  +At)  - A(t)  | becomes  infinite  with  the  old  algorithm.  To  compare  noise 
reduction  measurements  between  algorithms,  the  convergence  parameter 
values  [l  for  the  old  algorithm  need  to  be  multiplied  by  M/(M-1),  or  1.2,  to 
obtain  the  most  reasonable  equivalent  values  for  the  new  algorithm.  That  is 
to  say,  the  value  = 0.  05  for  the  old  algorithm  corresponds  to  the  value 
= 0.  06  for  the  new  algorithm.  At  M = 0.  5,  the  noise  reduction  for  the  new 
algorithm  is  6,  5 dB  and  is  still  increasing.  At  no  convergence  rate  did  the 

noise  reduction  for  the  old  algorithm  exceed  5 dB.  Roundoff  error  and  esti- 

— t _ 

mation  of  the  denominator  term  (X-X)  (X-X)  with  an  exponentially-smoothed 
power  average  may  account  for  the  lower  noise  reduction  values  with  the  old 
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TABLE  H-9 


ADAPTIVE  FILTERING  NOISE  REDUCTION  VERSUS 
CONVERGENCE  RATE  FOR  NEW  ALGORITHM 
(ALPA,  6 SITES,  1 970  DAY  238  0800- 1 145) 


Convergence  Rate 

Noise  Reduction 

(M)  ! 

(dB) 

0. 0010 

0„  641 

0. 0015 

0.  825 

0. 0020 

0.  969 

0.  0030 

1.  183 

0.  0040 

1.  334 

0. 0050 

1 o 447 

0. 0070 

1.60° 

0o  0100 

1.  772 

0. 0150 

1.960 

0. 0200 

2.  106 

0. 0300 

2„  352 

0. 0400 

2.  571 

0.  0500 

2.  775 

0. 0700 

3.  146 

0.  1000 

3.628 

0.  1500 

4.  292 

0. 2000 

4.  829 

0. 3000 

5.644 

0.  4000 

6.  193 

0.  5000 

6.  527 

-29 


algorithm.  At  high  convergence  rates,  the  major  roundoff-error  problem 

T 

seems  to  be  in  the  adaptive  filter  output  y(t)  = X A.  Error  in  the  data  vec- 
tor X may  change  the  sign  of  y(t).  In  addition,  rounding  the  adaptive  filter 
output  y(t)  to  the  nearest  integer  may  result  in  a zero  value  and  prevent  any 
update  of  the  filter  vector  A . A more  detailed  discussion  of  this  situation 
appears  in  Subsection  III-D  of  the  previous  report  (Barnard,  1973).  Other 
factors  contributing  to  the  difference  in  performance  for  the  two  algorithms 
may  be  the  clipping  of  the  filter  weights  at  an  absolute  value  of  one  half  (with 
the  old  algorithm)  and,  of  course,  the  algorithms  themselves. 

In  the  noise  sample  from  day  238  of  1970,  a weak  signal  reached 
ALPA  at  1021  from  an  azimuth  of  approximately  60°.  Figure  II- 10  shows 
the  beamsteer  output  and  adaptive -filter  output  with  the  new  algorithm  at 
fi=  0.  005  and  M=  0.  5 for  this  event.  Tick  marks  at  the  zero  levels  of  all 
three  beams  are  spaced  5 minutes  apart  and  correspond  to  the  time  immedi- 
ately above  the  middle  trace.  The  beamsteer  output  appears  in  the  top  trace 
the  adaptive -filter  output  for  M = 0.  005  in  the  middle  trace,  and  the  adaptive- 
filter  output  at  H = 0.  5 in  the  bottom  trace.  At  the  lower  convergence  rate, 
the  adaptive -filter  output  is  slightly  attenuated  relative  to  the  beamsteer  out- 
put. In  contrast,  the  adaptive -filter  output  at  H = 0.  5 in  the  bottom  trace  is 
negligible  in  comparison  with  the  beamsteer  output.  At  slow  convergence 
rates,  the  effect  of  the  beamsteer  and  adaptive  beamformers  is  reflected  in 
the  array  wavenumber  response 

i27r(fr.  - k • x.) 

A(f.i?)  = 2-j  Aj(f)  e J J 
j = l 

for  the  time-shift-and-sum  processor  and  the  adaptive  processor.  The  re- 
sponse A(f , lT)  is  the  complex  amplification  factor  for  a hypothetical  plane 
wave  of  frequency  f and  wavenumber 
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EVENT  FROM  60 

(DAY  238  1970,  STEER  DIRECTION  27 o' 
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= |- y sin  0,  - cos  0 j 


where  V is  the  apparent  velocity  across  the  array  and  0 is  the  azimuth 
from  which  the  plane  wave  arrives.  In  the  formula  for  A(f , k)  , the  vector 
= (x.  , y.)  is  the  position  of  the  j-th  sensor  relative  to  the  array  reference 
location,  f.  is  the  time  delay  (rounded  to  the  nearest  sample  interval)  re- 
quired to  time-align  a signal  from  the  look  direction  (as  seen  by  the  j-th  sen- 
sor) with  the  signal  as  seen  at  the  array  reference  location,  and  A (f)  *s  the 
Fou).  er  transform  of  the  convolution  filter  applied  to  the  j-th  sensor.  For 
the  amsteer  processor,  the  Fourier  transform  A^(f)  is  the  constant  value 
1/M,  where  M is  the  number  of  sensors  used  to  form  the  beam.  For  a fixed 
look  velocity  and  look  azimuth,  the  array  response  A(f  , k)  can  be  com- 
puted as  a function  of  azimuth  for  waves  propagating  at  the  look  velocity. 
Figure  II- 1 1 is  the  beamste  ^r  response  1 0 log ^ ^ | A ( f , k ) | computed  in 
this  way  for  the  six-sensor  noise  sample  fi  om  day  238  of  1970  when  the  steer 
azimuth  is  270°.  For  energy  arriving  from  the  west  (the  look  direction),  the 
beams teer  pattern  achieves  its  maximum  response  of  approximately  zero  dB. 
Unlike  the  beamsteer  array  response,  the  ABF  response  changes  with  time. 
Figures  11-12  through  11-15  display  the  adaptive  processor  response  at  1020, 
1025,  1030,  and  1035.  By  1035  a null  has  formed  at  the  60°  azimuth  of  the 
signal  presented  in  Figure  II- 10.  The  creation  of  the  null  is  an  automatic 
consequence  of  the  adaptive  algorithm,  which  attempts  to  minimize  the  filter 
output  power  subject  to  signal-preservation  constraints.  The  constraints  are 
reflected  at  an  azimuth  of  270°,  where  the  ABF  response  is  approximately 
zero  dB.  (Slight  variations  from  unity  response  at  the  look  azimuth  are  pos- 
sible because  the  sensor  time  delays  are  rounded  to  the  nearest  integer  sam- 
ple interval.  ) 

Signal  degradation  measurements  for  the  new  algorithm  again  use  the 
same  three  signals  contained  in  the  previously  processed  data  sample  frr-n 
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TIME-SHIFT  -FIND -SUM  BERM  PATTERN 
BERM  LOOK  VELOCITY  IS  3.6.  LOOK  RZIMUTH  270.0 
FREQUENCY  IS  0.04000  HZ.  PERIOD  25.0  SECONDS 


FIGURE  II- 11 

BEAMSTEER  RESPONSE  FOR  SIX-CHANNEL  NOISE  SAMPLE 

FROM  DAY  238  OF  1970 


ABF  BEAM  PATTERN  AT  238/10.20.00 
CONVERGENCE  RATE  IS  0, 0050 
BEAM  LOOK  VELOCITY  IS  3.6.  LOOK  AZIMUTH  270.0 
FREQUENCY  IS  0.04000  HZ.  PERIOD  25.0  SECONDS 


FIGURE  11-12 

ABF  RESPONSE  FOR  SIX-CHANNEL  N'  ISE  SAMPLE 
AT  1020  ON  DAY  238  OF  1970 


mm 


fIBF  BEAM  PATTERN  AT  238/  iO.  25. 00 
CONVERGENCE  RATE  IS  0.0050 
BEAM  LOOK  VELOCITY  IS  3.6.  LOOK  AZIMUTH  270.0 
FREQUENCY  IS  0.04000  HZ.  PERIOD  25.13  SECONDS 


FIGURE  II- I 3 

ABF  RESPONSE  FOR  SIX-CIIANNE  [.  NOISE  SAMPLE 
AT  102S  ON  DAY  238  OF  1870 


ALPfi 

ABF  BEAM  PATTERN  AT  230/10.30.00 
CONVERGENCE  RATE  !S  0.0030 
BEAM  LOOK  VELOCITY  IS  3.6.  LOOK  AZIMUTH  270.0 
FREQUENCY  IS  0. 04000  HZ . PER 1 00  25. 0 SECONDS 


FIGIJRK  11-14 


ABF  RESPONSE  FOR  SIX-CHANNEL  NOISE  SAMPLE 
AT  1030  ON  DAY  238  OF  1970 
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RBF  BERM  PATTERN  RT  238/10,35.00 
CONVERGENCE  RRTE  IS  0.0050 
BEAM  LOOK  VELOCITY  IS  3.6.  LOOK  AZIMUTH  270.0 
FREQUENCY  IS  0.04000  HZ,  PERIOD  25.0  SECONDS 


FIGURE  11-15 


ABF  RESPONSE  FOR  SIX  - CHANNE  L NOISE  SAMPLE 
AT  1035  ON  DAY  238  OF  1970 
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2000  to  2345  on  day  276  of  1971.  As  in  the  case  of  noise  reduction  measure- 
ments for  the  new  algorithm,  the  only  difference  in  processing  is  the  substi- 
tution of  3.  6 km/sec  for  the  3.  5 km/ sec  apparent  velocity  employed  with  the 
old  algorithm.  The  first  signal  examined  is  the  weak  signal  arriving  at  ap- 
proximately 2138  from  an  azimuth  between  300°  and  305°.  The  steer  direc- 
tion for  this  event  is  302.5°.  Table  II-  1 0 and  Figure  11-16  portray  signal 
degradation  as  a function  of  convergence  rate  for  this  event  with  the  new  ad- 
daptive  algorithm.  As  in  the  case  of  the  old  algorithm,  negative  signal  de- 
gradation is  obtained  at  the  low  convergence  rates.  The  amplification  of  the 
signal  relative  to  the  beamsteer  output  may  be  due  to  rounding  the  steer  de- 
lays to  the  nearest  sample  interval.  This  rounding  procedure  permits  slight 
deviations  from  the  unity-response  constraint  condition  for  the  look  direction. 
Figure  11-17,  which  is  a plot  of  the  ABF  beam  pattern  at  = 0.  005  shortly 
after  signal  onset,  confirms  the  measured  signal  degradation  of  -0.022  dB 
may  correspond  to  an  array  response  greater  than  0 dB  in  .’'e  look  direction. 
No  such  phenomenon  appears  in  the  time-shift-and-sum  beam  pattern  for  this 
array  configuration  and  look  direction  (Figure  11-18).  A comparison  of  the 
ABF  pattern  in  Figure  11-17  with  the  beamsteer  pattern  in  Figure  11-18  illus- 
trates the  capability  of  adaptive  processing  to  create  a superdirective  array 
response.  Using  the  noise  reduction  measurements  from  day  238  of  1970  for 
the  new  algorithm,  the  resulting  signal-to-noise  ratio  improvement  estimates 
as  a function  of  convergence  rate  appear  in  Table  II- 1 1 and  Figure  11-19.  The 
maximum  signal-to-noise  gain  of  1.47  dB  occurs  at  fi=  0.005.  This  gain  is 
about  0.2  dB  higher  than  the  maximum  gain  of  1.27  dB  for  this  weak  signal 
with  the  old  algorithm.  At  the  higher  convergence  rates,  the  new  algorithm 
outperforms  the  old  by  more  than  0.  5 dB.  A visual  examination  of  the  site 
vertical-component  traces  shows  poor  signal  similarity  across  the  array  for 
this  particular  weak  event.  This  fact  probably  explains  the  high  signal  degra- 
dation values  at  the  more  rapid  convergence  rates.  Later  in  this  section, 
when  utilization  of  the  full  array  and  a signal-enhancement  prefilter  permits 
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TABLE  11-10 

ADAPTIVE  FILTERING  SIGNAL  DEGRADATION  VERSUS 
CONVERGENCE  RATE  WITH  NEW  ALGORITHM  FOR 
A WEAK  SIGNAL  FROM  300°.  305° 

(ALPA,  6 SITES,  1971  DA Y 276  21.  37.  0 1 -2 1 . 45.  32) 


Convergence  Rate 

Noise  Reduction 

(M) 

(dB) 

0.  0010 

- 0.  233 

0.U015 

- 0.  26  3 

0.  0020 

- 0.  261 

0,0030 

- 0.  203 

0.  0040 

- 0.  116 

0,0050 

- 0.  022 

0.0070 

0,  161 

0.  0100 

0.408  ; 

0.0150 

0.  755 

0.  0200 

1, 042 

0.  0300 

1.490 

0.0400 

1.  829 

0,0500 

2.  101 

0.  0700 

2.  531 

0.1000 

3.  023 

0,  1500 

3.639 

0.  2000 

4,  1 15 

0.  3000 

4.  857 

0.  4000 

5.457 

0.  5000 

5.976 

Signal  Degra 


Convergence  Rate  /i 
FIGURE  11-16 

SIGNAL  DEGRADATION  AS  A FUNCTION  OF  CONVERGENCE 
RATE  WITH  NEW  ALGORITHM  FOR  A WEAK  SIGNAL  FROM 
300°-305°  (ALPA,  6 SITES,  1971  DAY  276  21.  37.  01-21. 45.  32) 
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ABF  BERM  PATTERN  RT  276/21.39.00 
CONVERGENCE  RATE  IS  0.0050 
BERM  LOOK  VELOCITY  IS  3.6,  LOOK  AZIMUTH  302.5 
FREQUENCY  IS  Q.0'4000  HZ,  PERIOD  25.0  SECONDS 


FIGURE  P-  17 

SIX - CHANNEL  ABF  RESPONSE  FOR  WEAK  EVENT  FROM  300°-305° 

AT  2139  ON  DAY  2 7f>  OF  1971 
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TIME-SHIFT -FIND-SUM  BEAM  PATTERN 
BEAM  LOOK  VELOCITY  IS  3.6,  LOOK  AZIMUTH  302.5 
FREQUENCY  IS  0.04000  HZ,  PERIOD  25.0  SECONDS 


TABLE  11-11 

ADAPTIVE  FILTERING  SIGNAL-TO-NOISE  GAIN  VERSUS 
CONVERGENCE  RATE  WITH  NEW  ALGORITHM  FOR 
A WEAK  SIGNAL  FROM  300°  -305° 

(ALPA,  6 SITES,  1971  DAY  276  21.37.01-21.45.  32, 
USING  NOISE  FROM  1970  DAY  238  0800-1145) 


Convergence  Rate 
00 

Signal-to-Noise  Gain 
(dB) 

0.0010 

0.  874 

0.0015 

1.  088 

0.0020 

1. 230 

0.  0030 

1.  386 

0.0040 

1.450 

0.0050 

1.469 

0.0070 

1.448 

0.  0100 

1.  364 

0.0150 

1.  205 

0. 0200 

1. 064 

0.0300 

0.  862 

0.0400 

0.  742 

0.0500 

0.674 

0.0700 

0.615 

0.1000 

0.  605 

; o.  1 500 

0.  653 

0.2000 

0.  714 

0. 3000 

0.  787 

0.4000 

0.  736 

0. 5000 

0.  551 

Convergence  Rate  M 


FIGURE  II-  1 9 

SIGNAL -TO -NOISE  GAIN  AS  A FUNCTION  OF  CONVERGENCE 
RATE  WITH  NEW  ALGORITHM  FOR  A WEAK  SIGNAL  FROM 
300°-305°  (ALPA,  6 SITES,  1971  DAY  276  21.37.01-21.45.32 
USING  NOISE  FROM  1970  DAY  238  0800-1145) 


detection  of  the  New  Britain  event,  the  processing  results  will  shed  further 
light  on  this  possibility.  Figure  11-20  is  a display  of  the  weak  event  from 
300°-305°  for  the  beamsteer  output  and  adaptive  filter  output  at  4 = 0.005 
and  4 = 0.5.  On  the  middle  trace  (where  4=  0.005),  the  signal  is  approxi- 
mately as  strong  as  on  the  beamsteer  output.  In  contrast,  all  but  the  closing 
section  of  the  event  is  greatly  attenuated  at  4=  0.5  (the  bottom  trace).  Even 
on  the  bottom  trace,  however,  the  signal  remnants  are  conspicuous  enough 
to  suggest  the  presence  of  the  weak  event.  In  fact,  since  the  noise  surround- 
ing the  event  is  considerably  weaker  on  the  bottom  trace,  the  signal  detect- 
ability for  the  two  ABF  traces  is  closer  than  it  seems  at  first  glance.  This 
fact  is  important  when  considering  high  convergence  rates  for  detection  pro- 
cessing. 

Several  events  in  the  data  sample  from  day  276  of  1971  come  from 
azimuths  different  from  the  302.  5 steer  azimuth  of  the  time-shift-and- sum 
and  adaptive  beams  which  provided  the  signal  degradation  measurements  for 
the  weak  event  from  300  -305  . The  ability  to  suppress  such  off-azimuth 
events  is  an  important  measure  of  ABF  performance,  particularly  in  detec- 
tion processing.  If  the  AJ3F  traces  are  significantly  weaker  than  the  corres- 
ponding beamsteer  traces  for  off-azimuth  events,  adaptive  processing  can 
reduce  the  number  of  false  alarms  from  off-azimuth  events  and  provide  better 
event  arrival-direction  estimates.  A further  potential  advantage  is  the  pos- 
sibility of  separating  interfering  events.  To  demonstrate  interfering -event 
separation  capability,  however,  it  must  be  shown  that  adaptive  beamforming 
will  not  manipulate  interfering  events  in  such  a way  that  they  cancel  each 
other.  Table  11-12  presents  the  power  reduccion  relative  to  beamsteering  at 
4 = 0.  005  and  4 = 0.  5 over  the  time  gates  indicated  fc  r five  events  with  azi- 
muths away  from  the  302.  5 steer  direction.  In  the  case  of  the  Panama  and 
Kamchatka  events,  the  latter  event  of  each  pair  is  the  stronger.  If  signal 
degradation  with  the  new  algorithm  can  be  kept  within  reasonable  bounds,  the 
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dramatic  event-suppression  values  at  4 = 0.  5 suggest  that  operation  at  high 
convergence  rates  may  result  in  greatly  improved  detection  capability.  Fig- 
ures 11-21  through  11-23  illustrate  the  ABF  performance  for  events  close  to 
the  302.  5°  steer  direction.  The  source  direction  is  27  3°  for  the  two  Kam- 
chatka events  and  237°  for  the  Andreanof  Islands  event.  At  the  convergence 
parameter  value  4 = 0.5,  all  three  events  are  virtually  obliterated.  At  the 
slower  convergence  rate  in  the  middle  traces,  however,  the  attenuation  for 
the  two  Kamchatka  events  (30°  from  the  steer  direction)  is  less  than  for  the 
Andreanof  Islands  event  (65°  from  the  steer  direction).  To  illustrate  the  way 
in  which  adaptive  beamforming  at'.rr  uates  events  close  to  the  steer  direction, 
Figures  11-24  through  11-26  present  ABF  beam  patterns  at  4 = 0.005  for  the 
Andreanof  Islands  event  at  2325,  2330,  and  2335.  The  adaptive  algorithm 
bends  the  main  lobe  of  the  beam  pattern  array  from  the  off-azimuth  evert: 
the  maximum  response  (which  is  greater  than  unity)  occurs  at  an  azimuth  of 
approximately  315°  instead  of  302.5°.  The  corresponding  beamsteer  response 
for  the  302.  5°  steer  direction  was  presented  earlier  in  Figure  11-18. 

The  last  two  signals  used  for  signal  degradation  measurements  with 
the  new  algorithm  arrive  from  Kamchatka  at  an  azimuth  of  273  . The  first 
of  these  two  Kamchatka  events  reached  ALPA  at  approximately  2108.  Table 
11-13  and  Figure  11-27  depict  signal  degradation  as  a function  of  convergence 
rate  for  this  event  with  the  new  adaptive  algorithm.  The  .measured  degrada- 
tion values  with  the  new  algorithm  are  much  lower  than  with  the  old  algorithm 
(Table  11-5  and  Figure  11-5).  When  the  signal  similarity  across  the  array  is 
good  (as  in  this  case),  the  differences  between  the  old  algorithm  and  the  new 
algorithm  emerge:  the  strength  and  similarity  of  the  signal  result  in  a much 
slower  adaptation  rate  for  the  new  algorithm;  the  fact  that  this  event  has  great- 
er signal  similarity  across  the  array  than  the  weak  event  from  300°-305°, 
moreover,  makes  it  more  resistant  to  attenuation  by  the  adaptive  algorithm. 
Figure  11-28  is  a Calcomp  plot  of  the  beamsteer  output  and  adaptive  filter 
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SIX-CHANNEL  BEAMS  FOR  MAGNITUDE  4.  3 
NDREANOF  ISLANDS  EVENT  (DAY  276  1971, 
STEER  DIRECTION  302.  5°) 
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RBF  BERM  PATTERN  RT  276/23. 25. CD 
CONVERGENCE  RATE  IS  0.005G’ 

BERM  LOOK  VELOCITY  IS  3.6.  LOOK  AZIMUTH  312.5 
FREQUENCY  IS  0.0<t000  HZ.  PERIOD  25. C SECONDS 


FIGURK  11-24 

SIX-CIIANNEI,  ABF  RESPONSE  FOR  AN  DR  HANOI'  IS  I WHS 

EVENT  AT  2325  ON  DAY  276  OF  1971 


ALPA 

ABF  BEAM  PATTERN  AT  276/23.30.00 
CONVERGENCE  RATE  IS  0.0050 
BEAM  LOOK  VELOCITY  IS  3.6,  LOOK  AZIMUTH  302.5 
FREQUENCY  IS  O.OAOOO  HZ,  PERIOD  25.0  SECONDS 
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FIGURE  11-25 

SIX-CHANNEL  ABF  RESPONSE  FOR  ANDREANOF  ISLANDS 
EVENT  AT  2330  ON  DAY  276  OF  1971 
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ALPA 

ABF  BEAM  PATTERN  AT  276/23.  35. 00 
CONVERGENCE  RATE  IS  0.0050 
BEAM  LOOK  VELOCITY  IS  3.6,  LOOK  AZIMUTH  302.5 
FREQUENCY  IS  0.04000  HZ,  PERIOD  25.0  SECONDS 


FIGURE  11-26 

SIX-CHANNEL  ABF  RESPONSE  FOR  ANDREANOF  ISLANDS 
EVENT  AT  2335  ON  DAY  276  OF  1971 
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TABLE  11-13 


ADAPTIVE  FILTERING  SIGNAL  DEGRADATION  VERSUS 
CONVERGENCE  RATE  WITH  NEW  ALGORITHM  FOR 
A STRONG  SIGNAL  FROM  KAMCHATKA 
(ALPA,  6 SITES,  1971  DAY  276  21 . 07.  1 8-2 1 . 1 1 . 33) 


Convergence  Rate 
(M) 

Signal  Degradation 
<dB) 

0.0010 

0.  031 

0.0015 

0.  047 

0.  0020 

0.062 

0.0030 

0.092 

0.0040 

0.  118 

0.0050 

0.  141 

0.0070 

0.  177 

0.0100 

0.  214 

0.0150 

0.  246 

0.  0200 

0.  262 

0.0300 

0.  281 

0.0400 

0.  304 

0.0500 

0.  336 

0.0700 

0.418 

0.  1000 

0.  551 

0.  1500 

0.  722 

0.2000 

0.  814 

0.  3000 

0.  853 

0.4000 

! 0.819 

0.  5000 

0.  777 
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FIGURE  11-27 

SIGNAL  DEGRADATION  AS  A FUNCTION  OF  CONVERGENCE 
RATE  WITH  NEW  ALGORITHM  FOR  A STRONG  SIGNAL  FROM 
KAMCHATKA  (ALPA,  6 SITES,  1 97 1 DAY  276  21.  07.  1 8- 21 . 1 1 . 33) 
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outputs  at  M = 0.  005  and  M=  0.  5 for  this  event  with  the  new  adaptive  algor- 
ithm. The  differences  between  the  beamsteer  output  and  the  ABF  output  at 
the  slower  convergence  rate  are  slight.  Even  at  the  more  rapid  rate,  the 
ABF  output  is  only  slightly  attenuated,  although  signal  distortion  is  clearly 
evident.  Table  11-14  and  Figure  11-29  present  the  estimated  signal-to-noise 
ratio  improvement  as  a function  of  convergence  rate  for  this  event  using 
noise  reduction  measurements  from  day  238  of  1970.  The  greatest  gain  re- 
lative to  beamsteering  occurs  at  H = 0.5,  where  the  estimated  improvement 
is  5.75  dB.  Despite  the  signal  distortion  in  the  bottom  trace  of  Figure  11-28, 
the  5.  75  dB  gain  estimate  is  definitely  relevant  for  detection  processing. 

The  second  Kamchatka  event  reached  ALPA  at  approximately  2207. 
Signal  degradation  as  a function  of  convergence  rate  for  this  event  with  the 
new  adaptive  algorithm  appears  in  Table  11-15  and  Figure  11-30.  Not  only  is 
the  degradation  slight  for  this  event,  but,  indeed,  even  negative  at  the  high- 
est convergence  rates.  The  amplification  of  the  second  Kamchatka  event  at 
M = 0.  5 is  illustrated  by  the  bottom  trace  in  Figure  11-31.  In  contrast,  the 
adaptive-filter  output  at  M = 0.  005  (the  middle  trace)  is  almost  a carbon  copy 
of  the  beamsteer  output  in  the  top  trace.  Table  11-16  and  Figure  11-32  give 
the  estimated  signal-to-noise  ratio  improvement  as  a function  of  convergence 
rate  for  this  event  with  the  new  adaptive  algorithm  (using  noise  reduction 
measurements  from  day  238  of  1970).  Maximum  estimated  gain  of  7.5  dB  is 
realized  at  the  highest  convergence  rate.  The  gain  estimates  for  the  new 
algorithm  are  distinctly  higher  than  for  the  old  algorithm  (see  Table  II- 8 and 
Figure  II-8).  The  sharply  improved  performance  is  due  to  the  greatly  re- 
duced effective  adaptation  rate  of  the  new  algorithm  and  the  good  signal  sim- 
ilarity of  the  second  Kamchatka  event  between  sites. 

Four  of  the  events  in  the  data  sample  from  day  276  of  1971  arrive  at 
azimuths  different  from  the  27  3°  steer  azimuth  used  for  the  two  Kamchatka 
events.  Table  11-17  gives  the  power  reduction  relative  to  beamsteering  at 
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TABLE  11-14 


ADAPTIVE  FILTERING  SIGNAL-TO-NOISE  GAIN  VERSUS 
CONVERGENCE  RATE  WITH  NEW  ALGORITHM  FOR 
A STRONG  SIGNAL  FROM  KAMCHATKA 
(ALPA,  6 SITES,  1971  DAY  276  21.  07.  18-21.  11.  33, 
USING  NOISE  FROM  1970  DAY  238  0800-1145) 


Convergence  Rate 
(M) 

0. 0010 
0.  0015 
0. 0020 
0.  0030 
0.  0040 
0. 0050 
0.  0070 
0. 0100 
0. 0150 
0. 0200 
0.  0300 
0. 0400 
0.  0500 
0.  0700 
0. 1000 
0. 1500 
0.  2000 
0.  3000 
0.4000 
0.  5000 


Signal-to-Noise  Gain 
(dB) 

0.610 

0.  778 

0.  907 
1.091 
1. 216 

1.  306 
1.432 

1.  558 
1. 714 
L 844 

2.  071 

2.  267 
2.439 

2.  728 

3.  077 

3.  570 

4.  015 

4.  791 

5.  374 

5.  750 
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FIGURE  11-29 


SIGNAL- TO -NOISE  GAIN  AS  A FUNCTION  OF  CONVERGENCE 
RATE  WITH  NEW  ALGORITHM  FOR  A STRONG  SIGNAL  FROM 
KAMCHATKA  (AL PA,  6 SITES,  1971  DAY  276  21.07.18-21.11.33 
USING  NOISE  FROM  1970  DAY  238  0800-1  145) 
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TABLE  II-15 

ADAPTIVE  FILTERING  SIGNAL  DEGRADATION  VERSUS 
CONVERGENCE  RATE  WITH  NEW  ALGORITHM  FOR 
A VERY  STRONG  SIGNAL  FROM  KAMCHATKA 
(ALPA,  6 SITES,  1971  DA  Y 276  22.  07.  02-22.  1 1.  17) 


Convergence  Rate 


Signal  Degradation 


0. 0010 
0.0015 
0. 0020 
0.0030 
0.0040 
0.0050 
0.0070 
0.0100 
0.0150 
0.0200 
0.  0300 
0.0400 
0.0500 
0.0700 
0. 1000 
0.  1500 
0.2000 
0. 3000 


0.4000 

0.5000 


0.066 
0.  093 
0.  115 
0.  145 
0.  164 
0.  177 
0.  197 
0.  223 
0.255 
0.276 
0.  289 
0.279 
0.255 
0.  185 
0.055 

- 0.  173 

- 0.  379 

- 0.688 

- 0.885 

- 1.018 
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SIGNAL  DEGRADATION  AS  A FUNCTION  OF  CONVERGENCE 
RATE  WITH  NEW  ALGORITHM  FOR  A VERY  STRONG  SIGNAL 
FROM  KAMCHATKA  (ALPA,  6 SITES,  1971  DAY  276  22.07.02-22.  11.17) 
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TABLE  11-16 

ADAPTIVE  FILTERING  SIGNAI  -TO -NOISE  GAIN  VERSUS 
CONVERGENCE  RATE  WITH  NEW  ALGORITHM  FOR 
A VERY  STRONG  SIGNAL  FROM  KAMCHATKA 
(ALPA,  6 SITES,  1971  DAY  276  22.  07.02-22.  11.  17, 
USING  NOISE  FROM  1970  DAY  238  0800-1145) 


Convergence  Rate 
(M) 

Signal-to-Noise  Gain 
(dB) 

0.0010 

0.575 

0.0015 

0.  732 

0.0020 

0.  854 

0.0030 

1.038 

0.0040 

1.  170 

0.0050 

1.  270 

0.0070 

1.  412 

0.0100 

1.549 

0.0150 

1.  705 

0.0200 

1. 830 

0.0300 

2.063 

0.0400 

2.  292 

0.0500 

2.  520 

0.  0700 

2.  961 

0,1000 

3.  573 

0. 1500 

4,465 

0.2000 

5.  208 

0. 3000 

6.  332 

0,4000 

7.  078 

0.  5000  , 

7.  545 

7.  545 


Signal  -to  -No 
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FIGURE  11-32 

SIGNAL -TO-NOISE  GAIN  AS  A FUNCTION  OF  CONVERGENCE 
RATE  WITH  NEW  ALGORITHM  FOR  A VERY  STRONG  SIGNAL  FROM 
KAMCHATKA  (ALPA,  6 SITES,  1971  DAY  276  22.07.02-22.11.17, 
USING  NOISE  FROM  1970  DAY  238  0800-1146) 


1^-  0.  005  and  M = 0.  5 for  these  four  events.  Except  in  the  case  of  the  weak 
event  from  300  -305  , these  events  are  attenuated  by  more  than  12  dB  at 
M = 0.  5.  Figure  11-33  shows  the  beamsteer  output  and  the  adaptive  filter 
output  at  M = 0.  005  and  M = 0.  5 for  the  first  Panama  event.  Figure  11-34  is 
a plot  of  the  corresponding  beams  for  the  other  Panama  event.  The  associated 
adaptive  beam  patterns  in  Figure  11-35  through  11-37  reflect  the  strong  sup- 
pression of  this  event  at  (1=  0.005.  In  the  case  of  this  event  (160°  away  from 
the  steer  azimuth),  the  beam  patterns  at  the  times  2120,  2124,  and  2128  il- 
lustrate the  ability  of  adaptive  processing  to  form  deep  nulls  in  the  direction 
of  off-azimuth  energy  sources.  For  purposes  of  comparison,  the  six-channel 
time-shift-and-sum  pattern  for  a 273°  steer  azimuth  appears  in  Figure  11-38. 

The  weak  event  from  300°-305°  is  portrayed  in  Figure  11-39.  At  an  azimu- 
. o 

thal  separation  of  30  from  the  steer  direction,  the  event  is  only  slightly 
diminished  relative  to  the  beamsteer  output  in  the  adaptive  filter  output  for 
M = 0.  005  (the  middle  trace).  Curiously  enough,  the  ABF  output  for  this 
event  at  H = 0.  5 is  stronger  near  the  beginning  than  near  the  end.  The  con- 
verse was  true  for  the  corresponding  adaptive  beam  at  a steer  azimuth  of 
302.5  (Figure  11-20).  Perhaps  the  arrival  azimuth  shifts  significantly  dur- 
ing this  event.  Another  possioility  is  that  two  distinct  earthquakes  have  ar- 
rived simultaneously  at  ALP  A in  the  time  period  shown  in  Figure  11-39.  A 
situation  such  as  this  might  explain  the  lack  of  signal  similarity  for  this  event. 
Figure  11-40  shows  the  beam  outputs  for  the  Andreanof  Islands  event.  Here, 
where  the  steer  azimuth  is  only  36  from  the  event  arrival  azimuth,  the  sig- 
nal is  much  stronger  than  in  the  beams  for  the  302,  5°  steer  azimuth  (Figure 
11-23).  where  the  azimuthal  separation  is  65°.  Even  at  M = 0.5,  the  signal 
is  not  completely  eliminated. 

Ihe  results  presented  so  far  in  this  section  demonstrate  the  clear 
superiority  of  the  new  adaptive  algorithm.  This  superiority  manifests  itself 
primarily  in  the  form  of  reduced  signal  degradation.  The  differences  are 
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FIGURE  II  - 3 3 

CHANNEL  BEAMS  FOR  MAGNITUDE  3.6  PANAMA  EVENT 
(DAY  276  1971,  STEER  DIRECTION  273  ) 


CHANNEL  BEAMS  FOR  MAGNITUDE  4.  7 PANAMA  EVENT 
(DAY  276  1971,  STEER  DIRECTION  273°) 
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RBF  BERM  PRTTERN  RT  276/21.20.00 
CONVERGENCE  RRTE  15  0.0050 
BERM  LOOK  VELOCITY  15  3.6,  LOOK  RZ1MUTH  273.0 
l-REQUENCY  15  0.0*1000  HZ,  PERIOD  25.0  SECONDS 


FIGURE  11-35 

SIX-CHANNEL  ABF  RESPONSE  FOR  MAGNITUDE 
4.  7 PANAMA  EVENT  AT  2120  ON  DAY  276  OF  1971 
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RBF  BERM  PRTTERN  RT  276/21.24.00 
CONVERGENCE  RRTE  15  0.0050 
BERM  LOOK  VELOCITY  15  3.6,  LOOK  RZIMUTH  273.0 
FREQUENCY  15  0.04000  HZ,  PERIOD  25.0  SECONDS 


FIGURE  11-36 

SIX-CHANNEL  ABF  RESPONSE  FOR  MAGNITUDE 
4.  7 PANAMA  EVENT  AT  Z1Z4  ON  DAY  Z76  OF  1971 
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RBF  BERM  PATTERN  RT  276/21.28.00 
CONVERGENCE  RATE  15  0.0050 
BERM  LOOK  VELOCITY  IS  3.6,  LOOK  AZIMUTH  273.0 
FREQUENCY  IS  0.04000  HZ,  PERIOD  25.0  SECONDS 


FIGURE  11-37 

SIX-CHANNEL  ABF  RESPONSE  FOR  MAGNITUDE 
4.  7 PANAMA  EVENT  AT  2128  ON  DAY  276  OF  1971 
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T I ME-SH I FT -RND-SUM  BERM  PATTERN 
BERM  LOOK  VELOCITY  15  3.6,  LOOK  AZIMUTH  273.0 
FREQUENCY  IS  0.0*4000  HZ,  PERIOD  25.0  SECONDS 


FIGURE  11-38 

BEAMSTEER  RESPONSE  FOR  SIX-CHANNEL  SAMPLE 
FROM  DAY  276  OF  1971 


SIX-CHANNEL  BEAMS  FOR  WEAK  EVENT  FROM  300  -305 
(DAY  276  1971,  STEER  DIRECTION  273°) 


-CHANNEL  BEAMS  FOR  ANDREANOF  ISLANDS  EVENT 
(DAY  276  1971,  STEER  DIRECTION  273°) 
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particularly  dramatic  for  the  two  Kamchatka  events,  which  are  strong  events 
with  good  signal  similarity.  For  the  weak  event  from  300° -305  , the  im- 
provement in  signal-to-noise  gain  with  the  new  algorithm  is  only  slight.  The 
probable  explanation  is  the  poor  signal  similarity  of  this  event.  With  the  new 
adaptive  algorithm,  the  best  signal-to-noise  gains  occur  at  the  most  rapid 
convergence  rate  in  the  case  of  the  two  Kamchatka  events.  Even  for  the  weak 
event  from  300°-305°,  the  gain  at  0.  5 is  only  one  dB  less  than  the  gain  at 
H = o.  005.  In  detection  processing,  where  signal  distortion  is  a relatively 
minor  consideration,  these  results  suggest  that  significantly  improved  detec- 
tion capability  for  events  with  good  signal  similarity  is  possible  at  high  con- 
vergence rates.  The  probable  price  to  be  paid  for  this  improved  detection 
performance  is  a somewhat  reduced  capability  with  signals  such  as  the  weak 
event  from  300°-305°,  where  the  signal  similarity  is  poor.  Another  effect 
of  operating  at  rapid  convergence  rates  is  a narrowing  of  the  adaptive  pro- 
cessor beamwidth,  as  indicated  by  the  greater  suppression  of  off-azimuth 
events  (see  Tables  11-12  and  11-17).  In  one  sense,  the  effect  is  beneficial 
since  the  direction  of  incoming  events  can  be  more  accurately  determined. 

In  an  operational  surveillance  mode,  however,  the  improved  resolution  at 
high  convergence  rates  raises  the  number  of  beams  required  to  scan  all  azi- 
muths without  gaps  in  coverage.  With  reduced  background  noise  levels  at 
fast  convergence  rates,  the  probability  of  lime -overlapped  events  is  greater 
than  the  0.  05  - 0.  10  figure  normally  given  for  ALP  A.  The  off-azimuth  event 
suppression  results  suggest,  at  least,  that  rapid  adaptation  can  reduce  the 
number  of  times  that  one  event  masks  another. 


D.  ADAPTIVE  PROCESSING  WITH  THE  FULL  ARRAY 

The  purpose  of  this  subsection  is  to  evaluate  adaptive  processing 
which  uses  all  available  sensors  at  ALPA  as  input  to  the  time-shift-and-sum 
and  adaptive  beamformers.  Of  particular  interest  are  any  differences 


J 
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between  the  full-array  results  and  the  results  obtained  with  closely-spaced 
six-channel  arrays  such  as  the  ones  employed  in  Subsection  C.  Since  the 
new  adaptive  algorithm  outperforms  the  old  adaptive  algorithm,  all  adaptive 
beams  in  this  subsection  and  the  remainder  of  this  report  are  formed  with 
the  new  adaptive  algorithm. 

During  1970,  a maximum  of  nine  sensors  were  operational  at  ALP  A. 

In  order  to  compare  full-array  results  with  six-channel  results,  a noise  sam- 
ple from  0415  to  0815  on  day  335  of  1972  was  selected  to  replace  the  noise 
sample  from  day  238  of  1970.  For  the  six-channel  noise  reduction  measure- 
ments, the  vertical  components  from  sites  1,  2,  3,  6,  8,  and  9 are  input  to 
the  time-shift-and-sum  and  adaptive  beams.  The  steer  direction  corresponds 
to  an  azimuth  of  270°  and  an  apparent  velocity  of  3.  6 km/ sec.  The  PDE, 

LASA,  and  NORSAR  bulletins  list  no  events  arriving  at  ALPA  between  0415 
and  0815  on  30  November  1972.  A Fisher  detector  scan  indicates,  however, 
that  a small  signal  reached  ALPA  at  0458  from  an  azimuth  of  approximately 
12°.  The  energy  from  this  event  should  affect  the  measured  noise  reduction 
values  only  minimally.  Table  11-18  and  Figure  11-41  present  the  noise  reduc- 
tion for  six  channels  as  a function  of  the  convergence  parameter  M . The 
noise  reduction  values  for  this  sample  are  lower  than  for  the  sample  from  day 
238  of  1970,  particularly  at  the  higher  convergence  rates. 

Signal  degradation  measurements  from  day  276  of  1972  are  borrowed 
from  Subsection  C to  obtain  estimates  of  the  signal-to-noise  gains  which 
would  have  occurred  if  a signal  had  arrived  at  ALPA  between  0415  and  0815 
on  day  335  of  1972.  Signal-to-noise  ratio  improvement  relative  to  beamsteer- 
ing is  depicted  as  a function  of  convergence  rate  in  Table  11-19  and  Figure  11-42 
for  the  weak  event  from  300°-305°,  in  Table  11-20  and  Figure  11-43  for  the 
strong  signal  from  Kamchatka,  and  in  Table  11-21  and  Figure  11-44  for  the 
very  strong  signal  from  Kamchatka.  These  gain  estimates  provide  a basis 
for  comparison  with  the  full-array  processing  results  which  appear  in  the  re- 
mainder of  this  subsection. 
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Convergence  Rate 

(M) 


Noise  Reduction 
(dB) 


0.0010 
0.  0015 
0.  0020 
0.  0030 
0.  C040 
0.  0050 
0.  0070 
0.  0100 
0.0150 
0.  0200 
0. 0300 
0.0400 
0.0500 
0.  0700 
0.  1000 
0.  1500 
0, 2000 
0.  3000 
0.4000 
0,5000 


0.677 
0.  852 

0.  991 
1. 208 

1.  372 
1.  502 
1.699 

1.  901 

2.  1 18 
2.264 
2.465 
2.61  1 

2.  730 
2.  928 

3.  172 

3.  507 
3.792 

4.  265 
4,  622 
4.  856 


TABLE  11-19 


ADAPTIVE  FILTERING  SIGNAL-TO -NOISE  GAIN 
VERSUS  CONVERGENCE  RATE  FOR  A WEAK 
SIGNAL  FROM  300°-305°  (ALPA,  6 SITES, 

1971  DAY  276  21.  37.  01-  21.45.  32,  USING  NOISE 
^ROM  1972  DAY  335  0415-08151 


Convergence  Rate 

(Ji> 

Signal-to-Noise  Gain 
(dB) 

0.  0010 

0.  910 

0.  0015 

1.  M5 

0.  0020 

1.  252 

0.  0030 

1.411 

0.  0040 

1.  488 

0.  0050 

1.  524 

0.  0070 

1.538 

0.  0100 

1.493 

0.  0150 

1.  363 

0.  0200 

1.  222 

0.  0300 

0.  975 

0. 0400 

0.  782 

0.  0500 

0.629 

0.  0700 

0.  397 

0.  1000 

0.  149 

0.  1 500 

-0.  1 32 

0.  2000 

-0.  323 

0.  3000 

-0.  588 

0.  4000 

-0.  835 

0.  5000 

-1. 120 
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Convergence  Rate  M 


FIGURE  11-42 

SIGNAL -TO-NOISE  GAIN  AS  A FUNCTION  OF  CONVERGENCE 
RATE  FOR  A WEAK  SIGNAL  FROM  300°-305°  (ALPA,  6 SITES,  1971 
DAY  276  21.  37.01-21.  45.  32,  USING  NOISE  FROM  1972  DAY  335 

0415-0815) 


TABLE  11-20 

ADAPTIVE  FILTERING  SIGNAL-TO-NOISE 
GAIN  VERSUS  CONVERGENCE  RATE  FOR 
A STRONG  SIGNAL  FROM  KAMCHATKA  (ALPA, 
6 SITES,  1971  DAY  276  2 1.  07.  1 8-2  1 . 1 1 . 33, 
USING  NOISE  FROM  1972  DAY  335  0415-0815) 


Convergence  Rate 


Signal-to-Noise  Gain 


0.  0010 

0.646 

0.0015 

0.  805 

0.  0020 

0.  929 

0.  0030 

1.  116 

0.  0040  ; 

1.254 

0.  0050  j 

1.  361 

0.  0070 

1.  522 

0.  0100 

1.687 

0.  0150 

1.  872 

0.  0200 

2.  002 

0.  0300 

2.  1 84 

0.  0400 

2.  307 

0. 0500 

2.  394 

0.  0700 

2.  510 

0. 1000 

2.621 

0.  1500 

2.  785 

0.  2000 

2.  978 

0. 3000 

3.412 

0. 4000 

3.803 

0.  5000 

4.079 

1 

-)  |i{  :jj;  •••; 

1 1 ■ 1 1 t+H  'tt 
J 


I’M  * 

Mi  'it 
* t • j Ml 
1 * • Mi* 


*|M 

ft**  ' [I’ 

: 


1 . 

I 

. 

1 

I ! it  : 

I I I 

lit1 


jj 

. i -■ 

1 j t 

, i i 

II 

‘ T 

Convergence  Rate  /j, 

FIGURE  11-43 

SIGNAL -TO -NOISE  GAIN  AS  A FUNCTION  OF  CONVERGENCE 
RATE  FOR  A STRONG  SIGNAL  FROM  KAMCHATKA  (ALPA,  6 SITES 
1971  DAY  276  21.07.  18-21.  1 1.  33,  USING  NOISE  FROM  1972  DAY  335 

0415-0815) 
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TABLE  11-21 

ADAPTIVE  FILTERING  SIGNAL-TO-NOISE  GAIN 
VERSUS  CONVERGENCE  RATE  FOR  A VERY 
STRONG  SIGNAL  FROM  KAMCHATKA  (ALPA, 

6 SITES,  1971  DAY  276  22.07.  02-22.11.17, 
USING  NOISE  FROM  1972  DAY  335  0415-0815) 


Convergence  Rate 
(M) 

Signal-to-Noise  Gain 
(dB) 

0.  0010 

0.611 

0.  0015 

0.  759 

0.  0020 

0.876 

0.  0030 

1.063 

0. 0040 

1.208 

0. 0050 

1.  325 

0. 0070 

1.  502 

0. 0 100 

1.678 

0.  0 150 

1.  863 

0. 0200 

1.  988 

0.  0300 

2.  176 

0. 0400 

2.  332 

0.  0500 

2,  475 

0.  0700 

2.  743 

0.  1000 

3.117 

0. 1500 

3.  680 

0.  2000 

4.  1 71 

0. 3000 

4.  953 

0.  4000 

5.  507 

0. 5000 

5.  874 
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FIGURE  11-44 


SIGNAL  - TO -NOISE  GAIN  AS  A FUNCTION  OF  CONVERGENCE 
RATE  FOR  A VERY  STRONG  SIGNAL  FROM  KAMCHATKA  (ALPA,  6 SITES 
1971  DAY  276  22.  17.02-22.  11.  17,  USING  NOISE  FROM  1972  DAY  335 

0415-0815) 


The  full-array  noise  reduction  measurements  utilize  the  same  noise 
sample  from  0415  to  0815  on  da/  335  of  1972  as  was  used  for  the  six-channel 
results.  The  only  difference  in  processing  is  that  fifteen  sites  (all  but  sites 
7,  16,  18,  and  19)  form  the  channel  inputs  to  the  time-shift-and-sum  and 
adaptive  beamformers.  Table  11-22  and  Figure  11-45  portray  the  noise  re- 
duction for  fifteen  channels  as  a function  of  convergence  rate. 

Signal  degradation  measurements  for  the  full  array  again  employ  the 
same  three  signals  in  the  previously  processed  data  sample  from  2000  to 
2345  on  day  276  of  1971.  All  processing  procedures  for  this  sample  are  iden- 
tical to  those  which  yielded  the  signal  degradation  measurements  in  Subsec- 
tion C for  the  new  adaptive  algorithm  except  that  sixteen  sites  (all  but  sites 
7,  '.l,  and  18)  constitute  the  input  to  the  beamformers.  The  first  signal  con- 
sidered is  the  weak  signal  from  300°- 305°.  The  look  azimuth  for  this  event 
is  302.  5°.  Table  11-23  and  Figure  11-46  depict  signal  degradation  as  a func- 
tion of  convergence  rate  for  this  event.  The  resulting  signal-to -noise  gain 
estimates  corresponding  to  the  noise  sample  from  day  335  of  1972  appear  in 
Table  11-24  and  Figure  11-47.  Maximum  signal-to-noise  ratio  improvement 
is  1.  15  dB  at  M = 0.  5.  There  is  also  a local  maximum  at  H = 0.  004,  where 
the  estimated  gain  is  0.26  dB.  The  shape  of  the  gain  curve  is  distinctly  dif- 
ferent from  the  corresponding  six-channel  gain  curve  in  Figure  11-42,  where 
maximum  improvement  relative  to  beamsteering  is  1.5  dB  at  U=  0.  007.  Fig- 
ure 11-48  is  a display  of  the  weak  event  from  300°-305°  for  the  beamsteer  out- 
put and  ABF  output  at  H = 0.  005  and  M=  0.  5.  Note  that  the  peak-to-peak  am- 
plitude at  the  end  of  the  event  in  the. bottom  trace  is  comparable  to  the  peak- 
to-peak  amplitude  of  the  event  in  the  beamsteer  output.  The  noise  level  be- 
fore and  after  the  event,  in  contrast,  is  greatly  attenuated  at  /U  = 0.  5.  Thus, 
if  signal  detection  were  based  on  peak-to-peak  amplitude  measurements,  the 
best  detection  results  for  this  event  would  occur  at  the  highest  convergence 
rate. 
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TABLE  11-22 

ADAPTIVE  FILTERING  NOISE  REDUCTION 
VERSUS  CONVERGENCE  RATE 
(ALPA,  15  SITES,  1972  DAY  335  0415-0815) 


Convergence  Rate 

Noise  Reduction 

(Mt 

(dD) 

0. 0010 

0.  357 

0. 0015 

0.482 

0. 0020 

0.587 

0.  0030 

0.759 

0. 0040 

0.  895 

0.  0050 

1. 007 

0.  0070 

1.  i 87 

0.  0100 

1.  386 

0. 0150 

1,622 

0. 0200 

1.  794 

0.0  300 

2.  046 

0. 0400 

2.  234 

0. 0500 

2.  385 

0. 0700 

2.  628 

0. 1000 

2.  914 

0 1500 

3.  295 

0.  2000 

3.617 

0. 3000 

4.  151 

0.  4 000 

4.  557 

0.  5000 

4.  819 

mk 


TABLE  11-23 

ADAPTIVE  FILTERING  SIGNAL  DEGRADATION 
VERSUS  CONVERGENCE  RATE 
FOR  A WEAK  SIGNAL  FROM  300°- 305° 
(ALPA,  16  SITES,  19'.  1 DAY  276  21.  37.  01-21.  45.  32) 


Convergence  Rate 
(M) 

Signal  Degradation 
(dB) 

0.0010 

0.  203 

0.0015 

0.  288 

0.  0020 

0.  367 

0.  0030 

0.  507 

0.0040 

0.6  32 

0.0050 

0.  744 

0.0070 

0.  042 

0.  0100 

1.  191 

0.0150 

1.  517 

0. 0200 

1.  766 

0.  0300 

2.  137 

0.  0400 

2.423 

0. 0500 

2.663 

0.  0700 

3.  057 

0.  1000 

3.485 

0.  1500 

3.  895 

0.  2000 

4.  077 

0.  3000 

4.  104 

0.4000 

3.  921 

0.  5000 

3.670 
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FIGURE  11-46 


SIGNAL  DEGRADATION  AS  A FUNCTION  OF  CONVERGENCE 
RATE  FOR  A WEAK  SIGNAL  FROM  300°  -305°  (ALPA,  16  SITES 
1971  DAY  276  21.37.01-21.45.32) 


TABLE  11-24 


ADAPTIVE  FILTERING  SIGNAL-TO  - NOISE  GAIN  VERSUS 
CONVERGENCE  RATE  FOR  A WEAK  SIGNAL  FROM  300°- 305 
(ALPA,  FULL  ARRAY,  1971  DAY  276  21.  37.  01-21.45.32, 
USING  NOISE  FROM  1972  DAY  335  0415-0815) 


Convergence  Rate 
(M> 

Signal -to-Noise  Gain 
(dB) 

0.  0010 

' 0.154 

0.  0015 

0.  194 

0.  0020 

0.  220 

0. 0030 

0.  252 

0.  0040 

0.  263 

0. 0050 

0.257 

0.  0070 

0.  245 

0„  0100 

0.  195 

0.  0 l 50 

0.  105 

0=  0200 

0.  028 

0.  0300 

-0.  091 

0. 0400 

-0.  189 

0. 0500 

-0.  278 

0.  0700 

- 0.429 

0. 1000 

- 0.  571 

0. 1500 

- 0.  600 

0.  2000 

- 0.460 

00  3000 

0.  047 

0.  4000 

0.636 

0.  5000 

1 . 149 

1-t — 

0.  001 


..  .1.  ..... 


• M M«T 
,1*  •••■ 


0.01 
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FIGURE  Ii-47 

SIGNAL -TO -NOISE  GAIN  AS  A FUNCTION  OF  CONVERGENCE 
RATE  FOR  A WEAK  SIGNAL  FROM  300°-305°  (ALPA,  FULL  ARRAY, 
1971  DAY  276  21.  37.01  -21.  45.  32,  USING  NOISE  FROM  1972  DAY  335 

0415-0815) 
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Table  11-25  presents  the  power  reduction  relative  to  beamsteering  at 
H = 0.  005  and  H = 0.  5 over  the  specified  time  gates  for  five  events  with  azi- 
muths away  from  the  302.  5°  steer  direction.  As  in  the  case  of  the  corres- 
ponding six-channel  results  for  the  302.  5°  steer  direction  (Table  11-12),  all 
off-azimuth  events  are  reduced  by  more  than  12  dB  at  n = 0.5.  Figures  II- 
49  through  11-53  show  the  beamsteer  output  and  ABF  output  at  = 0.  005  and 
/i  = 0.5  for  the  five  events  in  Table  11-25.  The  amplitude  of  the  beamsteer 
output  in  the  top  trace  of  each  figure  is  smaller  than  the  corresponding  six- 
channel  beamsteer-output  amplitudes  since  the  time-shift-and-sum  amplifi- 
cation factors  for  the  full  array  (Figure  11-54)  are  lower  than  the  correspond- 
ing factors  for  the  six-channel  array  (Figure  11-18)  at  all  azimuths  away  from 
the  302.5°  steer  azimuth.  Figures  11-55  and  11-56  present  adaptive  beam  pat- 
terns for  H = 0.  005  at  the  beginning  and  middle  of  the  latter  Panama  event. 

The  8 dB  suppression  of  this  event  is  reflected  by  the  deepening  null  at  the 
113°  azimuth.  This  null  persists  until  the  end  of  the  event  at  2128,  when  the 
beam  pattern  is  essentially  the  same  as  at  2124.  In  contrast,  the  adaptive 
beamformer  achieved  the  6 dB  suppression  of  the  second  Kamchatka  event  at 
H = 0.  005  by  deflecting  the  main  lobe  of  the  beam  pattern  away  from  the  273° 
arrival  azimuth  of  the  event.  Figures  11-57  through  11-59,  which  show  the 
adaptive  beam  pattern  at  the  beginning,  middle,  and  end  of  the  event,  demon- 
strate this  deflection  process. 

The  last  two  signals  used  for  full-array  signal  degradation  measure- 
ments are  the  two  Kamchatka  events  from  an  azimuth  of  273°.  Again,  all 
processing  procedures  for  these  two  events  are  the  same  as  those  which  pro- 
duced the  signal  degradation  measurements  in  Subsection  C for  the  new  adap- 
tive algorithm  except  for  the  sites  employed  (all  but  7,  11,  and  18).  The  first 
of  these  events  arrived  at  about  2108.  Table  11-26  and  Figure  11-60  present 
the  full-array  signal  degradation  for  this  event  as  a function  of  convergence 
rate.  The  corresponding  signal-to-noise  gain  values  based  on  the  full-array 


FULL- ARRAY  BEAMS  FOR  MAGNITUDE  4.  5 KAMCHATKA  EVENT 
(DAY  276  1971,  STEER  DIRECTION  302.  5°) 


T I ME -SH I FT -AND -SUM  BERM  PRTTERN 
BERM  LOOK  VELOCITY  IS  3.6,  LOOK  RZIMUTH  302*5 
FREQUENCY  IS  0*0^000  HZ,  PERIOD  25*0  SECONDS 


FIGURE  11-54 

BEAMSTEER  RESPONSE  FOR  FULL-ARRAY  SAMPLE 
FROM  DAY  276  OF  1971 
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ABF  BEAM  PATTERN  AT  276/21 .20.00 
CONVERGENCE  RATE  IS  0*0050 
BEAM  LOOK  VELOCITY  IS  3*6,  LOOK  AZIMUTH  302*5 
FREQUENCY  IS  0*04000  HZ,  PERIOD  25*0  SECONDS 


FIGURE  11-55 


FULL- ARRAY  ABF  RESPONSE  FOR  MAGNITUDE 
4.  7 PANAMA  EVENT  AT  2120  ON  DAY  276  OF  1071 


ABF  BERM  PATTERN  RT  276/21.24.00 
CONVERGENCE  RATE  15  0.0050 
BERM  LOOK  VELOCITY  15  3.6.  LOOK  AZIMUTH  302.5 
FREQUENCY  IS  0.04000  HZ.  PERIOD  25.0  SECONDS 


FIGURE  II -‘SI) 

FUEL-ARRAY  ABF  RESPONSE  FOR  MAGNITUDE 
4.7  PANAMA  EVENT  AT  2124  ON  DAY  276  OF  1971 
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ALPA 

ABF  BEAM  PATTERN  AT  276/22.07.00 
CONVERGENCE  RATE  IS  0.0050 
BEAM  LOOK  VELOCITY  IS  3.6,  LOOK  AZIMUTH  302.5 
FREQUENCY  IS  0.04000  HZ.  PERIOD  25.0  SECONDS 


FIGURE  11-57 


FULL- ARRAY  ABF  RESPONSE  FOR  MAGNITUDE 
4.  9 KAMCHATKA  EVENT  AT  2207  ON  DAY  276  OF  1971 
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RBF  BERM  PRTTERN  RT  276/22* 10*  30 
CONVERGENCE  RRTE  IS  0.0050 
BERM  LOOK  VELOCITY  IS  3.6.  LOOK  RZIMUTH  302.5 
FREQUENCY  IS  0.04000  HZ.  PERIOD  25.0  SECONDS 


FIGURE  11-58 


FULL- ARRAY  ABF  RESPONSE  FOR  MAGNITUDE 
4.  9 KAMCHATKA  EsfENT  AT  2210  1/2  ON  DAY  276  OF  1971 


ABF  BERM  PATTERN  AT  276/22. 14.00 
CONVERGENCE  RATE  IS  0.0050 
BEAM  LOOK  VELOCITY  IS  3.6.  LOOK  AZIMUTH  302.5 
FREQUENCY  IS  0.04000  HZ,  PERIOD  25.0  SECONDS 


FIGURE  11-59 

FULL- ARRAY  ABF  RESPONSE  FOR  MAGNITUDE 
4.  9 KAMCHATKA  EVENT  AT  2214  ON  DAY  276  OF  1971 


TABLE  11-26 


ADAPTIVE  FILTERING  SIGNAL  DEGRADATION 
VERSUS  CONVERGENCE  RATE  FOR  A STRONG 
SIGNAL  FROM  KAMCHATKA  (ALPA,  16  SITES, 
1971  DAY  276  21. 07.  18-21.  1 1.  33) 


Convergence  Rate 

(M) 

0. 0010 
0.  0015 
0.  0020 
0.  0030 
0.  0040 
0.  0050 
0.  0070 
0.  0 100 
0.  0 150 
0.  0200 
0.  0300 
0.  0400 
0.  0500 
0.  0700 
0.  1000 
0.  1 500 
0„ 2000 
0. 3000 
0.  4000 
0.  5000 


Signal  Degradation 
(dB) 

0.  0 1 8 
0U  028 
0.  040 
0.  064 
0.  09! 

0.’I8 
0.  1 74 
0.  255 
0.  372 
0.  469 
0.613 
0.  71  1 
0.  778 
0.  851 
0.  874 
0.  812 
0.  706 
0.  486 
0.  305 
0.  173 


noise  reduction  measurements  for  day  335  of  1972  appear  in  Table  11-27  and 
Figure  11-61.  The  beamsteer  tr  ices  and  adaptive  filter  traces  at  = 0.  005 
and  4 = 0.  5 are  shown  for  this  event  in  Figure  11-62.  Some  signal  distortion 
is  evident  in  the  bottom  trace  for  the  highest  convergence  rate.  The  second 
Kamchatka  event  reached  AL,P  A at  approximately  2207.  The  full-array  sig- 
nal degradation  for  this  event  is  given  in  Table  11-28  and  Figure  11-63  as  a 
function  of  convergence  rate.  Using  the  full-array  noise  reduction  measure- 
ments from  day  335  of  1972,  the  resulting  signal-to-noise  ratio  improvement 
values  are  depicted  in  Table  11-29  and  Figure  11-64.  The  estimated  full  array 
signal-to-noise  gains  for  this  event  are  lower  than  the  corresponding  six- 
channel  values  at  all  convergence  rates.  The  greatest  difference  of  about 
2 dB  at  /i  = 0.  5 is  due  to  the  greater  signal  degradation  with  the  full  array. 
Figure  11-65  displays  the  beam  outputs  for  the  second  Kamchatka  event. 

Slight  attenuation  occurs  in  both  adaptive  filter  beams. 

Table  11-20  gives  the  full  array  power  reduction  relative  to  beamsteer- 
ing at  n = 0.  005  and  M = 0.  5 for  four  events  away  from  the  273°  steer  direc- 
tion in  the  data  sample  used  to  obtain  signal  degradation  measurements  for 
the  two  Kamchatka  events.  Event  suppression  at  4 = 0.  5 is  greater  than  12 
dB  for  all  events  but  the  first  Panama  event.  Figure  11-66  is  a Calcomp  plot 
of  the  beamsteer  output  and  ABF  outputs  at  4 = 0.005  and  4 = 0.5  for  the 
weak  event  from  300  -305  . The  corresponding  plot  for  the  Andreanof  Islands 
event  is  displayed  in  Figure  11-67  . 

The  signal-to-noise  ratio  improvement  of  adaptive  beamforming  rela- 
tive to  beamsteering  is,  in  geneial,  lower  with  the  full  array  than  with  the 
six  channels.  Differences  between  the  full  array  and  the  closely-spaced  six- 
channel  array  in  measured  signal  degradation  are  highly  variable.  This  fact 
is  particularly  evident  at  the  more  rapid  convergence  rates,  where  2 dB  dif- 
ferences occur  for  the  weak  event  from  300°-305°  and  the  second  Kamchatka 
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TABLE  11-27 

ADAPTIVE  FILTERING  SIGNAL-TO-NOISE  GAIN 
VERSUS  CONVERGENCE  RATE  FOR  A STRONG  SIGNAL 
FROM  KAMCHATKA  (A LPA , FULL  ARRAY,  1971  DAV  276 
21.  07.  18-21.  11.  33,  USING  NOISE  FROM  1972  DAY  335  0415-0815) 


Convergence  Rate 

Signal-to-Noise  Gain 

(M) 

(dB) 

0.0010 

0.  339 

0.0015 

0.454 

0.0020 

0.  547 

0.0030 

0.695 

0.0040  ? 

0.  804 

0.  0050 

0.  889 

0.  0070 

1.  013 

0.  0100 

1.  131 

0.0150 

1.  250 

0.  0200 

1.  325 

0.  0330 

1. 433 

0.  0400 

1.523 

0.  0500 

1.607 

0.  0700 

1.  777 

0. 1000 

2.040 

0.  1500 

2.483 

0.  2000 

2.  911 

0. 3000 

3.665 

0.4000 

4.  252 

0. 5000 

4.646 

H 
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FIGURE  11-61 

SIGNAL  - TO -NOISE  GAIN  AS  A FUNCTION  OF  CONVERGENCE 
RATE  FOR  A STRONG  SIGNAL  FROM  KAMCHATKA  (ALPA,  FULL 
ARRAY,  1971  DAY  276  21 . 07.  1 8-21.  1 1 . 33,  USING  NOISE  FROM  1972 

DAY  335  041  5-081  5) 
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MAGNITUDE  4.5  KAMCHATKA  EVENT 
STEER  DIRECTION  27 


TABLE  11-28 

ADAPTIVE  FILTERING  SIGNAL  DEGRADATION 
VERSUS  CONVERGENCE  RATE  FOR  A VERY 
STRONG  SIGNAL  FROM  KAMCHATKA  (ALPA, 
16  SITES,  1971  DAY  276  22.  07.  02-22.  1 1.  17) 


I 

I 


C 


onvergence  Rate 


Signal  Degradation 


(M) 


(dB) 


0. 0010 
0.0015 
0. 0020 
0.  0030 
0. 0040 
0. 0050 
0.  0070 
0.  0100 
0. 0150 
0. 0200 
0.  0300 
0. 0400 
0.0500 
0.  0700 
0. 1000 
0. 1500 
0. 2000 


0.088 
0.  131 
0.  173 
0.253 
0.  325 
0.  392 
0.  506 
0.642 
0.  804 

0.  917 

1.  056 
1.  129 
1.  167 
1.  192 
1.  185 
1.  162 
1.  150 
1.  125 
1.068 
0.  973 


0. 3000 
0. 40uu 
0. 5000 


i 


*• 


Convergence  Rate  M 


FIGURE  11-63 


SIGNAL  DEGRADATION  AS  A FUNCTION  OF  CONVERGENCE 
RATE  FOR  A VERY  STRONG  SIGNAL  FROM  KAMCHATKA 
(ALPA,  16  SITES,  1971  DAY  276  22.  07.  02-22.  11.  17) 
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TABLE  11-29 

ADAPTIVE  FILTERING  SIGNAL-TO-NOISE 
GAIN  VERSUS  CONVERGENCE  RATE  FORA 
VERY  STRONG  SIGNAL  FROM  KAMCHATKA 
(ALPA,  FULL  ARRAY,  1 97  1 DA Y 276  22.  07.  02 -22.  1 1 . 1 
USING  NOISE  FROM  1972  DAY  335  04 '5-0815) 


Convergence  Rate 
(U) 

Signal-to-Noise  Gain 
(dB) 

0. 0010 

0.  269 

0. 0015 

0.  351 

0o  0020 

0.414 

0.  0030 

0.  506 

0. 0040 

0.  570 

0. 0050 

0.6  15 

0.  0070 

0.  681 

0 0 100 

0.  744 

0.  0 150 

0.  818 

0. 0200 

0.  877 

0.  0300 

0.  990 

0„  0400 

1.105 

0.  0500 

1.218 

0.  0700 

1. 4 36 

0.  1000 

l.  729 

0.  1500 

2.  i 33 

0. 2000 

2.467 

0.  3000 

3.  026 

0.  4000 

3.  489 

0o 5000 

3a  846 
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FIGURE  11-64 


SIGNAL  - TO -NOISE  GAIN  AS  A FUNCTION  OF  CONVERGENCE 
RATE  FOR  A VERY  STRONG  SIGNAL  FROM  KAMCHATKA 
(ALFA,  FULL  ARRAY,  1971  DAY  276  22.07.02-22.11.17, 
USING  NOISE  FROM  1972  DAY  335  0915-0815) 
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FULL -ARRAY  BEAMS  FOR  A WEAK  EVENT  FROM  300 
(DAY  276  1971,  STEER  DIRECTION  273°) 


event.  Differences  in  measured  signal  degradation  are  not  the  main  reason 
for  the  reduced  improvement  relative  to  beamsteering  with  the  full  array. 

On  the  contrary,  the  principal  factor  contributing  to  this  result  is  the  mea- 
sured noise  reduction  for  the  data  sample  from  0415  to  0815  on  day  335  of 
1972:  at  all  convergence  rates,  the  noise  reduction  is  lower  with  the  full 
array  than  with  six  channels.  One  of  the  reasons  for  the  lower  full-array 
noise  reduction  is  that  the  more  widely-spaced  sensors  of  the  full  array  have 
less  coherence  between  sensor  pairs  than  the  closely-packed  elements  of  the 
array  used  for  the  six-channel  measurements.  Another  reason  is  that,  for 
the  arrays  employed  in  the  noise  reduction  measurements,  the  ratio  of  the 
directivity  index 

/ 2 n 

|A(f,0)|2  d 6 

I 0 

for  the  fifteen-channel  array  to  that  for  the  six-channel  array  is  greater  than 
15/6  (*4  dB).  In  the  formula  above,  A(f,  0)  is  the  array  amplification  fac- 
tor for  a plane  wave  of  frequency  f arriving  from  an  azimuth  of  (36O/2tt)0 
degrees.  Figures  11-68  and  11-69,  which  picture  the  time-shift-and-sum  beam 
patterns  for  the  six-channel  array  and  fifteen-channel  array,  respectively,  at 
0.04  Hr  and  3.6  km/sec,  demonstrate  that  the  directivity- index  ratio  is  great- 
er than  4 dB.  This  fact  i’  '.ies  that  the  ratio  of  coherent  energy  to  spatially 
uncorrelated  energy  in  the  beamsteer  output  is  less  for  the  fifteen-channel 
array  than  for  the  six-channel  array  if  coherent  energy  is  equally  likely  at 
all  azimuths.  As  a result,  the  potential  signal-to-noise  gain  of  adaptive  beam- 
forming relative  to  beamsteering  is  less  for  fifteen  channels  than  for  six  chan- 
nels. 
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T I ME-SH I FT-RND-SUM  BERM  PATTERN 
BERM  LOOK  VELOCITY  IS  3.6.  LOOK  AZIMUTH  270.0 
i-REQUENCY  IS  O.OHOOO  HZ.  PERIOD  25.0  SECONDS 


FIGURE  11-68 

BEAMSTEER  RESPONSE  FOR  SIX-CHANNEL 
NOISE  SAMPLE  FROM  DAY  335  OF  1972 
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T I ME-SH I FT-^ND— SUM  BERM  PATTERN 
BERM  LOOK  VELOCITY  IS  3.6,  LOOK  AZIMUTH  270,0 
(-REQUENCY  IS  Q.QHQQQ  HZ.  PERIOD  25,0  SECONDS 


FIGURE  11-69 

REAMSTEER  RESPONSE  FOR  FULL- ARRAY 
NOISE  SAMPLE  FROM  DAY  336  OF  J 972 
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E.  FREQUENCY  FILTERING  FOR  IMPROVED  DETECTION 

So  far,  the  focus  of  attention  in  this  report  has  been  the  impioved  de- 
tection performance  obtained  from  the  directional-filtering  capability  of  adap- 
tive beamforming.  Another  important  factor  affecting  detector  performance 
is  the  relative  weighting  given  to  the  various  frequency  components  in  the  data. 
In  short,  what  frequency  filter  should  be  applied  to  the  data  in  order  to  achieve 
best  detection  results?  An  additional  question  is  whether  this  frequency  filter 
should  be  applied  to  the  vertical-component  data  of  the  individual  sites  before 
adaptive  beamforming  or  whether  it  should  be  applied  to  the  beam  outputs  after 
adaptive  filtering.  In  a time-domain  adaptive  filtering  system  such  as  the  one 
employed  in  this  study,  the  answer  to  this  question  is  that  the  frequency  filter 
jhould  be  applied  before  adaptive  processing  so  that  the  adaptive  beamformer 
can  concentrate  on  achieving  processing  gain  at  the  frequencies  most  useful 
for  detection.  If,  instead,  the  frequency  filter  is  applied  to  the  beam  outputs, 
the  adaptive  processing  gains  at  frequencies  relevant  to  detection  are  reduced 
by  the  mean- square-power  performance  criterion,  which  forces  the  adaptive- 
filtering  system  to  react  to  conditions  at  frequencies  of  little  use  in  detection 
processing.  The  precise  nature  of  the  prefilter  to  be  applied  to  the  site  ver- 
tical component  data  before  input  to  the  multichannel-filtering  system  has  yet 
to  be  discussed.  The  remainder  of  this  subsection  is  devoted  to  the  develop- 
ment of  a technique  for  determining  the  frequency  weightings  for  this  prefil- 
ter. 

One  frequently  employed  detection  system  is  a square-law  (or  power) 
detector,  in  which  the  average  power  of  a beam  over  some  specified  time  gate 
is  the  detector  output.  Since  noise  fluctuates  in  a random  fashion,  the  noise 
output  power  over  the  specifit  d gate  length  has  a characteristic  probability 
density  function.  If  a specific  signal  (uncorrelated  with  the  noise)  is  added  to 
the  noise,  the  resulting  power  probability  density  function  shifts  upward  by  an 
amount  equal  to  the  signal  power  over  the  detection  gate.  Figure  11-70 
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illustrates  this  situation.  If  the  noise  power  has  a Gaussian  probability  dis- 
tribution, the  best  detection  results  occur  when  the  ratio  of  the  signal  power 
to  the  noise-power  standard  deviation  is  a maximum.  The  best  frequency 
filter  for  detection  is  the  filter  which  maximizes  this  ratio.  Suppose  that  a 
discrete  Fourier  transform  is  taken  over  the  detection  gate.  Let  $ (f  ) 

S S K 

and  <t>  (f  ) denote  the  squared  amplitude  of  the  signal  and  noise,  respec- 
nn  k 

tively,  at  the  discrete-Fourier-transform  frequency  f^  . Let  W(f^)  denote 

the  power  weight  given  to  the  frequency  f . Then  the  signal  power  over  the 

T k T T 

detection  gate  is  W S,  and  the  noise  power  variance  is  W ZZ  W , where 

the  K-component  vectors  W,  S,  and  Z are,  respectively, 


W(fj) 

<t>  (f,)- 

nn  1 

<t>  (f.) 

nn  1 

W(f2) 

f2> 

4>nn(y  - 

nn  L 

<t>  (f.) 

nn  2 

• 

s = 

• 

and  Z = 

• 

• 

• 

• 

• 

• 

• 

> 

9 

Here  <t>  (f,  ) denotes  the  expectation  of  the  squared  noise  amplitude  $ (f  ) 

nn  k ' nn  k 

at  the  frequency  f . That  is  to  say,  ^^f^  is  the  noise  power  density 

T 

spectrum  at  f . The  matrix  expectation  ZZ  is  the  covariance  matrix  for 

T / / T T 

the  noise  power  density.  The  ratio  W S / v W ZZ  W or,  equiv  ilently, 

T T / T T 

W SS  W / W ZZ  W is  maximized  when  the  weight  vector  W is  propor- 

T -1  ... 

tional  to  (ZZ  ) S.  Hopefully  all  components  of  W are  positive.  This  con- 
dition is  met  when  the  quantities  4>  (f,  ) - <t>  (f,  ) are  independent  from 

nn  k nn  k 

'T* 

frequency  to  frequency,  so  that  ZZ  is  a diagonal  matrix.  To  design  a filter 

in  this  v/ay,  fourth-order  moments  must  be  measured.  No  software  is  cur- 

T 

rently  available  for  these  measurements,  so  that  matrices  ZZ  are  not 
available. 
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If  the  noise  is  colored  and  Gaussian,  simplifying  assumptions  are  pos- 
sible. As  an  artifice,  let  the  data  be  passed  through  a noise-whitening  filter, 

so  that  the  elements  <t>  (f.  ) of  S are  replaced  by  * (f  ) / <J>  (t~)  and  the 

nn  k y ^Ss'  k' / nn'  k' 

elements  ~ ^ by  [$nn(f) J - 1 . Since  the  noise 

is  now  white  and  Gaussian,  the  variables  $ (f,  ) / <$  (fTT  are  independent 

2 nn  k / ^ nn  k v 

and  each  have  a X distribution,  so  that  their  variance  is  proportional  to 

their  expectation,  which  is  one.  Then  the  quantity  WTSS W*1  / WTZZTW  to 
be  maximized  is  proportional  to  W SS  W / W W , or 


W 

W 


Maximization  is  achieved  when  the  vector  W points  toward  the  same  direc- 
tion as  the  vector  S . Thus  the  components  of  the  weight  vector  W are  pro- 


portional to  <t>  (f.)  / <J>  (f  ) , and  the  combined  power  response  of  the  noise- 

s s k / n n k 

whitening  filter  and  the  frequency  weighting  W(f  ) is  proportional  to 

K 


• ss<V  /[^nn^k^  J * The  assoc*atecl  amplitude  weighting  is  proportional 
to  the  signal  amplitude  divided  by  the  noise  power.  As  the  equivalent  number 
of  degrees  of  freedom  in  the  filtered  noise  power  spectrum  increases,  the 
noise  power  probability  distribution  approaches  a Gaussian  distribution,  so 


that  the  amplitude  weighting  /<j>  (fj~  / <J>  (f)  is  lear-optimum  for  detec - 

y ss  f nn 

tion.  The  equivalent  number  of  degrees  of  freedom  in  the  filtered  noise  power 
spectrum  can  be  increased  by  lengthening  the  specified  detection  gate  length 

Using  the  power  spectrum  of  the  weak  event  from  300°-305°  in  the  data 

sample  from  day  276  of  1971  and  some  noise  spectra  from  late  autumn  of  1971, 

a detection  filter  was  specified  with  amplitude  response  /$  (7)~  / <t>  (f) 

V ss  / nn 

A zero-phase  time-domain  convolution  filter  was  designed  to  approximate  this 
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amplitude  response.  Figure  11-71  pictures  the  response  of  the  time-domain 
filter  as  a function  of  frequency.  Note  the  slight  depression  at  periods  be- 
tween 16  and  18  seconds.  In  the  6 to  8-second  band,  the  desired  zero  re- 
sponse was  weighted  especially  heavily  to  produce  the  lower  sidelobes  there. 
The  prefilter  whose  response  is  shown  in  Figure  11-71  is  used  in  all  long- 
period  processing  for  the  remainder  of  this  report  unless  specifically  noted 
othe  rwise. 

I ADAPTIVE  PROCESSING  WITH  THE  NEW  PREFILTER 

The  purpose  of  this  subsection  is  to  determine  the  effect  of  the  pre- 
filter in  Figure  11-71  on  adaptive  beamforming  signal-to-noise  ratio  improve- 
ment relative  to  beamsteering.  Since  this  prefilter  is  applied  to  data  sam- 
pled once  every  two  seconds,  an  adaptive-filter  length  of  fifteen  points  per 
channel  is  sufficieit  to  span  all  possible  propagation  delays  across  ALPA  at 
the  3.  6 km/sec  apparent  velocity  of  Rayleigh  waves.  In  the  beginning  of  this 
subsection,  tne  adaptive-filter  length  remains  at  the  3 1 -point-per-channel 
value  used  previously  in  this  section.  The  final  part  of  this  subsection  ex- 
amines the  effect  of  cutting  the  adaptive-filter  length  to  fifteen  points  per 
channel. 

Noise  reduction  measurements  with  the  3 1 -point-long  adaptive  filter 
employ  the  data  sample  from  0415  to  0815  on  day  335  of  1972.  Processing  is 
exactly  the  same  as  in  Subsection  D except  for  a two-second  sample  interval 
and  the  application  of  the  new  prefilter  in  Figure  11-71.  Table  11-31  and  Fig- 
ure 11-72  present  the  noise  reduction  as  a function  of  the  convergence  para- 
meter jU  . A comparison  of  Table  11-3'  with  Table  II -2 _ in  Subsection  D re- 
veals that  the  noise  reduction  is  significantly  greater  with  data  passed  through 
the  new  prefilter.  The  probable  explanation  for  the  increased  noise  reduction 
is  the  additional  emphasis  given  to  the  more  coherent  energy  at  low  freauenc- 
ie  s. 
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FIGURE  11-72 


NOISE  REDUCTION  AS  A FUNCTION  OF  CONVERGENCE  RATE  WITH 
31-POINT-LONG  ADAPTIVE  FILTER  (ALPA,  15  SITES,  1972  DAY  335 

0415-0815) 
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The  same  three  events  in  the  data  sample  from  2000  to  2345  on  day 
276  of  1971  provide  signal  degradation  measurements  for  the  31 -point-long 
adaptive  filter.  The  first  of  these  is  the  weak  event  from  300°- 305°.  Table 
11-32  and  Figure  11-73  portray  signal  degradation  with  the  new  prefilter  as  a 
function  of  convergence  rate  over  a 512-second  gate.  For  this  event,  signal- 
to-noise  gain  estimates  based  on  the  noise  sample  from  day  335  of  1972  ap- 
pear in  Table  11-33  and  Figure  11-74.  The  improvement  relative  to  beam- 
steering is  slightly  higher  than  with  the  old  prefilter  (Table  11-24)  at  all  con- 
vergence rates. 

At  the  302.  5°  steer  direction  used  for  the  weak  event  from  30f'°-305°, 
there  are  five  off-azimuth  events  in  the  data  sample  from  day  276  of  1971. 
Table  11-34  gives  the  power  reduction  relative  to  beamsteering  at  p = 0.  007 
and  H - 0.  5 for  these  events.  For  the  two  Panama  events  and  the  Andreanof 
Islands  event,  the  suppression  of  these  events  is  comparable  to  the  suppres- 
sion obtained  with  the  old  prefilter  (Table  11-25).  For  the  two  Kamchatka 
events  30°  away  from  the  steer  direction,  however,  the  attenuation  is  signif- 
icantly lower  with  the  new  prefilter.  The  emphasis  given  to  low-frequency 
energy  by  the  new  prefilter  not  only  increases  the  effective  beamwidth  of  the 
adaptive  beamformer  but  also  concentrates  more  energy  near  the  beginning 
of  the  events,  so  that  the  adaptive  filter  set  cannot  react  as  effectively  to  the 
presence  of  these  events. 

These  two  Kamchatka  events  arriving  from  a 2 7 3°  azimuth  also  furnish 
the  remaining  signal  degradation  measurements  for  the  31 -point-long  adaptive 
filter.  The  first  arrived  at  about  2108.  For  this  event..  Table  11-35  and  Fig- 
ure 11-75  depict  signal  degradation  with  the  new  prefiller  as  a function  of  the 
convergence  parameter  g . The  resulting  signal-to-noise  gain  estimates 
corresponding  to  the  noise  sample  from  day  335  of  1072  appear  in  Table  11-36 
and  Figure  11-76.  Maximum  signal-to-noise  gain  of  6.  36  dl3  relative  to  beam- 
steering at  g = 0.  5 is  significantly  higher  than  the  corresponding  4.65  fill 
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FIGURE  11-73 


SIGNAL  DEGRADATION  AS  A FUNCTION  OF  CONVERGENCE  RATE  WITH 
31 -POINT-LONG  ADAPTIVE  FILTER  FOR  A WEAK  SIGNAL  FROM 
300°-305°  (ALPA,  16  SITES,  1 971  DAY  276  21 . 37.  01  - 21 . 45.  32) 


TABLE  11-33 

SIGNAL -TO -NOISE  GAIN  VERSUS  CONVERGENCE  RATE  WITH 
31 -POINT-LONG  ADAPTIVE  FILTER  FOR  A WEAK  SIGNAL  FROM 
300°-305°  (ALPA,  FULL  ARRAY,  1971  DAY  276  21 . 37.  01  -21 . 45.  32, 
USING  NOISE  FROM  1972  DAY  335  0415-0815) 
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SIGNAL -TO-NOISE  GAIN  AS  A FUNCTION  OF  CONVERGENCE 
RATE  WITH  31 -POINT-LONG  ADAPTIVE  FILTER  FOR  A WEAK 
SIGNAL  FROM  300°-305°  (ALPA,  FULL  ARRAY,  1971  DAY  276 
21.  32.  01-21. 45.  32,  USING  NOISE  FROM  1972  DAY  335  0415-0815) 
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TABLE  11-35 


SIGNAL  DEGRADATION  VERSUS  CONVERGENCE  RATE  WITH 
31 -POINT-LONG  ADAPTIVE  FILTER  FOR  A STRONG  SIGNAL 
FROM  KAMCHATKA  (ALPA,  16  SITES,  1971  DAY  276  21 . 07.  1 8-21 . 1 1 . 33) 


Convergence  Rate 

Signal  Degradation 

(M) 

(dB) 
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TABLE  11-36 


SIGNAL  -TO -NOISE  GAIN  VERSUS  CONVERGENCE  RATE  WITH 
31 -POINT-LONG  ADAPTIVE  FILTER  FOR  A STRONG  SIGNAL  FROM 
KAMCHATKA  (ALPA,  FULL  ARRAY,  1971  DAY  276  21.07.18-21.11.33 
USING  NOISE  FROM  1972  DAY  335  0415-0815) 
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gain  with  the  old  prefilter.  The  increased  adaptive  filtering  gain  due  to  the 
new  prefilter  is  an  added  benefit  over  and  above  any  detection  improvement 
obtained  through  frequency  filtering.  The  second  of  the  two  Kamchatka  events 
reached  ALPA  at  approximately  22  07.  Table  11-37  and  Figure  11-77  present 
signal  degradation  as  a function  of  convergence  rate  with  the  new  prefilter  for 
this  event.  The  signal-to-noise  gain  estimates  obtained  by  combining  the  sig- 
nal degradation  measurements  for  the  second  Kamchatka  event  with  noise  re- 
duction measurements  for  day  335  of  1972  are  given  in  Table  11-38  and  Figure 
11-78.  In  this  case,  the  signal-to-noise  ratio  improvement  relative  to  beam- 
steering at  H = 0.  5 is  more  than  four  dB  greater  than  with  the  old  prefilter 
(Table  11-29). 

Table  11-39  gives  the  power  reduction  relative  to  beamsteering  at 
= 0.  00  7 and  ^ = 0.  5 for  four  events  away  from  the  273°  azimuth  of  the  two 
Kamchatka  events  in  the  data  sample  from  day  276  of  1971.  Event  suppres- 
sion for  the  Andreanof  Islands  event  36°  from  the  look  direction  is  less  for 
the  new  prefilter  than  for  the  old  prefilter.  Otherwise  the  event-suppression 
results  are  comparable  for  both  prefilters.  The  fact  that  the  ABF  output  at 
- 0.007  is  slightly  stronger  than  the  beamsteer  output  for  the  first  Panama 
event  is  due  to  an  in  tial  ABF  response  greater  than  the  beamsteer  response 
toward  the  113  azimu'h  of  that  event:  the  new  prefilter's  enhancement  of  the 
low-frequency  energy  near  the  beginning  of  the  ever,  does  not  allow  sufficient 
time  to  attenuate  the  event  on  the  adaptive  beam  output. 

Past  experience  with  multichannel  filtering  systems  has  shown  that  a 
filler  length  greater  than  the  travel  time  across  an  array  seldom  produces 
any  significant  additional  processing  gain.  In  the  case  of  adaptive  beamform- 
ing, moreover,  signal  degradation  might  conceivably  increase  substantially 
with  a greater  filter  length  than  necessary.  Likewise,  another  reason  for  ex- 
amining the  affect  of  adaptive  filter  length  is  the  resulting  computational  load, 
which  is  proportional  to  the  number  of  components  in  the  adaptive  filter  vector. 
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TABLE  11-37 


SIGNAL  DEGRADATION  VERSUS  CONVERGENCE  RATE  WITH 
31 -POINT-LONG  ADAPTIVE  FILTER  FOR  A VERY  STRONG  SIGNAL  FROM 
KAMCHATKA  (ALPA,  16  SITES,  1971  DAY  276  22.  07.  02-22.  1 1 . 17) 


Convergence  Rate 
(/0 

Signal  Degradation 
(dB) 

0. 0010 

0.  059 

0.  0015 

0.  087 

0. 0020 

0.113 

0.  00  30 

0.  163 

0.  0040 

0.  209 

0.  00  50 

0.  252 

0. 0070 

0.  328 

0. 0100 

0.  423 

0.  0150 

0.  543 

0. 0200 

0.630 

0. 0300 

0.  747 

0. 0400 

0.  814 

0.  0500 

0.  851 

0. 0700 

0.  879 

0. 1000 

0.  868 

0. 1500 

0.  800 

0.  2000 

0.  707 

0.  3000 

0.  506 

0.  4000 

0.  314 

0.  5000 

0.  137 
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FIGURE  11-77 

SIGNAL  DEGRADATION  AS  A FUNCTION  OF  CONVERGENCE  RATE  WITH 
31 -POINT-LONG  ADAPTIVE  FILTER  FOR  A VERY  STRONG  SIGNAL  FROM 
KAMCHATKA  (ALPA,  16  SITES,  1971  DAY  276  22.07.02-22.11.17) 
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TABLE  II- 38 


SIGNAL -TO -NOISE  GAIN  VERSUS  CONVERGENCE  RATE  WITH 
31 -POINT-LONG  ADAPTIVE  FILTER  FOR  A VERY  STRONG  SIGNAL  FROM 
KAMCHATKA  (AL PA,  FULL  ARRAY,  1971  DAY  276  22.07.02-22.11.  17, 
USING  NOISE  FROM  1972  DAY  335  0415-0815) 


Convergence  Rate 


Signal-to-Noise  Gain 


(M) 

0. 0010 
0. 0015 
0.  0020 
0. 0030 
0.  0040 
0. 0050 
0. 0070 
0. 0100 
0. 0150 
0. 0200 
0.  0300 
0. 0400 
0. 0500 
0. 0700 
0.  1000 
0. 1500 
0. 2000 
0. 3000 
0.  4000 
0. 5000 


(dB) 
0.  440 
0.  592 

0.  717 
0.  912 
1. 059 

1.  175 
1.  350 
1.  535 
1. 754 

1.  925 

2.  212 

2.  469 

2.  712 

3.  164 

3.  772 

4.  645 

5.  387 

6.  568 

7.  435 

8.  067 


i 
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FIGURE  11-78 

SIGNAL -TO -NOISE  GAIN  AS  A FUNCTION  OF  CONVERGENCE  RATE  WITH 
31 -POINT-LONG  ADAPTIVE  FILTER  FOR  A VERY  STRONG  SIGNAL  FROM 
KAMCHATKA  (ALPA,  FULL  ARRAY,  1971  DAY  276  22.07.02-22.11.17, 
USING  NOISE  FROM  1972  DAY  335  0415-0815) 


If  there  are  no  meaningful  processing  losses  at  a reduced  filter  length,  then 
operational  considerations  favor  a shorter  filter.  For  these  reasons  just 
discussed,  the  remainder  of  this  subsection  presents  results  obtained  with 
an  adaptive  filter  length  of  fifteen  points  per  channel. 

The  data  sample  from  0415  to  0815  on  day  335  of  1972  provides  noise 
reduction  measurements  at  the  shorter  filter  length.  The  steer  azimuth  is 
again  270  . Table  11-40  and  Figure  11-79  give  the  noise  reduction  as  a func- 
tion of  convergence  rate.  At  H = 0.  04  and  below,  the  noise  is  actually  lower 
with  the  fifteen-point  filter  than  with  the  31-point  filter  (Table  11-31).  At 
/i=  0.5,  the  longer  filter  produces  only  about  0.  3 dB  more  noise  reduction. 

The  weak  event  from  300°-305°,  the  two  Kamchatka  events  from  273°, 
and  the  New  Britain  event  from  240°  in  the  data  sample  from  2 000  to  2345  on 
day  276  of  1971  furnish  signal  degradation  measurements  for  the  15-point- 
long  adaptive  filter.  Table  II-4  t and  Figure  11-80  depict  signal  degradation 
for  the  weak  event  from  300°-305°  at  the  shorter  filter  length.  For  this  event. 
Table  11-42  and  Figure  11-81  reveal  that  the  estimated  signal-to-noise  gain  us- 
ing the  noise  sample  from  day  335  of  1972  is  higher  at  all  convergence  rates 
with  a 15-point  filter  than  with  a 31 -point  filter  (Table  11-33).  In  fact,  the 
difference  is  more  than  1.8  dB  at  (i=  0.5.  For  this  event,  where  signal  sim- 
ilarity across  the  array  is  relatively  poor,  the  shortened  adaptive -filter  length 
seems  to  obstruct  the  signal-degradation  process . Figure  11-82  pictures  the 
beamsteer  output  and  adaptive  beams  at  H = 0.  007  and  = 0.  5 for  the  weak 
event  from  300  -305  . This  figure  provides  the  first  opportunity  to  examine 
the  time  traces  for  beams  formed  from  vertical-component  data  passed  through 
the  new  prefilter.  Compared  with  the  corresponding  traces  for  the  old  prefil- 
ter in  Figure  11-48,  the  beam  outputs  have  a much  smoother  appearance  due 
to  the  elimination  of  high-frequency  energy  by  the  new  prefilter.  The  lower 
signal  degradation  with  the  new  prefilter  at  the  maximum  convergence  rate  is 
reflected  in  the  bottom  trace,  where  more  cycles  of  the  event  are  preserved 
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TABLE  11-40 


NOISE  REDUCTION  VERSUS  CONVERGENCE  RATE  WITH 
15-POINT-LONG  ADAPTIVE  FILTER  (ALPA,  15  SITES,  1972  DAY  335 

0415-0815) 


I 


Convergence  Rate 

Noise  Reduction 

(M> 

(dB) 

0.0010 

0.  595 

. 

0.0015 

0.  789 

0.  0020 

0.  946 

0.  0030 

1.  193 

0.0040 

1.  383 

! , | 

- / 

0.00  50 

1.  539 

0.  0070 

1.  787 

0.  0100 

2.064 

0.  0150 

2.  394 

0.0200 

2.638 

j 

0.0300 

3.007 

0.0400 

3.  297 

- 1 

0.  0500 

3.  547 

- , 

0.  0700 

3.974 

- 

0.  1000 

4.  511 

0.  1500 

5.  248 

0.  2000 

5.  853 

0.  3000 

6.  780 

0.  4000 

7.  430 

0.  5000 

7.  876 

i I 


FIGURE  11-79 

REDUCTION  AS  A FUNCTION  OF  CO> 
NT-LONG  ADAPTIVE  FILTER  (ALPA 

0415-0815) 
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TABLE  11-41 

SIGNAL  DEGRADATION  VERSUS  CONVERGENCE  RATE  WITH 
15-POINT-LONG  ADAPTIVE  FILTER  FOR  A WEAK  SIGNAL  FROM 
300°-305°  (ALPA,  16  SITES,  1971  DAY  276  21.37.01  -21.45.32) 


Convergence  Rate 
(10 

Signal  Degradation 
(dB) 

0.  0010 

0.  078 

0.  0015 

0.  127 

0.  0020 

0.  180 

0.  0030 

0.293 

0.  004  0 

0.411 

0.  0050 

0.530 

0.  0070 

0.761 

0.  0100 

1. 081 

0.  0150 

1.  535 

0.  0200 

1.  902 

0.  0300 

2.448 

0.  0400 

2.844 

0.  0500 

3.  161 

0.  0700 

3.  667 

0.  1000 

4.248 

0.  1500 

4.  922 

0.  2000 

5.  325 

0.  3000 

5.519 

0.4000 

5.  166 

0.  5000 

4.  538 
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Convergence  Rate  M 


FIGURE  11-81 

SIGNAL -TO -NOISE  GAIN  AS  A FUNCTION  OF  CONVERGENCE  RATE  WITH 
15-POINT-LONG  ADAPTIVE  FILTER  FOR  A WEAK  SIGNAL  FROM 
300° -30 5°  (ALPA,  FULL  ARRAY,  1971  DAY  276  21.37.01-21.4  5.32, 
USING  NOISE  FROM  1972  DAY  335  0415-0815) 
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than  in  the  corresponding  trace  of  Figure  11-48.  The  ripples  at  2135  in  the 
upper  two  traces  correspond  to  a data  glitch  at  site  2 and  are  eliminated  in 
the  adaptive  beam  output  for  b = 0.  5. 

Table  11-43  presents  the  power  reduction  relative  to  beamsteering  at 
fl-  0.  007  and  b = 0.5  with  the  1 5 -point-long  adaptive  filter  for  the  five  events 
away  from  the  302.  5 steer  direction  used  for  the  weak  event  from  300  -305  . 

The  off-azimuth  event  suppression  is  close  to  that  obtained  with  the  31 -point- 
long  adaptive  filter  (Table  11-34).  At  /*  = 0.  007,  all  but  the  first  Panama 
event  are  attenuated  more  with  the  15 -point-long  adaptive  filter.  At  b - 0.5, 
three  events  are  reduced  more  with  the  longer  adaptive  filter,  two  events  less. 
Thus  off-azimuth  event  suppression  is  hurt  little,  if  any,  by  the  shorter  filter 
length.  Figures  11-83  through  11-87  are  Calcomp  plots  of  the  302.  5 beams 
obtained  from  data  passed  through  the  new  prefilter  for  the  five  events  in 
Table  11-43.  Compared  with  the  corresponding  beams  for  the  old  prefilter 
in  Figure  11-49  through  11-53  of  Subsection  D,  the  trace  outputs  contain  less  high- 
frequency  noise  and  tend  to  concentrate  energy  near  the  beginning  of  the  events. 
In  the  case  of  the  two  Kamchatka  events  30°  from  the  look  azimuth,  a cycle  or 
so  at  the  beginning  of  the  events  escapes  suppression  on  the  adaptive  beam 
outputs  for  b - 0.5. 

The  two  Kamchatka  events  in  the  data  sample  from  day  276  of  1971  sup- 
ply signal  degradation  measurements  for  the  15-point-long  adaptive  filter  at 
a steer  azimuth  of  273  • Signal  degradation  for  the  first  of  these  events  (ar- 
riving at  about  2108)  appears  in  Table  11-44  and  Figure  11-88.  The  resulting 
signal-to-noise  gain  corresponding  to  the  noise  sample  from  day  335  is  given 
in  Table  11-45  and  Figure  11-89.  At  the  maximum  convergence  rate,  the  esti- 
mated signal-to-noise  ratio  improvement  is  0.3  to  0.4  dB  lower  with  the  15- 
point  filter  length  than  with  the  31 -point  filter  length  (Table  11-36).  Figure 
11-90  displays  the  beamsteer  output  and  adaptive -filter  outputs  at  M = 0.  007 
and  b - 0.  5 for  this  event.  Signal  distortion  and  some  signal  degradation  are 
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BEAMS  FOR  MAGNITUDE  4.7  PANAMA  EVENT  WITH  1 5-POINT  - LONG 
ADAPTIVE  FILTER  (DAY  276  1971,  STEER  DIRECTION  302.5°) 


BEAMS  FOR  ANDREANOF  ISLANDS  EVENT  WITH  1 5 - POINT  - LONG 
ADAPTIVE  FILTER  (DAY  276  1971,  STEER  DIRECTION  302.  5°) 


TABLE  11-44 


SIGNAL  DEGRADATION  VERSUS  CONVERGENCE  RATE  WITH 
15-POINT-LONG  ADAPTIVE  FILTER  FOR  A STRONG  SIGNAL  FROM 
KAMCHATKA  (ALP A,  16  SITES,  1971  DAY  276  21.07.18-21.11.33) 


Convergence  Rate 

ial 

Signal  Degradation 
(dB) 

0. 0010 

0.  084 

0.  0015 

0.  121 

0. 0020 

0.  154 

0.  0030 

0.  216 

0. 0040 

0.  271 

0.0050 

0.  321 

0. 0070 

0.  411 

0.0100 

0.  529 

0.0150 

0.  694 

0. 0200 

0.  828 

0. 0300 

1. 036 

0. 0400 

1.  191 

0. 0500 

1.  310 

0.0700 

1. 476 

0.  1000 

1.  619 

0.  1500 

I.  728 

0.  2000 

1.  784 

0.  3000 

1. 857 

0.  4000 

1.  887 

0.  5000 

1. 857 
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FIGURE  11-88 

ION  AS  A FUNCTION  OF  CONVERGENCE  RATE  WITH 
ADAPTIVE  FILTER  FOR  A STRONG  SIGNAL  FROM 
LPA,  16  SITES,  1971  DAY  276  21.07.18-21.11.33) 


TABLE  11-45 


SIGNAL -TO-NOISE  GAIN  VERSUS  CONVERGENCE  RATE  WITH 
15-POINT-LONG  ADAPTIVE  FILTER  FOR  A STRONG  SIGNAL  FROM 
KAMCHATKA  (ALPA,  FULL  ARRAY,  1 97 1 DAY  276  21 . 07.  18-21 . 1 1 . 33, 
USING  NOISF  FROM  1972  DAY  335  0415-0815) 


Convergence  Rate 

oo 

Signal-to -Noise  Gain 
(dB) 

0. 0010 

0.  511 

0. 0015 

0.  668 

0.  0020 

0.  792 

0.  0030 

0.  977 

0. 0040 

1.112 

0. 0050 

1.  218 

0. 0070 

1.  376 

0. 0100 

1.  535 

0. 0150 

1.  700 

0.  0200 

1.  87  0 

0. 0300 

1. 971 

0.  0400 

2.  106 

0.  0500 

2.  237 

0.0700 

2.  498 

0. 1000 

2.  892 

0.  1500 

3.  520 

0. 2000 

4.  069 

0. 3000 

4.  923 

0. 4000 

5.  543 

0. 5000 

6.  019 

FIGURE  11-89 

SIGNAL -TO -NOISE  GAIN  AS  A FUNCTION  OF  CONVERGENCE  RATE  WITH 
15-POINT-LONG  ADAPTIVE  FILTER  FOR  A STRONG  SIGNAL  FROM 
KAMCHATKA  (ALPA,  FULL  ARRAY,  1971  DAY  276  21. 07.  18-21.  11.  33, 
USING  NOISE  FROM  1972  DAY  335  0415-0815) 
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visible  in  the  ABF  output  for  p = 0.  5 at  the  bottom.  For  the  second  Kamchat- 
ka event,  which  reached  ALP  A at  approximately  2207,  signal  degradation  with 
the  15 -point  filter  length  is  shown  in  Table  11-46  and  Figure  11-91.  This  event 
has  the  best  signal  similarity  across  channels  of  any  event  in  the  sample  from 
day  276  of  1971.  Table  11-47  and  Figure  II- 92  indicate  the  estimated  signal- 
to-noise  gain  for  the  second  Kamchatka  event.  As  in  the  case  of  the  earlier 
Kamchatka  event,  the  use  of  the  shorter  filter  length  reduces  the  improvement 
relative  to  beamsteering  by  0.  3 to  0.  4 dB  (see  Table  11-38).  Figure  11-93 
pictures  the  beam  outputs  for  this  event.  Signal  distortion  and  degradation 
are  considerably  less  than  for  the  first  Kamchatka  event  in  Figure  11-90. 

Table  11-48  lists  the  power  reduction  relative  to  beamsteering  at 

p = 0.007  and  fx  = 0.  5 with  the  15-point-locig  adaptive  filter  for  four  events 

away  from  the  273°  steer  azimuth  of  the  two  Kamchatka  events.  With  the  ex- 

o 

ception  of  the  first  Panama  event  at  M = 0.  007  and  the  weak  event  from  300  - 
305°  at  p = 0.5,  event  suppression  is  greater  for  the  shorter  filter  length 
than  the  31-point  filter  length  (Table  11-39).  At  the  273°  steer  direction,  off- 
azimuth  event  suppression  again  is  not  degraded  by  the  shorter  filter  length. 
Figure  11-94  exhibits  the  beam  outputs  for  the  weak  event  from  300°-305°. 
Interestingly,  some  energy  is  visible  on  the  adaptive  filter  output  for  4 = 0.  5 
at  the  time  where  signal  degradation  is  worst  for  the  302.5°  steer  azimuth 
(Figure  11-82).  Perhaps  the  apparent  arrival  azimuth  of  this  event  drifts  sig- 
nificantly with  time.  Figure  11-95  is  the  corresponding  beam-output  display 
for  the  Andreanof  Islands  event.  In  the  adaptive  filter  output  for  n = 0.  5 at 
the  bottom,  a cycle  or  so  of  this  event  (36°  from  the  look  azimuth)  escapes 
suppression. 


When  the  new  prefilter  is  applied  to  the  site  vertical-component  data 
before  input  to  the  beamforming  process,  the  New  Britain  event  becomes 
strong  enough  on  the  beam  outputs  to  provide  reasonably  reliable  signal  de- 
gradation measurements.  This  fact  provides  an  opportunity  to  determine 
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TABLE  11-46 


SIGNAL  DEGRADATION  VERSUS  CONVERGENCE  RA  TE  WITH 
15-POINT-LONG  ADAPTIVE  FILTER  FOR  A VERY  STRONG  SIGNAL  FROM 
KAMCHATKA  (ALPA,  16  SITES,  1971  DAY  276  22.07.02-22.11.17) 


Convergence  Rate 
(M) 

Signal  Degradation 
(dB) 

0. 0010 

0.  105 

0. 0015 

0.  150 

0.  0020 

0.  191 

0. 00  30 

0.  264 

0. 0040 

0.  328 

0.  0050 

0.  386 

0. 0070 

0.  484 

0.0100 

0.  600 

0.0150 

0.  742 

0.  0200 

0.  841 

0.  0300 

0.  962 

0. 0400 

1. 024 

0.  0 500 

1. 056 

0. 0700 

1. 088 

0. 1000 

1.  107 

0.  1500 

1. 090 

0. 2000 

1.019 

0.  3000 

0.  782 

0. 4000 

0.  486 

0. 5000 

0.  175 
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TABLE  11-47 


SIGNAL -TO -NOISE  GAIN  VERSUS  CONVERGENCE  RATE  WITH 
15-POINT-LONG  ADAPTIVE  FILTER  FOR  A VERY  STRONG  SIGNAL  FROM 
KAMCHATKA  (ALPA,  FULL  ARRAY,  1971  DAY  276  22.07.02-22.11.17, 
USING  NOISE  FROM  1972  DAY  335  0415-0815) 


E 


Comvergence  Rate 


0.  0010 


0.  0015 


0.  0020 


0.  0030 


0. 0040 


0.  0050 


0. 0070 


0. 0100 


0. 0150 


0.  0200 


0.  0300 


0. 0400 


0.  0500 


0. 0700 


0. 1000 


0.  1500 


0. 2000 


0. 3000 


0. 4000 


0. 5000 


Signal-to-Noise  Gain 
(dB) 


11-172 


0.  490 


0.639 


0.  755 


0.  929 


1.055 


1.  153 


1.  303 


1.  464 


1. 652 


1.  797 


2.045 


2.  273 


2.491 


2.  886 


3.  404 


4.  158 


4.  834 


5.998 


6.  944 


7.  701 


m 


■ 

■ 
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Gain  (dB) 


Rl 


Si 


Convergence  Rate  M 


FIGURE  11-92 

SIGNAL  - TO -NOISE  GAIN  AS  A FUNCTION  OF  CONVERGENCE  RATE  WITH 
15-POINT-LONG  ADAPTIVE  FILTER  FOR  A VERY  STRONG  SIGNAL  FROM 
KAMCHATKA  (ALPA,  FULL  ARRAY,  1 97 1 DAY  276  22.  07.  02-22.  1 1 . 17, 
USING  NOISE  FROM  1972  DAY  335  0415-0815) 


m m 

. N 

o m 


f-H  h- 

00  O'  00 

in  0s  O' 

in  o ro 


II 


rt 

W 

Pu 

D 

oo  in 


in  _ f-i 
w Z h 
w O M 
b;  rH  w 

n ^ ^ 
l-M  f 1 t— < 

PU  H H 


>-'  fV*  ^ 

in  « < 
H P P 
2 P5  < 

w w o 
> w z 
who 

e w J 
ho  ' 

on  H 

d n-  Z 
5 N S 

2 g 8 

< p ' 

i W in 

w ^ 

w x 

o H 


o <->  ,-t  m 

N N N N 


rt  rt 

B B 

rt  rt 
C C 
rt  rt 

P P 


BEAMS  FOR  A WEAK  EVENT  FROM  300  -305  WITH  1 5 -POINT-LONG 
ADAPTIVE  FILTER  (DAY  276  1971,  STEER  DIRECTION  273°) 
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whether  the  inferior  signal-to-noise  gain  for  the  weak  event  from  300°-305° 
is  due  to  poor  signal  similarity  or  low  signal-to-noise  ratio.  Table  11-49 
and  Figure  11-96  portray  signal  degradation  for  the  New  Britain  event  as  a 
function  of  the  convergence  parameter  n . The  degradation  is  significantly 
lower  than  for  the  weak  event  from  300°-305°,  so  that  poor  signal  similarity 
is  the  apparent  reason  for  the  results  with  that  event.  Table  11-50  and  Figure 
11-97  present  the  estimated  signal-to-noise  gain  for  the  New  Britain  event  us- 
ing the  noise  sample  from  day  335  of  1972.  The  6.  6 dB  gain  at  H = 0.  5 is 
even  higher  than  the  corresponding  6.  0 dB  gain  for  the  first  Kamchatka  event, 
which  is  considerably  stronger.  Figure  11-98  displays  the  beamsteer  output 
and  adaptive  beam  outputs  at  M = 0.  007  and  \i  - 0.  5 for  the  New  Britain  event 
with  the  new  prefilter  and  an  adaptive -filter  length  of  fifteen  points.  A clear 
detection  is  evident  on  the  adaptive  beam  output  for  n « 0.  5 in  the  bottom 
trace.  To  demonstrate  the  effectiveness  of  the  new  prefilter  relative  to  the 
old  prefilter  in  detecting  signals.  Figure  11-99  pictures  the  beamsteer  outputs 
and  ABF  outputs  at  n = 0.  005  and  H - 0.  5 for  the  New  Britain  event  with  the 
old  prefilter.  A"  the  maximum  convergence  rate  in  the  bottom  trace,  the 
event  is  barely  visible.  The  difference  in  signal-to-noise  ratio  between  the 
new  prefilter  and  the  old  prefilter  seems  to  be  about  6 dB  on  the  adaptive 
beams  at  H = 0.  5. 

Five  events  in  the  data  sample  from  day  276  of  1971  arrive  from  azi- 
muths away  from  the  240°  azimuth  of  the  New  Britian  event.  Table  11-51  re- 
cords the  power  reduction  relative  to  beamsteering  with  the  new  prefilter  for 
these  events.  In  general,  the  off-azimuth  event  suppression  is  considerably 
less  at  the  240°  look  azimuth  than  at  the  302.  5°  and  273°  steer  directions.  One 
contributing  factor  is  the  excellent  240  time-shift-and-sum  beam  pattern, 
which  appears  in  Figure  II— 100.  The  potential  for  improvement  with  adaptive 
amforming  is  substantially  limited  in  this  situation.  Figure  II- 1 01  presents 
the  beam  outputs  for  the  first  Panama  event.  At  H = 0.  007,  the  adaptive  beam 
is  unable  to  suppress  the  event  as  well  as  the  time-shift-and-sum  beam  in  the 
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TABLE  11-49 

SIGNAL  DEGRADATION  VERSUS  CONVERGENCE  RATE  WITH 
15-POINT-LONG  ADAPTIVE  FILTER  FOR  A WEAK  SIGNAL  FROM 
NEW  BRITAIN  (ALPA,  16  SITES,  1971  DAY  276  23.  10.  30-23.  18.  00) 


Convergence  Rate 
00 

Signal  Degradation 
(dB) 

0. 0010 

0.  079 

0. 001 5 

o.  no 

0. 0020 

0.  134 

0. 0030 

0.  171 

0. 0040 

0.  195 

0.  00  50 

0.  211 

0. 0070 

0.  228 

0. 0100 

0.  234 

0.  0150 

0.  226 

0.  0200 

0.  212 

0.  0300 

0.  187 

0.  0400 

0.  170 

0.  0500 

0.  161 

0. 0700 

0.  164 

0. 1000 

0.  206 

0.  1 500 

0.  323 

0. 2000 

0.  454 

0.  3000 

0.  715 

0.  4000 

0.993 

0.  5000 

1.  304 

I 


TABLE  11-50 

SIGNAL -TO -NOISE  GAIN  VERSUS  CONVERGENCE  RATE  WITH 
15-POINT-LONG  ADAPTIVE  FILTER  FOR  A WEAK  SIGNAL  FROM 
NEW  BRITAIN  (ALPA,  FULL  ARRAY,  1971  DAY  276  23.  10.  30-23.  18.  00 
USING  NOISE  FROM  1972  DAY  335  0415-0815) 


I* 


Convergence  Rate 
(M) 

Signal-to-Noise  Gain 
(dB) 

0.0010 

0.  516 

0.  0015 

0.  679 

0. 0020 

0.  812 

0.0030 

1. 022 

0. 0040 

1.  188 

0.  0050 

1.  328 

0.0070 

1.559 

0. 0100 

1. 830 

0. 0150 

2.  168 

0.  0200 

2.  426 

0.0300 

2.  820 

0.0400 

3.  127 

0.0500 

3.  386 

0. 0700 

3.  810 

0. 1000 

4.  305 

0. 1500 

4.  925 

0. 2000 

5.  399 

0. 3000 

6.065 

0. 4000 

6.  437 

0. 5000 

6.  572 

se 


FIGURE  11-97 

SIGNAL -TO-NOISE  GAIN  AS  A FUNCTION  OF  CONVERGENCE  RATE 
15-POINT-LONG  ADAPTIVE  FILTER  FOR  A WEAK  SIGNAL  FROM  NEW 
(ALPA,  FULL  ARRAY,  1971  DAY  276  23.  10.  30-23.  18.00,  USING  NOIS 

1972  DAY  335  0415-0815) 
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FULL-ARRAY  BEAMS  FOR  NEW  BRITAIN  EVENT  WITH 
nr,n  PREFILTER  (DAY  276  1971,  STEER  DIRECTION  240 


T I ME-SH I FT-AND-SUM  BEAM  PATTERN 
BEAM  LOOK  VELOCITY  IS  3,6,  LOOK  AZIMUTH  240,0 
FREQUENCY  IS  0,04000  HZ,  PERIOD  25,0  SECONDS 


FIGURE  II- 100 

BEAMSTEER  RESPONSE  FOR  FULL-ARRAY 
DATA  SAMPLE  FROM  DAY  276  OF  1971 
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FIGURE  II- 101 

BEAMS  FOR  A MAGNITUDE  3.6  PANAMA  EVEN.T  WITH  15-POINT-LONG 
ADAPTIVE  FTLTER  (DAY  276  1971,  STEER  DIRECTION  240°) 


time  interval  2036-2039.  Figures  11-102  through  11-104  show  the  ABF  beam 
patterns  at  2033,  2036  1/2,  and  2040.  The  ABF  response  for  113°  at  25  sec- 
onds is  initially  greater  than  the  beams  tee  r response  and  does  not  drop  until 
2040.  Unlike  the  adaptive  beam  for  P = 0.  007,  the  adaptive  beam  for  p = 0.  5 
virtually  annihilates  the  first  Panama  event.  Figures  11-105  through  11-108 
display  the  remaining  four  off-azimuth  events.  In  the  case  of  the  two  Kamchatka 
events  33°  from  the  steer  direction,  considerable  low-frequency  energy  at  the 
beginning  of  the  events  leaks  into  the  adaptive  beams  at  the  highest  convergence 
rate.  The  beam  outputs  for  the  second  Panama  event  are  similar  to  those  for 
the  first  Panama  event,  although  event  suppression  is  better  for  the  second. 

After  the  event  has  subsided  on  the  beamsteer  output,  it  persists  on  the  adap- 
tive beam  output  for  a short  time  after  212’.5  at  the  slower  convergence  rate. 

For  the  weak  event  from  300  -305°,  event  suppression  relative  to  beamsteer- 
ing is  worse  at  the  240°  steer  direction  than  at  the  273°  steer  direction.  The 
probable  explanation  is  the  improved  ability  of  the  time-shift-and-sum  beam- 
former  to  attenuate  the  event  when  the  look  azimuth  is  240°.  Compared  with 
the  ABF  outputs  at  the  273°  steer  azimuth  (Figure  11-94),  the  adaptive  beams 
at  240°  seem  to  be  slightly  weaker  for  this  event. 

The  Andreanof  Islands  event  from  2 37°  is  only  three  degrees  from  the 
240  steer  direction.  Figure  11-109  is  a Calcomp  plot  of  the  beamsteer  and 
ABF  outputs  at  p = 0.007  and  p = 0.  5 for  this  event.  At  the  slower  conver- 
gence rate,  the  adaptive  beam  output  closely  resembles  the  beamsteer  output. 

At  the  maximum  convergence  rate,  however,  the  ABF  output  is  3.4  dB  weak- 
er than  the  beamsteer  output  in  the  top  trace.  Since  a three-degree  shift  in 
azimuth  has  only  a minimal  effect  on  the  steer  delays,  it  is  extremely  unlikely 
that  the  3.4  dB  attenuation  is  due  to  the  high  resolution  of  adaptive  beamform- 
ing. More  likely  explanations  are  poor  signal  similarity  across  the  array  and 
an  apparent  arrival  azimuth  different  from  the  true  237°  azimuth  of  the  event. 

A visual  examination  of  the  site  vertical-component  data  reveals  that  the 
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flLPR 

ABF  BERM  PATTERN  AT  276/20.33.00 
CONVERGENCE  RATE  IS  0,0070 
BEAM  LOOK  VELOCITY  IS  3.6.  LOOK  AZIMUTH  240.0 
FREQUENCY  IS  0.04000  HZ.  PERIOD  25.0  SECONDS 


FIGURE  11-102 

FULL-ARRAY  ABF  RESPONSE  FOR  MAGNITUDE 
3.6  PANAMA  EVENT  AT  20  33  ON  DAY  276  OF  1971 


ABF  BEAM  PATTERN  AT  276/20,36,30 
CONVERGENCE  RATE  IS  0,0070 
BEAM  LOOK  VELOCITY  IS  3,6,  LOOK  AZIMUTH  240,0 
FREQUENCY  IS  0,04000  HZ,  PERIOD  25,0  SECONDS 


FIGURE  11-103 

FULL-ARRAY  ABF  RESPONSE  FOR  MAGNITUDE 
3.6  PANAMA  EVENT  AT  20  36  1/2  ON  DAY  276  OF  1971 
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ABF  BEAM  PATTERN  AT  276/20.40.00 
CONVERGENCE  RATE  IS  0.0070 
BEAM  LOOK  VELOCITY  IS  3.6.  LOOK  AZIMUTH  240,0 
FREQUENCY  IS  0.04000  HZ.  PERIOD  25.0  SECONDS 


FIGURE  II-  104 

FULL-ARRAY  ABF  RESPONSE  FOR  MAGNITUDE 
3.6  PANAMA  EVENT  AT  2040  ON  DAY  276  OF  1971 
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BEAMS  FOR  MAGNITUDE  4.5  KAMCHATKA  EVENT  WITH  15-POINT-LONG 
ADAPTIVE  FILTER  (DAY  276  1971,  STEER  DIRECTION  240°). 


Adaptive  Filter  Output  (/i-  0.  5) 

FIGURE  11-106 

BEAMS  FOR  MAGNITUDE  4.7  PANAMA  EVENT  WITH  15 -POINT-LONG 
ADAPTIVE  FILTER  (DAY  276  1971,  STEER  DIRECTION  240°) 


BEAMS  FOR  MAGNITUDE  4.9  KAMCHATKA  EVENT  WITH  15 -POINT -LONG 
ADAPTIVE  FILTER  (DAY  276  1971,  STEER  DIRECTION  240°) 


I 

I 
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signal  similarity  is  not  perfect,  but  not  much  worse  than  that  of  the  two  Kam- 
chatka events.  Perhaps  a combination  of  less-than-perfect  signal  similarity 
and  a discrepancy  between  the  source  azimuth  and  arrival  azimuth  can  account 
for  the  partial  suppression  of  this  event. 

The  results  of  this  subsection  clearly  demonstrate  the  effectiveness  of 
the  new  prefilter  depicted  in  Figure  11-71.  Detection  improvement  from  the 
new  prefilter  is  not  limited  to  that  achieved  by  frequency  filtering.  The  new 
prefilter's  emphasis  on  the  more  coherent  energy  at  low  frequencies  raises 
adaptive  filtering  noise  reduction  relative  to  beamsteering  by  more  than  3 dB 
compared  with  the  old  prefilter  of  Figure  II- 1 . This  additional  noise  reduction 
occurs  when  vertical-component  data  passed  through  the  new  prefilter  are  in- 
put to  an  adaptive  beamformer  having  either  31  or  15  points  per  channel.  The 
difference  in  noise  reduction  for  the  two  adaptive  filter  lengths  is  between  0.  3 
and  0.4  dB  with  the  new  prefilter.  Signal  degradation  for  the  weak  event  from 
300°-305°,  however,  is  more  than  2.  1 dB  less  with  the  shorter  filter  length 
For  this  event,  which  has  poor  signal  similarity,  the  signal-to -noise  gain  is 
more  than  1.8  dB  greater  with  the  shorter  adaptive  filter  length.  For  the  two 
Kamchatka  events  with  good  similarity,  the  difference  in  signal  degradation  be- 
tween the  two  adaptive  filter  lengths  is  less  than  0.  1 dB.  The  loss  in  estimated 
signal-to-noise  gain  is  between  0.  3 and  0.4  dB  with  the  shorter  filter  length. 

The  significantly  reduced  computation  load  associated  with  the  shorter  adaptive- 
filter  length  probably  justifies  the  use  of  the  1 5-point-long  filter  in  any  on-line 
implementation  of  adaptive  beamforming.  Estimated  signal-to-noise  gain  re- 
lative to  beamsteering  for  the  New  Britain  event  is  6.6  dB  at  the  maximum 
convergence  rate  with  the  new  prefilter  and  the  1 5 -point-long  adaptive  filter. 

The  highest  signal-to-noise  gains  were  obtained  at  //  = 0.  5 for  the  New  Britain 
event  and  the  other  three  events  used  for  signal  degradation  measurements  with 
the  new  prefilter.  Using  the  1 5 -point-long  adaptive  filter,  the  estimated  signal- 
to-noise  gains  were  3.3  dB,  6.0  dB,  and  7.7  dB,  respectively,  for  the  weak 
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event  from  300  -305  , the  first  Kamchatka  event,  and  the  second  Kamchatka 
event.  Thus  the  signal-to-noise  ratio  improvement  relative  to  beamsteering 
ranged  from  3 to  8 dB  for  the  four  events  examined.  An  improvement  of  6 dB 
can  apparently  be  expected  under  ordinary  conditions.  This  6 dB  improvement 
over  beamsteering  occurs  in  addition  to  detection  gains  obtained  from  the  fre- 
quency response  of  the  new  prefilter. 


G.  SUMMARY 

The  new  adaptive  algorithm 

M t +At)  = A + A(X:X) 

xTx 

yields  dramatic  signal-to-noise  ratio  improvement  relative  to  the  old  adaptive 
algorithm 

A(  t + A t)  = A + 

(X-X)T(X-X) 

employed  in  the  previous  adaptive -processing  study.  With  the  old  adaptive 
algorithm,  signal  degradation  increases  as  the  signal-to-noise  ratio  increases 
and  the  maximum  likelihood  constraints  do  not  guarantee  preservation  of  the 
signal.  With  the  new  adaptive  algorithm,  on  the  other  hand,  signal  degradation 
seems  to  depend  more  on  signal  similarity  across  the  array  than  on  signal-to- 
noise  ratio,  so  that  the  two  Kamchatka  events  are  well-preserved,  although  in 
distorted  form,  even  at  the  maximum  convergence  rate.  For  the  weak  event 
from  300  -305  , which  has  poor  signal  similarity  between  sites,  the  estimated 
signal-to-noise  gain  based  on  the  noise  sample  from  day  238  of  1970  is  never 
much  more  than  one  dB  higher  with  the  new  adaptive  algorithm  than  with  the  old 
algorithm  when  six  sites  and  the  old  prefiltcr  of  Figure  II- 1 are  employed  in  the 
beamforming  process.  In  the  case  of  the  weak  event  from  300°-305°,  highest 


signal-to-noise  improvement  relative  to  beamsteering  occurs  with  both  adap- 
tive algorithms  at  slow  convergence  rates  near  0.004  and  g = 0.  005.  The 
essential  difference  between  the  two  adaptive  algorithms  emerges  in  the  pro- 
cessing results  for  the  two  Kamchatka  events.  With  the  old  adaptive  algorithm, 
the  greatest  signal-to-noise  gains  of  0.  965  dB  and  0.  136  dB  arise  at  p = 0.  004 
and  p = 0.  002,  respectively,  for  the  first  and  second  Kamchatka  events.  At 
p = 0.25,  processing  gains  relative  to  beamsteering  are  either  minor  (for  the 
first  Kamchatka  event)  or  negative  (for  the  second  Kamchatka  event).  With 
the  new  adaptive  algorithm,  however,  the  best  signal-to-noise  gains  for  the 
two  Kamchatka  events  are  at  the  maximum  convergence  rate.  At  p = 0.  5,  the 
estimated  signal-to-noise  ratio  improvements  based  on  the  noise  sample  from 
day  238  of  1970  are  5.75  dB  and  7.5  dB,  respectively,  for  the  first  and  second 
Kamchatka  events  with  six  sites  and  the  old  prefilter.  In  detection  processing, 
where  signal  distortion  is  a relatively  minor  consideration,  these  sharply  im- 
proved gains  demonstrate  that  significantly  improved  detection  capability  for 
events  with  good  signal  similarity  is  possible  at  rapid  convergence  rates.  An- 
other effect  of  operating  at  high  convergence  rates  is  a narrowing  of  the  adap- 
tive processor  beamwidth,  as  indicated  by  the  greater  suppression  of  off- 
azimuth  events  at  p = 0.  5 than  at  lower  values  of  the  convergence  parameter. 

The  effectiveness  of  the  new  adaptive  algorithm  hinges  on  slowing  the  effective 
adaptation  rate  in  the  presence  of  on-azimuth  signals.  In  the  extreme  case  of 
an  ideal  signal  (where  X = X),  the  filter  vector  A stops  completely  with  the  new 
algorithm,  whereas  the  step  size  | A(t  + \t)  - A(t)  | becomes  infinite  with  the 
old  algorithm.  Consequently,  replacing  the  old  adaptive  algorithm  with  the  new 
adaptive  algorithm  reduces  signal  degradation,  particularly  for  events  with  good 
signal  similarity. 

When  beams  are  formed  from  data  passed  through  the  old  prefilter,  the 
estimated  signal-to-noise  gain  of  adaptive  beamforming  relative  to  beamsteer- 
ing is  normally  lower  with  the  full  array  than  with  six  channels.  Differences 


in  measured  signal  degradation  between  the  full  array  and  the  closely-spaced 
six-channel  array  are  highly  variable,  particularly  at  the  higher  convergence 
rates,  where  2 dB  differences  occur  for  the  weak  event  from  300°-305°  and 
the  second  Kamchatka  event.  The  main  reason  for  the  reduced  improvement 
relative  to  beamsteering  with  the  full  array,  however,  is  not  the  change  in 
measured  signal  degradation  but  instead  the  measured  noise  reduction  for  the 
data  sample  from  day  335  of  1972:  at  all  convergence  rates,  the  noise  reduc- 
tion is  lower  with  the  full  array  than  with  six  channels.  One  of  the  reasons 
for  the  lower  full-array  noise  reduction  is  that  the  more  widely-spaced  sen- 
sors of  the  full  array  have  less  coherence  between  sensor  pairs  than  the  close- 
ly packed  elements  of  the  array  used  for  the  six-channel  measurements.  An- 
other reason  is  that,  for  the  arrays  used  in  the  noise  reduction  measurements, 
the  ratio  of  the  beamsteer  directivity  index  for  the  fifteen-channel  array  to  that 
for  the  six-channel  array  is  greater  than  15/6  ( = 4 dB).  Thus  the  ratio  of  co- 
herent energy  to  spatially  uncorrelated  energy  in  the  beamsteer  output  is  less 
for  the  fifteen-channel  array  than  for  the  six-channel  array  (assuming  coherent 
energy  is  equally  likely  at  all  azimuths).  As  a result,  the  potential  noise  re- 
duction of  adaptive  beamforming  relative  to  beamsteering  is  less  for  fifteen 
channels  than  for  six  channels. 

The  new  signal-enhancement  prefilter  shown  in  Figure  11-71  is  a definite 
improvement  over  the  old  prefilter.  Not  only  is  detection  gain  achieved  by  the 
frequency  response  of  the  new  prefilter,  but  also  the  new  prefilter's  emphasis 
on  the  more  coherent  energy  at  low  frequencies  increases  ABF  noise  reduc- 
tion relative  to  beamsteering  by  more  than  3 dB  compared  with  the  old  prefilter. 
This  additional  noise  reduction  occurs  with  both  the  31 -point  and  15-point  adap- 
tive filter  lengths.  The  difference  in  noise  reduction  with  the  new  prefilter  is 
only  0.  3 - 0.4  dB  for  the  two  adaptive-filter  lengths.  For  the  two  Kamchatka 
events  with  good  signal  similarity,  the  difference  in  signal  degradation  is  less 
than  0.  1 dB,  so  that  the  loss  in  estimated  signal-to-noise  gain  due  to  reduced 
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noise  reduction  is  0.  3 - 0.  4 dB  with  the  shorter  filter  length.  On  the  other 
hand,  signal  degradation  for  the  weak  event  from  300°-305°,  which  has  poor 
signal  similarity,  is  more  than  2.  1 dB  lower  with  the  1 5 -point-long  adaptive- 
filter  length,  so  that  the  estimated  signal-to-noise  gain  is  actually  more  than 
1.8  dB  higher  with  the  shorter  filter  length.  Estimated  signal-to-noise  gain 
relative  to  beamsteering  for  the  New  Britain  event  is  almost  6.  6 dB  at  the 

maximum  convergence  rate  with  the  new  prefilter  and  the  15 -point-long  adap- 
tive filter. 


Table  11-52  summarizes  the  estimated  signal-to-noise  gains  relative  to 
beamsteering  for  on-azimuth  events  processed  with  the  new  adaptive  algorithm. 
Noise  reduction  estimates  for  this  summary  table  are  based  in  all  cases  on  the 
noise  sample  from  0415  to  0815  on  day  335  of  1972.  Both  of  the  blocks  showing 
results  for  the  new  prefilter  apply  to  ALP  A data  in  which  all  available  sites 
were  processed.  With  the  exception  of  the  weak  event  from  300°-305°  when 
processed  with  the  old  prefilter  using  six  sites,  the  best  signal-to-noise  ratio 
improvement  relative  to  beamsteering  occurs  at  the  most  rapid  convergence 
rate.  Evidently  the  adaptive  processor  can  take  advantage  of  transient  correla- 
tion between  successive  data  vectors  X in  such  a way  that  noise  reduction  in- 
creases more  than  signal  degradation  as  the  convergence  parameter  p approach- 
es its  maximum  value  0.  5.  When  vertical-component  data  passed  through  the 
old  prefilter  form  the  beam  inputs,  signal-to-noise  gain  is  normally  less  with 
the  full  array  than  with  six  channels  because  of  the  lower  noise  reduction  ob- 
tained with  the  more  widely- spac ed  sites  of  the  full  array.  The  new  prefilter 
provides  increased  signal-to-noise  ratio  improvement  relative  to  beamsteering 
over  and  above  any  detection  gain  from  its  frequency  response.  The  use  of  a 
15-point-long  adaptive  filter  instead  of  a 3 1 -point-long  adaptive  filter  increases 
the  signal- to-noise  ratio  for  the  weak  event  from  300°-305°  more  than  1.8  dB 
by  reducing  the  signal  degradation  for  this  event,  which  has  poor  signal  simi- 
larity between  sites.  At  the  same  time,  however,  the  lower  noise  reduction 


TABLE  11-52 

SUMMARY  OF  ESTIMATED  SIGNAL- TO-NOISE  RATIO  IMPROVEMENT 
RELATIVE  TO  BEAMSTEERING  AT  SELECTED  CONVERGENCE 
PARAMETERS  ^ UNDER  VARIOUS  CONDITIONS 


Signal-to-Noise  Gain  in  dB  with  6 Sites  and  Old  Prefilter 


Event 

U=  0.  004 

0.  007 

M = 0.  5 

Weak  Event  from  300°-305° 

1.488 

1.  538 

-1.  120 

First  Kamchatka  Event 

1.254 

1.  522 

4.  079 

Second  Kamchatka  Event 

1. 208 

1.  502 

5.  874 

Signal-to-Noise  Gain  in  dB  with  Full  Array  and  Old  Prefilter 

Event 

H = 0.  5 

Weak  Event  from  300°-305° 

0.  263 

0.  245 

1.  149 

First  Kamchatka  Event 

0.  804 

1.  013 

4.  646 

Second  Kamchatka  Event 

0.  570 

0.  681 

3.  846 

Signal-to-Noise  Gain  in  dB  with  New  Prefilter  and  31 -Point  ABF 

Event 

H=  0.  007 

/i=  0.  5 

Weak  Event  from  300°-305° 

0.  903 

1.  512 

First  Kamchatka  Event 

1. 464 

6.  358 

Second  Kamchatka  Event 

1.  350 

8.  067 

Signal-to-Noise  Gain  in  dB  with  New  Prefilter  and 

1 5-Point  ABF 

Event 

H = 0.  007 

H = 0.  5 

Weak  Event  from  300°-305° 

1. 026 

3.  338 

First  Kamchatka  Event 

1.  376 

6.  019 

Second  Kamchatka  Event 

1.  303 

7.  701 

New  Britain  Event 

1.  559 

6.  572 
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with  the  shorter  filter  length  results  in  decreased  signal-to-noise  gain  for  the 
two  Kamchatka  events,  which  have  good  signal  similarity.  The  significantly 
diminished  computational  load  for  the  shorter  adaptive-filter  length  probably 
justifies  the  use  of  the  1 5 -point-long  filter  in  any  on-line  implementation  of 
adaptive  beamforming.  In  any  event,  the  signal-to-noise  ratio  improvement 
for  the  full  array  with  the  new  prefilter  and  an  adaptive -filter  length  of  1 5 
points  per  channel  falls  within  a range  of  3 to  8 dB.  Normally,  a gain  of  6 dB 
might  be  expected  in  this  mode  of  operation.  Thus  adaptive  beamforming 
could  conceivably  lower  the  Rayleigh-wave  detection  threshold  at  ALPA  by  0.  3 
magnitude  units. 

Off-azimuth  events  are  strongly  suppressed  by  adaptive  beamforming 
at  the  maximum  convergence  rate,  where  the  greatest  improvement  over  time- 
shift-and-sum  beamforming  is  achieved  for  on-azimuth  events.  Some  off- 
azimuths  signals  are  almost  completely  annihilated.  Ordinarily,  the  strongest 
off-azimuth  events  are  attenuated  the  most.  For  some  signals  30°  to  60°  from 
the  steer  direction,  however,  considerable  energy  leaks  into  the  adaptive  beams, 
even  in  the  case  of  the  strongest  events.  The  use  of  the  new  signal-enhancement 
prefilter  produces  some  reduction  in  off-azimuth  event  suppression.  The  ap- 
parent explanation  for  this  result  is  the  emphasis  given  to  low-frequency  energy 
by  the  new  prefilter.  This  feature  not  only  increases  the  effective  beamwidth 
of  the  adaptive  processor  but  also  concentrates  more  energy  near  the  beginning 
of  events,  so  that  the  adaptive  filter  set  cannot  react  as  effectively  to  the  pre- 
sence of  off-azimuth  events.  Even  with  the  new  prefilter,  however,  the  atten- 
uation of  off-azimuth  events  is  considerable.  This  fact  demonstrates  that  great- 
er directional  resolution  at  ALPA  and  other  similar  long-period  arrays  is  pos- 
sible through  the  beam- narrowing  capability  of  adaptive  processing.  Such  ar- 
rays, in  effect,  become  superdirective  antennas  when  adaptive  beamforming  is 
employed. 
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SECTION  III 


INTERFERING-EVENT  RESULTS 

A.  DISCUSSION 

The  previous  section  of  this  report  focused  on  the  ability  of  adaptive 
beamforming  to  detect  individual  unmixed  seismic  events  in  the  presence  of 
background  noise.  This  section  presents  adaptive  processing  results  for  mix- 
ed events  in  which  both  an  on-azimuth  signal  and  an  off-azimuth  interfering 
event  are  present  in  time-overlapped  form  on  the  data  traces.  With  data  as 
actually  recorded,  there  is  no  reliable  method  to  isolate  the  contribution  from 
each  separate  event  on  the  beam  output.  The  most  realistic  method  for  eval- 
uating the  effect  of  adaptive  multichannel  filtering  on  mixed  events  seems  to  be 
a simulation  where  two  data  samples,  each  containing  one  signal,  are  summed 
to  create  a composite  sample  containing  an  interfering  event.  Such  a simula- 
tion technique  permits  considerable  variety  in  the  situations  which  can  be  ex- 
amined. The  relative  strength  of  the  events,  the  time  interval  between  their 
arrivals,  and  their  azimuthal  separation  are  all  subject  to  control  with  the  pro- 
cedure described.  In  assessing  the  potential  detection  gains  of  adaptive  beam- 
forming over  conventional  time-shift-and-sum  processing,  the  beam  outputs 
for  the  composite  data  sample  are  the  most  important  performance  criterion. 
As  long  as  something  closely  resembling  a seismic  event  appears  on  the  com- 
posite sample  beam  output  during  the  time  period  when  the  on-azimuth  signal 
is  present,  a detection  can  be  declared.  Since  only  the  composite-sample 
beam  is  available  in  practice,  phenomena  on  the  beam  outputs  for  the  two  in- 
dividual samples  can  only  shed  light  on  the  processes  involved  in  adaptive 
beamforming,  but  cannot  serve  as  the  deciding  factor  in  its  evaluation.  The 
relative  on-azimuth  to  off-azimuth  event  strength  at  which  a detection  is 
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possible  on  the  beamsteer  or  ABF  output  for  the  composite  sample  determines 
the  detection  gain  of  adaptive  multichannel  filtering  over  time-shift-and-sum 
be amf arming.  A second  but  still  important  consideration  in  adaptive  process- 
ing performance,  once  a signal  is  detected,  is  the  r.ccuracy  with  which  the 
peak-to-peak  amplitude  of  the  composite-sample  beam  reflects  the  peak-to- 
peak  amplitude  of  the  on-azimuth  event.  Measurement  of  this  amplitude  fur- 
nishes a surface- wave  magnitude  estimate  for  the  detected  on- azimuth  event. 

As  mentioned  earlier,  the  basic  mixed-event  simulation  procedure  sums 
two  data  samples,  one  with  an  on-azimuth  signal  and  another  with  an  off-azimuth 
event,  to  create  a composite  data  sample  containing  an  interfering  event.  Fig- 
ure III- 1 is  a schematic  diagram  of  the  interfering-event  simulation  procedure. 
Before  the  summation  which  creates  the  composite  data  sample  used  in  design- 
ing the  adaptive  filter  set,  individual  sites  of  the  sample  with  the  on-azimuth 
signal  are  individually  scaled,  passed  through  a zero-phase  prefilter  identical 
for  all  sites,  and  individually  time-shifted  to  time-align  the  on-azimuth  signal 
at  all  sites.  In  the  case  of  the  sample  with  the  off-azimuth  event,  preprocess- 
ing of  the  individual  sites  is  the  same  except  that,  befo  *e  the  on-azimuth  signal- 
alignment  time  delays  are  applied,  the  off-azimuth  event  can  be  time-shifted 
from  its  true  direction  so  as  to  appear  to  come  from  a different  direction. 

Thus  the  individual  sites  of  the  off-azimuth  event  sample  seem  to  be  time- 
shifted  so  that  the  tin.-'  delays  for  the  on- azimuth  signal  are  applied  to  an  event 
coming  from  the  specified  fictitious  arrival  direction.  Immediately  before  the 
adaptive -filter  update,  the  two  data  samples  are  summed  so  that  the  new  adap- 
tive filter  update  algorithm 
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utilizes  the  composite  data  sample  for  adaptive  filter  design.  To  eliminate  re- 
dundant computations,  the  computer  program  used  in  this  study  can  simultane- 
ously update  twenty  adaptive  filter  sets,  each  corresponding  to  ore  particular 
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FIGURE  III-  1 

BLOCK  DIAGRAM  OF  INTERFERING-EVENT  SIMULATION  PROCEDURE 


value  of  the  convergence  parameter  /!  . To  observe  the  effect  of  adaptive 
beamforming  on  each  of  the  two  data  samples  used  to  form  the  composite  sam- 
ple, the  adaptive  filter  sets  are  applied  separately  to  each  sample.  Subsequent- 
ly, the  beam  outputs  for  both  samples  are  summed  to  create  the  composite- 
sample  beam  output  for  each  adaptive  filter  set.  Similarly,  time-shift-and- 
sum  outputs  are  computed  for  both  data  samples  and  later  added  to  produce 
the  comp-' site-sample  beamsteer  output.  The  composite-sample  beam  output, 
for  beamsteering  or  for  any  adaptive  filter  set,  represents  the  only  beam  out- 
put available  in  the  simulated  real-world  situation.  However,  to  facilitate  an 
understanding  of  the  processes  involved  in  adaptive  beamforming  (or  in  beam- 
steering), Calcomp  plots  of  the  individual-sample  beam  outputs  are  displayed 
along  with  the  composite-sample  beam  output  for  any  selected  time-shift-and- 
sum  or  adaptive  beamformer.  The  series  of  three-trace  plots  produced  in 
this  way  form  the  principal  basis  for  evaluating  adaptive-processing  perform- 
ance relative  to  beamsteering.. 

The  results  presented  in  this  section  are  preliminary  and  encompass 
only  a limited  range  of  the  situations  possible  in  practice.  Subsection  B dis- 
cusses adaptive-processing  results  obtained  using  long -period  data  from 
ALP  A,  This  subsection  employs  three  different  pairs  of  interfering  events  to 
compare  the  relative  performance  of  adaptive  beamforming  and  time-shift-and- 
sum  processing.  In  one  case,  the  data  sample  with  the  off-azimuth  event  is 
formed  from  the  on-azimuth  event  sample  by  shifting  the  apparent  arrival 
azimuth  of  the  on-azimuth  signal  by  180°.  The  other  two  event  combinations 
utilize  two  distinct  data  samples.  Subsection  C is  a brief  discussion  of  inter- 
fering event  results  based  on  data  from  the  Korean  short-period  array.  Sub- 
section D summarizes  the  tentative  long -period  interfering -event  results  given 
in  Subsection  B. 
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B.  ADAPTIVE  BEAMFORMING  APPLIED  TO  LONG-PERIOD  INTERFER- 
ING EVENTS  AT  ALPA 

The  principal  goal  of  this  subsection  is  to  examine  the  detection  improve- 
ment of  adaptive  beamforming  relative  to  beamsteering  for  long-period  interfer- 
ing events  simulated  by  adding  two  suitably  preprocessed  data  samples  from  the 
Alaska  Long-Period  Array.  A secondary  objective  is  to  evaluate  the  reliabil- 
ity of  the  peak-to-peak  amplitude  on  the  composite-sample  adaptive  filter  out- 
put as  an  estimate  of  the  on-azimuth  signal's  peak-to-peak  amplitude.  Through- 
out this  subsection,  the  vertical  components  for  the  individual  sites  pass  through 
the  new  prefilter  shown  in  Figure  11-71  before  input  to  the  new  adaptive  update 
algorithm.  In  all  cases,  the  adaptive -filter  length  is  15  points  per  channel.  At 
various  points  in  this  subsection,  four  events  (all  recorded  during  1972)  are 
used  to  simulate  mixed  events.  Table  111  — 1 lists  the  pertinent  data  fo  * these 
events.  Each  of  these  events  is  considerably  stronger  than  the  background 
noise  level  and  has  good  signal  similarity  across  ALPA.  As  in  Section  II,  the 
sites  at  ALPA  are  referenced  in  terms  of  their  transmission  order.  Figure 
1-1  and  Table  1-1  in  the  introductory  section  give  the  array  geometry  at  ALPA. 

To  simulate  the  first  mixed  event,  the  data  samples  for  the  two  Tonga 
Islands  events  are  summed.  The  event  from  day  29b  serves  as  the  on-azimuth 
signal,  while  the  evet  t from  day  335  serves  as  the  off-azimuth  interfering 
event.  Start  times  for  the  two  data  samples  were  chosen  so  thas  the  on-azimuth 
event  from  day  296  arrives  about  three  minutes  later  than  the  interfering  event 
from  day  335.  The  interfering -event  azimuth  was  shifted  from  206°  to  25°  to 
achieve  180°  azimuthal  separation  between  the  two  events  used  in  the  simula- 
tion. Likewise,  the  event  from  day  335  was  scaled  so  that  the  peak-to-peak 
amplitude  of  the  on-azimuth  event  from  day  296  is  18  dB  less  than  the  corres- 
ponding amplitude  of  the  interfering  event  at  the  single-sensor  level.  Fourteen 
sites  (all  but  sites  7,  16,  17,  18,  and  19)  constitute  the  channel  inputs  to  the 
beamformers.  In  this  first  mixed-event  simulation,  the  signal-to-noise  ratio 


III- 5 


I 

A 

Z 


remains  fixed  at  -18  dB  and  the  convergence  parameter  H for  the  adaptive- 
filter  update  varies  in  order  to  observe  the  effect  of  the  convergence  rate  on 
adaptive  processing  detection  gain  relative  to  beamsteering. 

Figure  III-2  displays  the  0.  04  Hz  beamsteer  response  corresponding 
to  the  205°  steer  azimuth  used  to  time-align  the  on-azimuth  signal  at  the  14 
sites  processed.  At  the  25°  azimuth  of  the  interfering  event,  the  beamsteer 
response  is  almost  21  dB  lower  than  at  the  205°  look  azimuth.  The  beamsteer 
output  traces  for  the  first  interfering-event  simulation  appear  in  Figure  III- 3 . 
The  upper  trace  is  the  beamsteer  output  for  the  data  sample  from  day  296  with 
the  on-azimuth  signal.  The  middle  trace  is  the  beamsteer  output  for  the  data 
sample  from  day  335  with  the  interfering  event  from  25°.  The  bottom  trace 
is  the  composite-sample  beamsteer  output  formed  by  summing  the  upper  and 
middle  traces.  Even  though  the  interfering  event  is  attenuated  by  approxi- 
mately 21  dB,  no  clear  indication  of  the  on-azimuth  signal  in  the  top  trace  is 
visible  on  the  composite-sample  beamsteer  output.  The  ratio  of  the  maximum 
peak-to-peak  amplitude  during  the  time  gate  where  the  on-azimuth  signal  is 
present  to  the  maximum  n»ak-to-peak  amplitude  before  its  arrival  is  a com- 
monly used  measure  for  the  reliability  of  a detection.  To  avoid  confusion  with 
the  term  signal-to-noise  ratio,  which  in  a strict  sense  would  apply  to  amplitude 
or  power  ratios  of  the  on-azimuth  signal  in  the  top  trace  to  the  interfering  event 
in  the  middle  trace,  the  term  amplitude  rise  is  used  in  referring  to  the  com- 
posite sample  beam  output.  For  the  composite-sample  beamsteer  trace  in 
Figure  III- 3 , the  amplitude  rise  is  -2  dB,  so  that  no  detection  is  possible. 

For  a convergence  parameter  value  of  0.  007,  Figure  III-4  shows  the 
ABF  response  at  0646  corresponding  to  the  beamsteer  response  in  Figure  111-2. 
At  the  25°  azimuth  of  the  interfering  event,  the  ABF  response  is  about  4 dB 
lower  than  the  beamsteer  response.  The  adaptive  beam  outputs  for  H=  0.007 
appear  in  Figure  III-5.  At  this  slow  convergence  rate,  the  ABF  traces  are 
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FIGURE  III- 2 

14-CHANNEL  BEAMSTEER  RESPONSE 
FOR  THE  SIMU  LATED  MIXED  EVENT 
FROM  DAYS  296  AND  335  OF  1972 
(INTERFERING  EVENT  AZIMUTH  25°1 
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ABF  BEAM  PATTERN  AT  296/06.46.00 
CONVERGENCE  RATE  15  0.0070 
BEAM  LOOK  VELOCITY  IS  3.6,  LOOK  AZIMUTH  205.0 
FREQUENCY  IS  0.04000  HZ,  PERIOD  25.0  SECONDS 

FIGURE  III  - 4 

14-CHANNEL  ABF  RESPONSE  AT  0646 
FOR  THE  SIMULATED  MIXED  EVENT 
FROM  DAYS  296  AND  335  OF  1972 
(INTERFERING  EVENT  AZIMUTH  25°) 
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similar  to  the  time-shift-and-sum  traces.  The  principal  differences  a slightly 
greater  attenuation  of  the  interfering  event  in  the  middle  trace.  As  before,  the 
composite-sample  beam  does  not  reflect  the  presence  of  the  on-azimuth  signal: 
the  amplitude  rise  on  the  composite-sample  trace  is  about  0 dB. 

As  the  convergence  rate  increases,  adaptive  beamforming  is  able  to 
suppress  the  interfering  event  more,  and  the  measured  amplitude  rise  on  the 
composite  trace  is  greater.  Figure  111-6  presents  the  adaptive  beams  at 
M=  0.  05.  The  amplitude  rise  on  the  composite  trace  is  approximately  4 dB, 

6 dB  higher  than  on  the  beamsteer  output.  The  greatest  amplitude  rise  occur- 
red at  a convergence  rate  of  0.2.  Figure  111-7  is  a plot  of  the  adaptive  filter 
outputs  at  this  rate.  The  amplitude  rise  on  the  composite  trace  is  about  6 dB. 
Even  though  there  is  some  signal  distortion  at  this  rapid  convergence  rate,  a 
detection  can  definitely  be  made.  In  this  case,  adaptive  processing  produced 
an  8 dB  detection  gain  over  beamsteering.  The  maximum  peak-to-peak  ampli- 
tude on  the  composite  trace  at  /i  = 0.  2 is  about  2 dB  lower  than  the  corres- 
ponding amplitude  for  the  on-azimuth  signal  in  the  top  trace.  This  difference 
is  equivalent  to  0.  1 magnitude  units. 

An  interesting  aspect  of  Figure  111-7  is  that  part  of  the  interfering-event 
waveform  in  the  middle  trace  closely  resembles  the  negative  of  the  on-azimuth 
signal  in  the  top  trace.  In  the  maximum  likelihood  adaptive  filter  update  algo- 
rithm, the  adaptive  processor  attempts  to  create  from  the  data  sample  con- 
taining the  off-azimuth  event  a filter  output  180°  out  of  phase  with  the  filter  out- 
put from  the  data  sample  containing  the  on-azimuth  event.  Even  though  the  on- 
azirruth  event  is  passed  virtually  undistorted,  the  addition  of  the  filtered  off- 
azimuih  event  tends  to  produce  mutual  cancellation.  This  phenomenon  is  con- 
sistent with  the  maximum  likelihood  filter  design  specifications  (minimization 
of  the  filter  output  power  subject  to  unity-response  constraints  in  the  beam  lo  jk 
direction).  Even  with  partial  cance  llation,  however,  adaptive  beamforming  sup- 
presses the  interfering  event  so  strongly  that  a detection  gain  of  8 dB  is  still 
achieved. 
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ABF  OUTPUT  FOR  THE  MIXED  EVENT  FROM  DA  fS  2?^  AND  335  OF  1972  (M  = 0.05, 
SIGNAL  AZIMUTH  205°,  INTERFERING  EVENT  AZIMUTH  25°,  INTERFERING 

EVENT  1 8dB  ABOVE  THE  ON-AZIMUTH  SIGNAL) 
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The  second  mixed-event  simulation  employs  the  Tonga  Islands  event 

from  day  296  for  both  the  on-azimuth  signal  and  the  off-azimuth  interfering 

event.  For  the  data  sample  with  an  on-azimuth  signal,  the  simulation  uses 

the  data  exactly  as  recorded.  The  apparent  arrival  azimuth  of  the  event  is 
o o 

shifted  from  205  to  25  to  form  the  data  sample  serving  as  the  interfering 
event.  Neither  data  sample  is  scaled,  so  that  both  events  are  equally  strong 
at  the  single-sensor  level.  Likewise,  the  start  times  for  both  samples  are 
identical,  so  that  both  events  arrive  simultaneously  at  site  8,  the  central  site 
at  ALP  A.  Fifteen  sites  (all  but  sites  7,  17,  18,  and  19)  are  input  to  the  beam- 
forming process.  Figure  III-8  displays  the  adaptive  filter  outputs  at  M = 0.5. 

The  beam  output  for  the  off-azimuth  event  is  almost  a perfect  replica  of  the 
on-azimuth  signal  except  for  a 180°  phase  shift.  The  result  on  the  composite 
sample  trace  is  a near-zero  output  due  to  mutual  cancallation  of  the  events. 

Just  to  demonstrate  that  near-perfect  cancellation  does  not  happen  at  all  con- 
vergence rates.  Figure  III- 9 shows  the  coi  responding  adaptive  beams  at 

0.002.  This  second  mixed-event  simulation  is  presented  for  illustrative 
purposes  only.  When  the  same  event  is  used  for  both  samples  in  simulating 
an  interf'  0 event  situation,  each  channel  is  the  sum  of  two  perfectly  cor- 
related but  time-shifted  time  series.  In  accordance  with  its  design  objective, 
the  maximum  likelihood  adaptive  beamformer  easily  reduces  the  composite 
sample  energy  to  a bare  minimum  at  high  convergence  rates.  The  perfect  cor- 
relation between  corresponding  channels  of  the  two  data  samples  utilized  makes 
this  simulation  extremely  unrealistic.  In  practice,  differences  in  event  wave- 
forms, dispersion  characteristics,  relative  arrival  times,  and  relative  strengths 
combine  to  prevent  such  near-perfect  cancellation  of  interfering  events. 

The  last  mixed-event  simulation  in  this  subsection  sums  the  Kamchatka 
event  from  day  276  and  the  New  Guinea  event  from  day  007.  The  Kamchatka 
event,  at  an  azimuth  of  273°  from  ALP  A,  serves  as  the  on-azimuth  signal, 
while  the  New  Guinea  event  serves  at  the  off-azimuth  interfering  event.  The 


(M=  0.  002,  SIGNAL  AZIMUTH  205°,  INTERFERING  EVENT  AZIMUTH  25  , 
INTERFERING  EVENT  EQUAL  IN  STRENGTH  TO  THE  ON-AZIMUTH  SIGNAL) 


azimuth  of  the  New  Guinea  event  was  shifted  from  2 53°  to  93°  to  obtain  180° 
separation  between  the  two  events.  With  the  start  times  selected,  the  inter- 
fering event  arrives  eleven  minutes  before  the  on-azimuth  signal.  Twelve 
sites  (all  but  3,  4,  5,  6,  7,  11,  and  18)  are  used  in  forming  the  beams.  The 
beamsteer  response  corresponding  to  a steer  azimuth  of  273°  appears  in  Fig- 
ure III- 10.  At  the  93  azimuth  of  the  interfering  event,  the  attenuation  is  15 
dB,  close  to  the  maximum  for  all  azimuths.  In  this  last  mixed-event  simula- 
tion, the  relative  strength  of  the  two  events  is  varied  in  order  to  observe  its 
effect  on  adaptive-processing  performance. 


Initially,  the  interfering  event  is  scaled  so  that  it  is  12  dB  above  the 
on-azimuth  signal  at  the  single-sensor  level.  Figure  II.I-ll  presents  the  beam- 
steex-  output  traces  for  this  case.  The  amplitude  rise  is  approximately  8 dB 
on  the  composite-sample  beam.  A visual  inspection  of  the  bottom  trace  in  Fig- 
ure HI- 1 1 confirms  that  a detection  is  possible  with  time-shift-and-sum  pro- 
cessing. Figure  III-12  displays  the  corresponding  adaptive  beam  outputs  for 
0.2.  The  amplitude  rise  on  the  composite  sample  trace  is  nearly  15  dB 
at  this  convergence  rate.  Thus  the  detection  gain  over  beamsteering  is  ap- 
proximately 7 dB.  The  drop  in  maximum  peak-to-peak  amplitude  from  the  top 
trace  to  the  bottom  trace  is  slightly  more  than  2 dB,  or  0.  1 magnitude  units. 


When  the  off-azimuth  event  is  scaled  so  that  it  is  18  dB  above  the  on- 


azimuth  signal  at  the  single-sensor  level,  no  detection  is  possible  on  the  com- 
posite sample  beamsteer  trace.  Figure  III- 13  pictures  the  beamsteer  traces 
in  this  situation.  The  amplitude  rise  on  the  composite-sample  beamsteer  trace 
is  only  2 dB.  Figure  III-14  is  a plot  of  the  corresponding  adaptive  beams  for 
M=  0.  007.  At  this  convergence  rate,  the  adaptive  filter  set  attenuates  the  in- 
terfering event  enough  to  produce  a 6 dB  amplitude  rise  on  the  composite  beam. 
As  the  convergence  rate  increases,  the  composite-trace  amplitude  rise  im- 
proves until  it  reaches  a maximum  of  about  12  dB  at  V = 0.2.  In  this  case, 
the  ABF  detection  gain  over  beamsteering  is  10  dB.  The  adaptive  beam  outputs 
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FIGURE  III-  10 


12-CHANNEL  BEAMSTEER  RESPONSE 
FOR  THE  SIMULATED  MIXED  EVENT 
FROM  DAYS  276  AND  007 
(INTERFERING  EVENT  AZIMUTH  93°) 
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BEAMSTEER  OUTPUT  FOR  THE  MIXED  EVENT  FROM  DAYS  276  AND  007 
(SIGNAL  AZIMUTH  273°,  INTERFERING  EVENT  AZIMUTH  93°,  INTERFERING  EVENT 

1 8dB  ABOVE  THE  ON-AZIMUTH  SIGNA  t ' 


for  this  convergence  rate  appear  in  Figure  III-15.  A clear  detection  is  possible 
in  the  bottom  trace.  Here  the  decrease  in  the  maximum  peak-to-peak  amplitude 
from  the  on-azimuth  signal  beam  to  the  composite  sample  beam  is  slightly  less 
than  2 dB.  The  detection  gains  at  high  convergence  rates  are  possible  in  the 
interfering-event  situation  chiefly  through  the  dramatic  suppression  of  the  off- 
azimuth  event  before  the  arrival  of  the  on-azimuth  signal.  Immediately  after 
the  arrival  of  the  signal,  the  adaptive  beamformer  is  unable,  at  least  for  a 
while,  to  effect  mutual  cancellation  of  the  interfering  events  when  the  events 
are  not  highly  correlated.  The  burst  of  energy  on  ihe  composite  sample  beam 
at  the  arrival  time  of  the  on-azimuth  signal  permits  the  signal  to  be  detected. 

At  an  off-azimuth  event  strength  24  dB  above  the  on-azimuth  signal, 
time-shift-and-sum  processing  is  even  less  effective  than  at  the  18  dB  level. 

For  purposes  of  comparison,  however,  Figure  III  16  presents  the  beamsteer 
outputs  for  the  24  dB  case.  The  corresponding  adaptive  beams  at  U = 0.  2 are 
shown  in  Figure  111-17.  On  the  composite  sample  beam,  the  amplitude  rise  is 
slightly  more  than  8 dB,  the  same  as  for  the  beamsteer  processor  at  an  off- 
azimuth  event  strength  only  12  dB  above  the  on-azimuth  signal.  Thus  adaptive 
beamforming  is  able  to  produce  detection  results  comparable  to  those  obtained 
from  beamsteering  at  a 12-dB-lower  event  separation.  In  the  24  dB  case,  the 
difference  in  maximum  peak-to-peak  amplitudes  between  the  top  trace  and  the 
composite  trace  is  negligible  for  the  ABF  outputs  at  H = 0.2. 

To  determine  the  maximum  event  separation  at  which  a detection  is  pos- 
sible with  adaptive  beamforming,  the  event  separatic  n is  increased  to  30  dB 
and  36  dB,  respectively,  in  the  last  two  cases  examined,  figures  III- 1 8 and 
HI-  1 ° disPlaY  the  ABF  outputs  at  a convergence  rate  of  0.  5 for  the  30  dB  and 
16  dB  event  separation  levels,  respectively.  At  a 30  dB  separation,  the  am- 
plitude rise  is  about  6 dB.  A borderline  detection  is  possible  here.  At  a 36 
dB  separation,  the  amplitude  rise  is  only  4 dB,  so  tha"  no  detection  could  be 
declared  without  additional  information.  An  example  of  such  information  could 
be  the  expected  arrival  time  determined  from  a short-period  event  bulletin. 
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An  examination  of  the  ABF  outputs  at  various  levels  of  separation  re- 
veals that  adaptive  beamforming  is  able  to  suppress  the  off-azimuth  interfer- 
ing event  by  about  33  dB  maximum  (1.65  magnitude  units)  before  the  arrival 
of  the  on-azimuth  signal.  Another  interesting  result  is  that  the  convergence 
rate  corresponding  to  the  greatest  amplitude  rise  increases  as  the  event  sep- 
aration level.  At  the  30  dB  and  36  dB  levels,  the  best  convergence  rate  is 
0.5,  the  same  rate  at  which  best  signal-to-noise  gain  is  achieved  in  the  pre- 
sence of  background  noise  alone. 

The  results  of  this  subsection  demonstrate  that  adaptive  processing 
produces  significant  detection  gains  over  conventional  time-shift-and-sum 
beamforming  for  mixed  long-period  events.  The  results  at  different  levels  of 
event  separation  for  the  interfering  event  simulated  by  summing  the  Kamcha- 
tka event  from  day  276  and  the  New  Guinea  event  from  day  007  indicate  a de- 
tection gain  over  beamsteering  of  more  than  12  dB. 

C.  ADAPTIVE  BEAMFORMING  APPLIED  TC  SHOP  T-PERIOD  INTER- 
FERING EVENTS  AT  KOREA 

Some  processing  of  simulated  interfering  events  has  been  performed 
with  data  from  the  Korean  short-period  array.  The  preliminary  results  de- 
pend to  a great  extent  on  the  frequency  response  of  the  prefilter  applied  to  the 
data  used  as  input  to  the  beamformers.  With  a l-to-3-Hz  bandpass  filter  on 
the  data  from  the  individual  sites  at  Korea,  adaptive  beamforming  produced 
some  detection  gain  over  beamsteering.  The  replacement  of  that  prefilter  with 
a prefilter  passing  energy  between  1 and  2 Hz  improved  the  detection  perform- 
ance of  adaptive  processing  relative  to  beamsteering.  To  what  extent  the 
improvement  is  due  to  the  lowered  frequency  band  or  to  the  reduced  bandwidth 
of  the  prefilter  is  unknown  at  this  time.  Until  it  is  possible  to  specify  an  ap- 
propriate frequency  response  for  the  prefilter  applied  to  the  individual  sensor 
outputs,  a presentation  of  results  is  not  likely  to  provide  a reliable  indication 
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of  the  potential  processing  gains  of  adaptive  filtering  over  time-shift-and-sum 
r ’•ocessing.  For  this  reason,  no  detailed  short-period  interfering-event  re- 
sults appear  in  this  report.  In  future  work,  a better  understanding  of  the  fac- 
tors influencing  short-period  adaptive -beamforming  performance  will  be  a 
major  topic  of  investigation. 

D.  SUMMARY 

The  results  of  Subsection  B demonstrate,  for  the  two  cases  where  mix- 
ed events  are  simulated  by  summing  two  distinct  data  samples,  the  adaptive 
beamforming  produces  significant  detection  gains  over  beamsteering  for  mix- 
ed long-period  events.  These  gains,  which  occur  at  high  convergence  rates, 
are  achieved  in  a roundabout  way.  Before  the  arrival  of  the  on-azimuth  signal, 
attenuation  of  the  off-azimuth  interfering  event  is  pronounced.  Immediate’/ 
after  the  signal  arrival,  adaptive  beamforming  is  unable,  at  least  for  a while, 
to  produce  mutual  cancellation  of  the  interfering  events.  The  burst  of  energy 
on  the  composite-sample  beam  a.  the  same  time  as  the  arrival  of  the  on- 
azimuth  signal  permits  a signal  detection. 

In  the  first  mixed-event  simulation,  where  the  Tonga  Islands  events 
from  days  296  and  335  are  summed,  the  beamsteer  response  at  a 25-second 
period  in  the  direction  of  the  interfering  event  is  -21  dB,  considerably  lower 
than  normal.  When  the  off-azimuth  event  is  scaled  so  that  it  is  18  dB  stronger 
than  the  on-azimuth  signal  at  the  single- sensor  level,  the  amplitude  rise  on 
the  composite  sample  adaptive  beam  at  H = 0.  2 is  6 dB  as  compared  with  -2 
dB  on  the  corresponding  time-shift-and-sum  beam. 

For  the  final  mixed  event,  simulated  by  adding  the  Kamchatka  event 
from  day  276  to  the  New  Guinea  event  from  day  007,  the  0.  04  Hz  beamsteer 
response  in  the  direction  of  the  off-azimutu  event  is  -15  dB.  At  an  off-azimuth 
event  strength  18  dB  above  the  on-azimuth  signal,  the  amplitude  rise  on  tht 


composite-sample  ABF  output  is  12  dB,  10  dB  higher  than  the  2 dB  rise  on  the 
corresponding  beamsteer  trace.  At  an  event  separation  level  of  12  dB,  the 
amplitude  rise  on  the  time-shift-and-sum  beam  is  8 dB,  and  the  on-azimuth 
signal  is  just  detectable.  At  an  event  separation  level  of  24  aB,  the  amplitude 
rise  on  the  composite-sample  adaptive  beam  at  M = 0.  2 is  also  8 dB.  Thus, 
in  the  final  mixed-event  simulation,  adaptive  multichannel  filtering  exhibits 
the  same  detection  capability  as  beamsteering  when  the  event  separation  is  12 
dB  greater. 


In  those  cases  where  detection  is  possible  with  adaptive  processing, 
the  maximum  peak-to-peak  amplitude  on  the  composite-sample  beam  provides 
a reasonably  accurate  estimate  of  the  on-azimuth  signal's  maximum  peak-to- 
peak  amplitude.  The  maximum  amplitude  on  the  composite-samnie  adaptive 
beam  is  never  much  more  than  2 dB  lower  than  that  of  the  on-azimuth  signal. 


The  results  of  this  section  show  definite  promise  in  lowering  detection 
thresholds  for  Rayleigh-wave  arrivals  when  off- azimuth  interfering  events  are 
present.  Time-shift-and-sum  processing  is  seldom  able  to  reduce  the  detec- 
tion threshold  more  than  0.6  to  0.  9 magnitude  units  (12  to  18  dB)  below  that  of 
a single  seismometer.  Adaptive  beamforming  may  further  reduce  the  detec- 
tion threshold  more  than  0.  6 magnitude  units  (12  dB)  below  the  beamsteer  de- 
tection threshold.  However,  these  preliminary  results  span  only  a limited 
range  of  the  situations  possible.  Before  a definitive  assessment  of  adaptive- 
filtering  gains  over  beamsteering  can  be  made  for  the  mixed-event  problem, 
many  additional  simulations  are  needed  to  determine  the  factors  influencing 
adaptive-beamforming  performance  in  the  interfering-event  situation.  In  par- 
ticular, the  azimuthal  separation,  he  waveform  similarity,  the  relative  ar- 
rival times,  and  the  relative  strengths  of  interfering  events  need  to  be  studied 
to  determine  their  effect  on  the  mutual-cancellation  phenomenon.  For  example, 
when  the  azimuthal  separation  is  small  enough  that  the  off-azimuth  events  falls 
within  the  beamsteer- response  main  lobe,  adaptive-nrocessing  performance 
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relative  to  beamsteering  should  improve.  In  addition,  adaptive  beamforming 
needs  to  be  tested  on  smaller  long-period  arrays  such  as  the  Iran  Long-Period 
Array. 

A maximum  likelihood  adaptive-filtering  algorithm  produced  all  of  the 
results  in  this  section.  The  design  goal  of  this  algorithm  is  to  minimize  the 
filter  output  power  subject  to  unity- response  constraints  in  the  beam  look  di- 
rection. Better  detection  performance  might  occur  with  a Wiener  filter,  where 
the  design  goal  is  to  minimize  the  mean  square  difference  between  the  filter 
output  and  the  on-azimuth  signal.  In  any  case,  the  detection  gains  illustrated 
in  this  section  furnish  a useful  baseline  from  which  further  attempts  at  im- 
provement can  be  made. 
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SECTION  IV 
CONCLUSIONS 


For  unmixed  long-period  seismic  events  in  the  presence  of  background 
noise,  the  new  adaptive  algorithm 
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produces  much  higher  signal-to-noise  ratio  improvement  than  the  old  adaptive 


algorithm 
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used  in  the  previous  adaptive-beamforming  study.  The  improved  results  with 
the  new  adaptive  algorithm  are  due  to  a slowing  of  the  effective  adaptation  rate 
in  the  presence  of  on-azimuth  signals.  For  events  with  good  signal  similarity 
between  sites,  signal  degradation  with  the  new  adaptive  algorithm  is  dramat- 
ically reduced  at  high  convergence  rates.  Since  noise  reduction  increases 
significantly  as  a function  of  convergence  rate,  the  best  signal-to-noise  gains 
occur  at  high  convergence  rates  with  the  new  adaptive  algorithm.  In  contrast, 
severe  signal  degradation  at  high  convergence  rates  causes  the  maximum 
signal-to-noise  gains  for  the  old  adaptive  algorithm  to  occur  at  low  convergence 
rates  near  4=0.  005. 

When  the  full  ALPA  array  is  used  as  input  to  the  beamformers  instead 
of  a closely-spaced  six-channel  array,  the  measured  noise  reduction  relative 
to  beamsteering  is  lower  at  all  convergence  rates.  Differences  in  measured 
signal  degradation,  however,  are  highly  variable.  At  rapid  convergence  rates, 
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signal  degradation  is  2 dB  higher  for  one  event  and  2 dB  lower  for  another 
event.  In  most  cases,  the  estimated  signal- to-noise  gain  relative  to  beam- 
steering is  less  with  the  full  array  because  of  the  decreased  noise  reduction. 

The  new  signal-enhancement  prefilter  shown  in  Figure  11-71  is  a defi- 
nite improvement  over  the  old  prefilter  pictured  in  Figure  II- 1.  Not  only  is 
detection  gain  achieved  by  the  frequency  response  of  the  new  prefilter,  but 
also  the  new  prefilter's  emphasis  on  the  more  coherent  energy  at  low  frequen- 
cies increases  ABF  noise  reduction  relative  to  beamsteering  for  the  full  array 
by  more  than  3 dB  compared  with  the  old  prefilter.  This  additional  noise  re- 
duction occurs  with  both  31 -point  and  15-point  adaptive -filter  lengths.  The 
difference  in  noise  reduction  is  only  0.  3-0.4  dB  for  the  two  adaptive -filter 
lengths.  For  two  events  with  good  signal  similarity,  the  difference  in  signal 
degradation  is  less  than  0.1  dB,  so  that  the  loss  in  estimated  signal-to-noise 
gain  is  0.  3-0.4  dB  with  the  shorter  filter  length.  On  the  other  hand,  signal 
degradation  for  one  event  with  poor  signal  similarity  is  more  than  2.  1 dB  low- 
er with  the  15-point-long  adaptive  filter,  so  that  the  estimated  signal-to-noise 
gain  is  actually  1.8  dB  higher.  The  significantly  diminished  computational 
load  probably  justifies  the  use  of  the  shorter  adaptive  filter  length  in  any  on- 
line implementation  of  adaptive  beamforming. 

When  vertical-component  data  from  the  full  ALPA  array  are  passed 
through  the  new  prefilter,  the  signal-to-noise  gain  over  beamsteering  falls 
within  a 3-to-8  dB  range  for  an  adaptive  filter  length  of  15  points  per  channel. 
Normally,  a 6 dB  gain  might  be  expected  in  this  operational  mode.  Thus  adap- 
tive beamforming  could  conceivably  lower  the  Rayleigh  wave  detection  threshold 
at  ALPA  by  0.  3 magnitude  units. 

Off-azimuth  events  are  strongly  suppressed  by  adaptive  beamforming 
at  = 0.5,  the  same  convergence  rate  where  the  maximum  gain  over  beam- 
steering is  achieved  for  on-azimuth  events.  Some  off-azimuth  events  are  al- 
most totally  annihilated.  Ordinarily,  the  strongest  off-azimuth  events  are 
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attenuated  the  most.  For  some  signals  30°  to  60°  from  the  steer  direction, 
however,  considerable  energy  leaks  into  the  adaptive  beams,  even  for  the 
strongest  events.  The  new  prefilter  reduces  off-azimuth  event  suppression 
somewhat.  Even  with  the  new  prefilter,  however,  the  attenuation  of  off- 
azimuth  events  is  considerable.  Thus  the  beam-narrowing  capability  of  adap- 
tive processing  can  provide  greater  directivity  at  ALPA  and  other  similar  long- 
period  arrays. 

Until  now,  this  section  has  discussed  adaptive-processing  results  for 
single  seismic  events  in  the  presence  of  background  noise.  The  remainder  of 
this  section  deals  with  mixed-event  simulations  where  two  data  samples,  each 
containing  a signal,  are  summed  to  create  a composite  sample  containing  an 
interfering  event.  The  mixed-event  results  in  this  report  are  preliminary  and 
extend  over  only  a limited  range  of  the  situations  possible. 

For  the  two  mixed-event  simulations  where  two  different  data  samples 
are  added,  adaptive  beamforming  achieves  significant  detection  gains  over 
beamsteering.  These  gains,  which  are  realized  at  high  convergence  rates,  oc- 
cur in  a peculiar  way.  Before  the  on-azimuth  signal  arrival,  the  off-azimuth 
interfering  event  is  strongly  attenuated.  After  the  signal  arrival,  adaptive 
beamforming  cannot  produce  immediate  mutual  interfering  - event  cancellation. 
The  energy  burst  coinciding  with  the  on-azimuth  signal  arrival  permits  a sig- 
nal detection. 

In  the  first  mixed  event  sin  ulation,  the  time-shift-and- sum  beam  pat- 
tern in  the  off-a^imuth  event  direction  is  much  lower  than  normal  (-21  dB). 

The  off-azimuth  event  is  18  dB  stronger  than  the  on-azimuth  signal  at  the  sin- 
gle sensor  level.  The  beamsteer  amplitude  rise  after  signal  onset  is  -2  dB. 

For  ihe  adaptive  beam  ac  a 0.  2 convergence  rate,  the  amplitude  rise  is  6 dB. 
This  8 dB  difference  occurs  despite  the  excellent  benmsteer  response  toward 
the  off-azimuth  evert. 
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The  beame.teer  response  toward  the  off-azimuth  event  is  -15  dB  in  the 
final  mixed-event  simulation.  At  a 12  dB  event  separation  level,  the  beam- 
steer  amplitude  rise  is  8 dB,  and  the  on-azimuth  signal  is  just  detectable.  At 
a 24  dB  event  separation  level,  the  amplitude  rise  on  the  composite-sample 
adaptive  beam  at  M = 0.  2 is  also  8 dB.  Thus  adaptive  beamforming  produces 
the  same  detection  results  as  beamsteering  when  the  event  separation  is  12  dB 
greater. 

In  those  cases  where  detections  are  made  with  adaptive  beamforming, 
the  maximum  oeak-to-peak  amplitude  on  the  composite-sample  beam  furnishes 
a good  estimate  of  the  on-azimuth  signal’s  maximum  peak-to-peak  amplitude. 
The  maximum  amplitude  on  the  composite-sample  adaptive  beam  is  never  much 
more  than  2 dB  lower  than  that  of  the  on-azimuth  signal. 

The  mixed-event  results  show  definite  promise  in  lowering  detection 
thresholds  for  Rayleigh-wave  arrivals  masked  by  off-azimuth  events.  Beam- 
steering seldom  reduces  the  detection  threshold  more  than  0.6  to  0.  9 magni- 
tude units  (12  to  18  dB)  below  that  of  a single  seismometer.  Adaptive  beam- 
forming may  further  reduce  the  detection  threshold  more  than  9.6  magnitude 
units  (12  dB)  below  the  beamsteer  detection  threshold.  However,  these  re- 
sults are  quite  limited  in  scope.  Many  additional  simulations  are  needed  to 
determine  the  factors  affecting  adaptive-processing  performance  in  the 
mixed-event  situation.  Furthermore,  adaptive  beamforming  needs  to  be  tried 
at  smaller  long-period  arrays  such  as  the  Iran  Long-Period  Array. 

A Wiener  filter,  where  the  design  goal  is  to  minimize  the  mean  square 
difference  between  the  filter  output  and  the  on-azimuth  signal,  might  provide 
better  detection  performance  than  the  maximum  likelihood  adaptive  algorithm, 
which  sometimes  produces  mutual  cancellation  of  the  interfering  events. 
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